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Preface 


Business Mathematics is a branch of applied mathematics that uses cal- 
culus, algebra and mathematical programming, to mention but a few. It finds 
applications in finance and economics. The mathematical analysis of stock 
markets, portfolios and the management of risks may also be covered by this 
subject. 

I began working on this book in October 2005 for the course Business 
Mathematics that I taught at Mahidol University at Kancanaburi. This vol- 
ume is the collection of hand-outs which were given to students before each 
lecture. Apart from examples and exercises, exam papers are also included. 
They are put in the appendix and represent in a way more examples. 

These hand-outs were used hand in hand with another set of printed 
pages whose contents are in larger letter, which were used on the camera 
projector. It would be good if in the future we could compile these hand-outs 
into lecture notes to be published as a book. 

There were 66 students in the class. From after the midterm exam my 
main concern has been for the students to do the exercises and problems by 
themselves. Since in practice this was hardly the case with the homework, we 
turned to holding within our class sessions of exercise practice and quiz. This 
has worked to a certain degree, all the time one main problem remain being 
the large size of the class, I do hope the students have learnt the importance 
of problem-oriented approach, and that each of them will find it useful in the 
future. I hope they have learnt from me as much as I have learnt for and from 
them. 


Kit Tyabandha, PhD 
Bangkok, 14" January, 2007 
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Graph and derivative 
25%" October 2005 


Definition 1. Let n be a positive integer. Then x” means that x is multi- 
plied to itself n times. 


§ 


Definition 2. A monomial is an expression cosisting a real-numbered co- 
efficient times one or more variables each raised to the power of a positive 
integer. Adding and subtracting monomials to one another give us a polyno- 
mial. A univariate polynomial is a polynomial of one variable, of the form 
Ant" +...+ a,x + a9, and its degree is the highest power of that variable, 
that is n. The degree of a polynomial is sometimes known as its order. A 
polynomial can be simplified into a product of two polynomials through the 
process called factoring. A polynomial equation is an equation of the form 
p(-) = 0, where p(-) is a polynomial. 

8 


Definition 3. A quadratic equation is a second-degree polynomial equation, 
that is one of the form az? + bz +c =0, where a, b and c are constants and 


a#0. 
8 


Example 1. The factors of mx? + nz +p are ax + b and cx + d, where 
ac = m, bd = p and ad+ bc = n. Quadratic equations may be solved by 
factoring. 


Example 2. A quadratic polynomial of the form x? +bz can be transformed 
into the form (x +a)” by adding (6/2)? and then factor the result. Completing 
the square in the quadratic equation ax? + bx +c = 0 leads one to 


fase 2 4a 
7 oa} ~ 4a? 


which gives us the quadratic formula 


—b+ Vv b2 — dac 


2a 


c= 


Definition 4. At a point on the graph where the first derivative is zero, if 
the second derivative is positive, then that point is a local minimum; if the 
second derivative is negative, it is a local maximum; on the other hand, if the 
second derivative is either zero or is undefined, it is an inflection point. An 
asymptote is a straight line to which a non-linear curve smoothly approach as 
it goes towards infinity, never reaching it. At any point on a graph, if the first 
derivative is positive the function at that point is increasing, and if negative 
it is decreasing. At any point on a graph, if the second derivative is positive 
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it is said that the function is convex at that point, and if negative the latter 
is said to be concave. 


§ 


global maximum 
elas 


or 


global minimum 


(a) (b) 


Figure 1 (a) global maximum-, and (b) global minimum points of a 
function. 


y’=0 
y’<0 


local maximum 


ee. 


local minimum 
y’=0 
y’>0 


Figure 2 Local minimum- and local maximum points of a function. 
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concave up 


curves down 


y’>0 


y’<0 


concave down 


curves up 


(a) (b) 
Figure 3 The two curvature types, namely (a) concave up (y" > 0) and, 
(b) concave down (y" < 0). 


i 


Figure 4 The four possible curvatures in two dimensions, considering 
both y' and y", (a) y" >0 and, (b) y" <0. 


Figure 5 Inflection points, where y” = 0, (a) with y” ee and, 
(b) with y” decreasing. 


y’=0 y =0 


y’’<0 
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(a) (b) 


Figure 6 Stationary inflection points, where both y' = 0 andy" = 0, (a) 
with y" increasing and, (b) with y" decreasing. 


Example 3. The cubic function of the form y = ax? + bz? + cx +d has 
local maximum and minimum, when these exist, at the point where the first 
derivative is zero, that is when 3ax? + 2bx +c =0. At such points, 


—bt Vb? — 3ac 


3a 


Depending on whether b? — 3ac is positive or negative, the local minimum and 
maximum either exist or do not exist. As examples, when a=1,b=2,c=3 
and d = 4 this value is —5, and therefore the value of « becomes complex. In 
the space of real number this means that there is neither a maximum- or a 
minimum point. Instead, in this case we have an inflection point, which is the 
point where the second derivative becomes zero, that is to say, 6ax + 2b = 0, 
or @ = —z£. Figure 7 shows a plot of this case. 


y = ax? + bx? + ex +d, where a=1,b=2,c=3 andd=4 
10 T T 


8 EI 


6 a 


4h 4 


~10 L Hl Ll L L L H 


Figure 7 The cubic function y = x? + 227 +3244. 


Figure 8 shows the case wherea= 1, b=5,c=4 andd=3 
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y = ax? + ba? + cx +d, where a=1,b=5,c=4 andd=3 
16 T T T T T 


5h 4 


Figure 8 The cubic function y = x? + 5a? + 4a +3. 


Example 4. The hyperbolic function of the form (# — a)(y — 6) = c has 
as its asymptotes the lines x = a and y = b. Figure 9 shows the graph when 
a=1,b=2andc=3. Here the asymptotes are x = 1 and y = 2. 


(« —a)(y —b) =c, where a= 1, b =2 and c=3 
10 T T T T 74 


Bh 


| 2 


Figure 9 The hyperbolic function (x — 1)(y — 2) =3. 


Definition 5. The straight line has a form ax + by +c=0, or y=mar+d 
where m = —a/b and d = —c/b. The y-intercept is the point where x = 0, 
that is y = —c/b or y = d. The az-intercept is the point where y = 0, that is 
x = —c/a or = —d/m. 


§ 


Example 5. Given a sum of money, an isocost line represents the different 
combinations of two inputs of production that can be bought. The general 
formula is p,k + pil = e, where k and | are capital and labour, p, and ~ their 
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respective price, and e the allotted expenditures. In other words, pzk + pil — 
e = 0, or k = e/px — (pi /p,)t. All the values concerned are positive, which 
means that the graph we are interested in is in the first quadrant only. 


Example 6. At the simplest case where there are only two activities to 
choose from, if x and y are the two activites, then t,p, + typy = m, where tz 
and ty the number of hours spent on activities 2 and y respectively, pz and 
py are the price per unit hour of x and y, and m is the limit of income. 


Definition 6. The general demand function is of the form 


da = £(p, Y, Ps, De, ta, a,.- ) 


where qq is the quantity demand of good 2, p the price of x, y the income of 
the consumer, p, the price of substitute goods, p,. the price of complementary 
goods, t, the taste or fashion of the consumer, and a the advertisement level. 
In its simplest case where all other factors are constant, the demand equation 
takes the form p = c, — c2qa, where p is the price-, while qq the quantity 
demanded of the good z, and c; and cz are positive constants. The general 
supply function is of the form q, = f(p, c, po, te, n,0,---), where q, is the quan- 
tity supplied of good z, p the price of x, c the cost of production, po the price 
of other goods, ¢, the available technology, n the number of producers in the 
market, and o other factors, for example tax and subsidies. The simplified 
relation for the supply is p = c, + cogs, where q, is the quantity of x supplied, 
and c, and c2 are positive constants. 


§ 


Example 7. When a tax of ¢ per unit is imposed, the supply function 
becomes p — t = c, + coq, where c, and cp are positive constants, and the 
total cost function becomes c; = ce + (k+t)q, where k is the cost of producing 
each unit. Here also cy = kq. 


Example 8. Revenue is the amount of money received when a firm sells its 
output. The relation is r, = pq, where r; is the total revenue, p the price and 
q the quantity. 


Definition 7. Elasticity « of x with regard to y means the ratio of the 
change in x to the change in y. Here y could be some economic variable, for 


example price or income. Hence, the price elasticity of demand is €g = aa 2 
the price elasticity of supply €, = aie - The point elasticity is the elasticity 


value calculated at a point, the arc elasticity or the midpoint elasticity is 
the same averaged over an interval. The arc price elasticity of demand or 
supply is then €g = M4 Pit p2 | The arc income elasticity of demand or supply 


: P +492 
; — dg yitye 
18 Ey = dy qitga” 


§ 


Elasticity measures the sensitivity or responsiveness of a certain quantity 
to changes in some variable. 
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Definition 8. A function of one independent variable is a relation in the 
form y = f(x) such that there exists one and only one value of y in the range 
of f for each real number x in the domain of f. The variable y is called the 
dependent variable. 


§ 


Definition 9. A function is called multivalued function if the opposite to 
Definition 8 is true, that is there exist more than one values of y for some z. 


§ 


Definition 10. An implicit function is a function in which both dependent- 
and independent variables appear on the same side. An explicit function is one 
where the dependent variable is on the left hand side-, and the independent 
variable on the right hand side of the equation. 


§ 


Example 9. Examples of explicit functions are y = 52, y = 2° +27 —7 
and y = e* + (c + 1)Inz. Examples of implicit functions are x + y = 1, 


x? + 3ay-yty= eT and x7 Iny = e¥(x +2°). 


Definition 11. The exponential function has the general form y = a®, where 
the base a is a constant and z is called the indez-, power-, or the exponent of the 
exponential function. The logarithmic function is then log, y = x, providing 
that both a and y are positive real numbers. We call a” the 2" power of a, 
and call log, y the logarithm of y to the base a. When the logarithmic base 
is the natural number, we write log, n = Inn. 

8 


The logarithmic function is the inverse of the exponential function and 
vice versa. 


Theorem 1. Let p and gq be real numbers, a and b positive numbers, and 
m and n positive integers. Then, 


a? - al = aPt4 


1 

a?P=— 

aP 
(ab)? = aPbP 

1. 

Va=ar 
vy ail” Sa: 
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[a Ya 


b 


§ 
Theorem 2. 
log, mn = log, m+log,n 
log, “= log, m — log, n 
n 
log, m? = plog,m 
I 
log, a 
§ 


Both Definition 12 and Definition 13 are a definition of limits in calculus. 
Definition 13 is a more formal and definitive one. 


Definition 12. If f(x) is a function which draws closer to a unique finite 
real number J for all values of x as the latter draws closer to a, but x # a, 
then / is called the limit of f(x) as x approaches a. In notation this is, 


lim f(#) =1 
na 


§ 


Definition 13. For a function f(x), lim,_,, f(x) =! if and only if for every 
€ > 0, there exists 6 > 0 such that |f(x) —1| < € whenever 0 < |x —al <6. 


By Definition 13 we mean that one can get f(x) to be as close to J as one 
wish, provided that x gets close enough to a. 


Theorem 3. Providing that both lim,_,, f(z) and limz_,, g(z) exist, and 
let k be a constant and n a positive integer, then the rules of limits are the 
following. 

lmk=k 


za 


lim x” =a” 
“za 


lim kf(x) =k lim f(x) 


Jim [f(z) + g(@)] = lim f(a) + lim g(x) 
lim [f(2) - g()] = linn f(a) - Timm g(a) 
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, provided that lim g(x) 4 0 
za 
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n 
lim [f(2)]” = [lim f(2)] , provided that n > 0 
za xzra 


§ 


Definition 14. Let y = f(z). Then, the derivative of y with respect to z is, 


dy _.. f(a +h) — f(z) 
Gt h 


The various notations for the derivative include df(x)dz, as f' (x), y’, Dy 
and D(f(a)). The process for obtaining this is called differentiation. 


§ 


Theorem 4. Let y = x”. Then the power rule for differentiation is the 
following. 

dy = n-1 

eo 


§ 


Theorem 5. The derivatives of exponential and logarithmic functions are 
respectively 42 = e® and “22 = 1. 
§ 
It turns out that e” is the only function which is changed by neither dif- 
ferentiation nor integration. Leonhard Euler (1707-83) found that the natural 
number € = limp soo (1 + i”) = 2.71828 .... 


Exercise 1. Find the first derivative of the following: 


1. y = 2e” 2. y = 2In(z) 3. Cy = A +Ine 

4. p= 100(1 — e*) 5.qg=VI4+4 6. p=1+1lngq 
n 3t 

7.q= 3-H 8. ca =Ing+} 9. p= 98.55 


10.c=2?-—e74+Inz 11. c= at +0.8le% 12. + = 34.3e% 


Theorem 6. Let y be a function of u, and u a function of x. Then the 
chain rule states that, 

dy _ dy du 

dx dudz 


Exercise 2. Differentiate the following: 
l.p=5+ 2. y = soar 

3 4. y = (82 — 5)" 

5. ¢ = 0.9572)? + 31 6. = fg +42 

7. y =(2—-0.72)-8 8. ¢ = 100 + 30e°8Y 

9. y=lMe tbh + oH —In(3a? -22+5) 10. p= a + 21n(q) 

1 1 


1. s = 123(1 — 3e~°-628) i 


2. Y= 
(52+6) 3 
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Theorem 7. If y = u(zx)v(x), then 


Similarly, if y = u(z)v(#) w(x), then 


dy —s-_ dw a we! gga 
dx de dx dx 
8 
Exercise 3. Find the first derivatives of the following: 
1. c=3q/24+¢ 2. p= Bare" 
8. y= 23/4 4. y= (t)e*% +5 
5. s = 10ye°'l9 6. cy = 7 —qlng+V@ 
7.c2=yln(y?)+/3yny 8.I= 5 In(3y + 2) 
9. p= q’e% 3 In(q? — 8) 10. ¢ = a +q/q—1—q+100 
11. y =#e' +10Int 12, 5 = olny 
Theorem 8. Ify= Or then 
dy _ ae ~ “ae 
da v? 
8 


Exercise 4. Find the first derivative of the following and simplify your 
answer: 


_ t1 = it = 
l.y=H Stes er ey, 
-1.6 
4,y= net 5. p= 100-4 6. e= 500 (1- ) 
a 7 2,6—4 — yt-l 
7 p=aa 8 = Gig 8 YR 


Definition 15. The total cost c; comprises a fixed cost and a variable cost, 
that is c; = cf + Cy. The total revenue r; is price times output, r; = pq. 


§ 


Definition 16. The marginal cost cm is the change in total cost caused by 
the production of an additional unit. The marginal revenue rm is the change 
in total revenue coming from the sale of an extra good. In other words, 


= ott and rm, = oe, where q is the output. 

8 
Definition 17. The average cost c, is the total cost per unit output, cg = a? 
The average revenue rq is the total revenue per unit output, rg = ve 

8 
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Example 10. From Definition 15, ry = pg, and from Definition 17, ry = 
raq, therefore p— ra. From Definition 17, cy = re and from Definition 15, 
Ca = Cf + Cy. Therefore cg = Caf + Cav, where the average fixed cost cas = Sf 


q 
and the average variable cost cay = ae 


Procedure 1 calculates and draws graphs of the total-, average-, and 
marginal costs, which are respectively c, cg and Cy. Similar graphs of other 
quantities may be drawn likewise, for example ones for the total-, average-, 
and marginal products. 


Procedure 1 Total-, average-, and marginal graphs 


Given: c;(z) 
Ca 
Cm — ¢} 
for each f € {ct,Ca,Cm} do 
find the critical values x, for f’ =0 
find f" 
n+ |X¢| 
for i=1 ton do 
if f"(z¢) > 0 then 
f(x) is convex and is the relative minimum of f 
elseif f" (x) < 0 then 
f(x§) is concave and is the relative maximum of f 
endif 
endfor 
find inflection points x» from f” =0 
endfor 
plot c¢(z), then cg(x) and cm(x) 


§ 


Problem 1. Let the total cost function be c = gq? — 40q? + 800. Sketch the 
graph to show the relationship between total-, average-, and marginal costs. 
If c¢ = f(q), then Cm = de /dgq, and if r; = f(g), then rm = dr¢/dq, where c is 
the total cost, r; the total revenue, and q the level of output. 


Problem 2. The demand function for a monopolist is gq = 100 — 4p. Find 
Tt, Tm and fq, then evaluate these at g = 10 and explain the results. Further, 
find g when r, = 0, then find rm, when rz = 0, and then explain whether 
one should sell at this value of g. Draw graphs of r;, rm and ra on the same 
diagram. 


Problem 3. The demand function for a good is p = 200 — q''8. Find rz, 
Tm and rg, and then compare the slopes of the r,- and rg curves. Evaluate 
Tt, Tm and rq at g = 11 and at q = 26, and give an explanation of the values 
obtained. Find g when r,, = 0, and find gq when rz = 0. When does the sales 
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of more units begin to reduce the total revenue? Draw graphs of r;, r,, and 
rq on the same diagram. 


Problem 4. The fixed costs of a firm are 2000 and the variable costs are 
4.7q. Write the equation for c; and evaluate the c¢ at q = 35. Write the 
equation for c,, and evaluate the c,, at q = 35. Explain the meaning of c,, 
for this function. 


Problem 5. The average cost function of a firm is given by c, = q? —10q+ 
20 + 66. Write expressions for c¢ and c,,, and then calculate c; at q = 17. 
Write the equations for cr and Cy. 


Problem 6. The average cost function of a firm is c, = 97+*+. Differentiate 
Cq with respect to q, then describe and explain how the former changes with 
the change in the latter. 


Definition 18. The relationship between input and output is called a pro- 
duction function, q = f(l,k,r,te,s,e,-.-), where / is labour, k phical capital 
such as buildings and machines, r raw materials, t. technology, s land, and 
e enterprise. Assuming a short period of time, then / becomes the only in- 
dependent variable and the other remaining factors are parameters, that is 


fixed, and therefore g = f(1). Then the marginal product of labour is pum = a 
and the average product of labour is pia = +. 
8 
The value pj, in Definition 18 is a measure of productivity. 
Definition 19. The marginal propensity to consume is pem = oo The 


marginal propensity to save iS Dsm = aa The average propensity to consume 
1S Dea = oe The average propensity to save is Dsa = 3 Here y is the income, 
c the consumption, and s the saving. 


§ 


Example 11. Since y = c+ 5, it follows that pea + Psq = 1. Further, 


dy _ d d = 
io et ay) and hence 1 = pom + Psm: 


Example 12. The consumption function may be described by c = co + by, 
where cp and 0 are positive constants. Then peg > Pem and Psm > Dsa- 


Definition 20. The profit is 7 = ry — Gq. At the break-even point x = 0, 
that is 7, = G. 


§ 


Problem 7. Find the elasticity of demand when gq = i where a and c are 
constants. 


Problem 8. The demand function for train journeys is q = oe where gq is 


the number of passengers in thousands. Find the elasticity of demand in this 
case. Determine the percentage change in demand when the fare increases or 
decreases by 10 per cent. 
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Problem 9. The demand for mineral water is given by g = 100 — 2p, where 
p is the price per bottle and gq the number of bottles demanded. Find the 
expressions for rz, 7m and rz. Find the quantity and price that maximise the 
revenue. Find an equation for the price elasticity of demand firstly in terms 
of p, and then again in terms of g. Show that r; is maximum and r,, is zero 
when €g = —1, and find q at this éq. 


Problem 10. A demand for a wine is given by p = 1200e~°-°7, where 
p is the price per bottle and gq the number of bottles demanded. Write the 
equations for r;, rm and rg. Find the price and quantity when the revenue is 
maximised. Show that ¢g = —1 when r; is maximised. 


Problem 11. Derive the equation 


1 


Problem 12. Show that the profit becomes maximised when 


p=7— 
Oa) 


Problem 13. A shop selling shirts has a demand function p = 300 — 10q 
and a total cost function ¢ = 150+ 8q. Find the equations for r; and 7. 
Find the number of shirts which must be sold firstly in order to maximise 
the profit, and then again to maximise the total revenue. Show that rm = re 
when the profit is maximised. 


Problem 14. The average cost function of a mobile phone is cg = 10+ 2, 


and the average revenue function for the same is rg = 35. Find r:, ¢, rm 
and c,,. How many mobile phones must be made and sold in order to break- 
even? Find the profit function, and show that neither profit nor revenue 
has a maximum. Explain using r,, and c», why there is no maximum. Plot 
the graphs for r;, c; and m7 on one diagram, and plot r,, and c on another 
diagram. Comment on the graphs. 
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Calculus of multivariable functions 
1%* November 2005 


Definition 21 talks about functions of n independent variables, and Def- 
inition 22 the partial derivatives of these. In Theorem 9 we have the product 
rule, in Theorem 10 the quotient rule, and in Theorem 11 the generalised 
power function rule. 


Definition 21. A function y = f(21,...,2,) is called a function of n inde- 
pendent variables if there exists one and only one value of y in the range of f 
for each tuple of real number (21,...,%,) in the domain of f. Here y is called 
the dependent variable while x;, i= 1,...,n, the independent variables. 


§ 


The word tuple in Definition 21 means an ordered list. It is also known 
as n-tuple, where n is the size of the list. 


Definition 22. Let a multivariable function be y = f(#1,...,2n). The 
partial derivative of y with respect to z;, where 1 <i <n, is a measure of the 
instantaneous rate of change of y with respect to 2; while x; is held constant 
for all 7 Ai, where 1 <j <n. This partial derivative is defined as 


Oy lim f(...,a; + Ag;,...) —f(a1,..-,¢n) 
Ox;  Ax;—>0 Aa; 


and can be written in either one of the following forms. 


Oy Of 
—, fa, (@1,---,; Un), fe;, OV Ya; 
Fs Goes Bos 1s Bn) favs OF 
§ 
Theorem 9. Let z= g(z,y)-h(z,y). Then, 
Oz Oh Og 
SEs gt er Ae 
dx © Oa Ox 
and 
Oz Oh Og 
—~_ = g er h eras 
oy oy oy 
§ 
Theorem 10. Let z= oe and h(z,y) £0. Then, 
dz _h- ge aes gh 
Ox h? 
end a adh 
Og: lepee ea 
Oy h2 
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8 
Theorem 11. Let z= [g(x,y)]”. Then, 
Oz n_-y OF 
ac 
and ; é 
<signet 28 
dy "ay 
8 
Exercise 5. Find the first-order partial derivatives of the following: 
1. c= 100(1+ e714) 2. q = 5I°-7k93 — A — 3k +10 
3. q=In3a+2eIny 4.2=83(1+y+y’) 
5. p = 150e° 4 6.qg=Inx+lIny 
zeae t+art+etQeytay? 8. 2=27y? 
9. q = 219-77 KO4 10. u= 82? y? 
11. z=ay+% 12. u= 7k 
8 


In Definition 23 we find the meaning of second-order partial derivatives. 
Theorem 12 is about critical points, and Procedure 2 is a procedure for de- 
termining critical points. 


Definition 23. Let z = f(x,y). Then, the second-order direct partial deriva- 
tives are 

ME: gs Se eee 

Ox \ Ox Oy \ Oy 


Oo? 6) 
fox, (fn)a; xa and respectively fyy, (fy)y, oe 


These are also written 


The cross partial derivatives are 


These are also written as 


072 


2 
fry, (fay, seis and respectively fyz, (fy)a, dndy 


Exercise 6. Find the second-order partial derivatives of the following: 
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. z = 19-29-38 _ \(10 — 1 — 3k) 
. c= 100(1 + e149) 

b. q=Inz+lIny g = 51°-4449-2 + 50 — 71 —5k 
7. gq = 919-77 49-4 : q=Ind5¢+alny 

9. p=110e8* = =10. z=ay+ £ 
i1.u=10rty* = 12. z= 23(1+y+y?*) 


1.z=a2+2y 
32% = 274? 


Co AW 


§ 


When the second derivative is negative, the curve is concave towards the 
origin. 


Theorem 12. For a multivariable function z = f(x,y) to be a relative 
maximum at (a,b) necessarily f,,f, = 0, and frz,fyy <0 and fre -fyy > (fey)? 
at that point. For the same at the same to be a relative minimum, necessarily 
f,,f, = 0, and fre, fyy > 0 and fre -fyy > (fey)? there. Moreover, an inflection 
point is a point (a,b) at which fro - fyy < (fey)?, and both fy, and f,, have 
the same sign. On the other hand, a saddle point is a point (a,b) at which 
fox -fyy < (fey)?, but fr and fy, are of different signs. 

8 


Procedure 2 Procedure for determining a critical point of a function 
with two independent variables 


Given z = f(x,y) and a point (a,b), at this point, 
if f, = 0 and fy =0 then 
(a, b) is a critical point 
if fz -fyy > (fey)? then 
if fp, <0 and f,, <0 then 
(a,b) is a relative maximum of z 
elseif f,, > 0 and f,, > 0 then 
(a,b) is a relative minimum of z 
else { 
endif 
elseif fry -fyy < (fry)? then 
if fp, -fyy > 0 then 
(a,b) is an inflection point 
elseif f,, -fyy <0 then 
(a,b) is a saddle point 
else { 
endif 
else 
test inconclusive 
endif 
else 
(a, b) is no critical point 
endif 


§ 
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Problem 15. There are two dead ends in Procedure 2. The first one (t+) is 
the case where fzz-fyy > (fry)? and either (f22 = 0, fyy = 0), (fee = 0, fyy < 0), 
(fea = 0,fyy > 0), fea < 0,fyy = 0), (fax > O,fyy = 0), (fea < 0, fy, > 9), 
or (fr2 > 0,fy, < 0). The second one (f) is where fz, - fy, = 0. Find out 
what happen in these cases, and thus complete the missing lines of logic in 
Procedure 2. 


§ 


Definition 24 gives the meaning of derivative and differential. Example 
13 looks at the derivative & differential of functions of one variable & two 
variables. 


Definition 24. By derivative gu we mean an infinitesimally small change in 
y with respect to an infinitesimally small change in x. By differential dy and 
dx we mean an infinitesimally small change in the values of y and respectively 


z. 


8 
Example 13. For a function of one variable y = f(z), the total derivative 
is 
dy 
dz 


and the differential of y is 


d 
dy = (34) dx 


For a function of two variables z = f(x,y) partial derivatives are, the first- 
order partial derivatives 


Oz a Oz 
Ox Oy 
and the second-order partial derivatives 
Oe 072 0? 072 


One = 225 ay? = Zyy; ayoe = tay 0d 5 ay = Zyx 


The total differential of z is 


of of 
dz = (=) dx + (F) dy 


and for small changes which are not infinitesimal, dz becomes Az and the 
incremental change formula is 


of of 
Azer {JA —)|A 
: (é) sal @ J 
Definition 25 gives the general production function. In Example 14 we 
look at the Cobb-Douglas production function in more details. Theorem 13 


gives the law of diminishing returns to labour and the proof thereof, while 
similarly does Theorem 14 the law of diminishing returns to capital. 
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Definition 25. The general production function is q = f(1,k), where q is 
output of the production, J labour and k capital. The Cobb-Douglas produc- 
tion function in its general form is 


gq = al®k? (1) 


where a is a constant and0O<a<1,0<8<1,l>0Oandk>0. 


§ 


Example 14. With the Cobb-Douglas production function, the marginal 
product of labour is, 


0 = 
Pim = QM = a = aal? 1k? (2) 


and the marginal product of capital 


0 2” 
Pkm = qk = an =apr ke © (3) 


From this we see that pj, > 0 and ppm > 0. 


Theorem 13. From the Cobb-Douglas production function we have the 
law of diminishing returns to labour, which states that qy < 0. 


Proof. From Equation 1 in Definition 25, 


aq _ 8 (Oq\ _ Opim aq 
aC e—— ll _—_ =e CO _ 1 — 
war” al (3) a ep 
Since 0 < a < 1, therefore q < 0. q 


Theorem 14. Using the Cobb-Douglas production function, the law of 
diminishing returns to capital states that qzp < 0. 


Proof. From Equation 1 in Definition 25, 


_ 4a _ 9 (% = OPkm _ (g _ 184 
akk = OR2 ~ Ok \ Ok Ok ke 
which, with 0 < 6 < 1, tells us that azz < 0. q 


In Example 15 we see the changes in marginal product values. 


Example 15. Using the Cobb-Douglas production function, 
gkt = Uk = aaBl?—* ko 


Therefore, qi, > 0 and qx; > 0. In other words, pj increases as capital input 
k increases, and respectively pm increases as labour input / increases. 


Example 16 shows us the average functions of labour and capital, Ex- 
ample 17 the marginal functions of labour and capital, and Example 18 the 
comparison between marginal and average functions. 
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Example 16. For the Cobb-Douglas production function in Equation 1 the 
average product of labour is 


a= ; =o ke (4) 
and the average product of capital is 


Pha = ; = al?k?—" (5) 
Example 17. Again using the Cobb-Douglas production function of Equa- 
tion 1, the marginal product of labour is 


0 
Pim = si = aal?—1 KF (6) 


and the marginal product of capital is 


0 = 
Pkm = aE = afleke (7) 


Example 18. From the APL equation, Equation 4, and the MPL equation, 
Equation 6, and since 0 < a < 1, therefore pry < pig. Similarly from the APK 
equation, Equation 5, and the MPK equation, Equation 7, since 0 < 6 < 1, 
we have prem < Dka- 


In Example 19 one sees the conditions for using labour, Equation 8, and 
the conditions for using capital, Equation 9. 


Example 19. A producer likes to have a positive marginal function, which 
means that the productivity increases as the input increases. But the second 
derivative is negative, which means that this rate of increase slows down as 
time goes by. In practice, the conditions for using labour are, 


= Oq dpim = a?q 
Pim = 5 > 0; “GF = Be <0; and Pim < Pla (8) 


The conditions for using capital are similarly, 


a Oq dpm = d7q 
Pkm = ak > 0, Thi Ok? <0, and prem < Pra (9) 


Definition 26 and Theorem 15 deal respectively with production function 
graphs and slope of an isoquant. 
Definition 26. An isoquant is a graph in two dimensions, k = f(1), plotted 
to represent a production function q = f(l,k). The slope oe is called the 
marginal rate of technical substitution. The value of this slope at (Io, ko) is 
denoted by $f ie 
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§ 


Theorem 15. The slope of an isoquant is the ratio of the marginal products. 


Proof. The total differential of g = f(I, k) is 


_ (94 Oq 
dq = (Ft) dl + (3) dk 


Along any isoquant, dq = 0, therefore, 


0= (3) dl + (52) dk (10) 


This directly yield, after some manipulation, 


Definition 27. In the Cobb-Douglas production function equation, Equa- 
tion 1, let both inputs | and k change by the same proportion, and let » be 
the constant of this proportionality. Then gz = a(Al)°(Ak)°, which leads to 
go = A°+8q,. When a+ 8 = 1, the case is described as constant returns to 
scale, when a+ 8 <1 as decreasing returns to scale, and when a+ 6 >1 as 
increasing returns to scale. 


8 
Definition 28. A homogeneous Cobb-Douglas production function of order 
r is, 
f(Al, Ak) = ATE, k) 
where r = (a, 8). 
8 


The utility function in its general form is given in Definition 29. Defi- 
nition 30 gives the Cobb-Douglas utility function, Definition 31 the marginal 
utility, and Definition 32 the meaning of indifference curves. 


Definition 29. A utility function expresses utility as a function of goods 
consumed. In its general form this is, 


u=f(z,y) 


where z and y are the quantities of goods X and respectively Y consumed. 


§ 
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Definition 30. The Cobb-Douglas utility function is in its general form, 
u = ary? 


where a is a constant, and0<a<1,0<8<1,2>0andy>0. 


§ 


Definition 31. The marginal utility for a utility function with one variable, 
u = f(x), is su = Uy = Ugm. The marginal utility for a utility function with 
two variables, u = f(x,y), is 8 = ue = Uem and oe iy = tage 

8 
Definition 32. The indifference curve is a graph y = f(x) drawn to repre- 
sent a utility function u = f(x,y). Its slope gu is called the marginal rate of 
substitution. Setting the total differential equal to zero, 


Ou Ou 
0=du= ($) dx + (5) dy 


we find 
dy _ Uz 
da Uy 
and 
dy _ _Uam 
dx Uym 


§ 


Partial elasticities are describe in Definition 33, partial elasticities of 
demand, and Example’s 20 and 21, respectively the partial elasticity with 
respect to labour and the partial elasticity with respect to capital. 


Definition 33. Let a demand function be 


da = £(pa, Y; Po) (11) 


where qq is the quantity demanded of good a, p, the price of a, y consumer’s 
income, and pp, the price of another good b. Then, the price elasticity of 
demand is, 


a= Oda Pa 
ODa Wa 
The income elasticity of demand is, 
Soe Oda y 
“OY qa 


And the cross-price elasticity of demand is, 


— 240 Po 
° Ops Ta 
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§ 


Example 20. With the demand function as in Equation 11, the partial 
elasticity with respect to labour is, 


Oql 
Eql = al q 


And from Equation’s 6 and 4, this leads to, 


For the Cobb-Douglas production function, Equation 1, then ¢q = a. 


§ 


Example 21. Again, with the demand function as in Equation 11, the 
partial elasticity with respect to capital is, 


Oqk 
Eqk = dk q 
Then, from Equation’s 7 and 5, 
gps Pkm 
? Pka 


For the Cobb-Douglas production function, Equation 1, we have €,; = . 


§ 
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Exponential, logarithmic and nonlinear functions 
8'® November 2005 


Definition 34 gives a definition of a function. Example 22 gives some 
examples of this. In Definition 35 we talk about variables and parameters of 
a function. 


Definition 34. A function is an operator or a procedure which accepts a 
permissible input and transforms it into a unique output. The input is some 
nonempty set. If a function is defined to be y = f(x), then x is the input 
vector, y the output, and f(-) the function itself. 

8 


Example 22. If f(-) is the function of dressing, then its input is possibly a 
person and its output a dressed person. If f(-) is the function of making up, 
then the input is perhaps a girl and the output a made-up girl. 


Definition 35. Let y = f(a,x) be a function, where a is a set of all its 
parameters, and x a set of all its variables. Then y is its dependent variable 
and x;, for all i € x, are its independent variables. In other words, 2; vary, 
y follows, and a; could assume any value within the range of its permissible 
ones, but its value must be constant. 


§ 


Definition 36 and Example 23 address inverse function, and respectively 
operator and inverse operator. 


Definition 36. An inverse function is an expression of the independent 
variable in terms of the dependent variables. The inverse of the function 
f(-) is denoted by f~+(-). If £(-) is a function which admits one independent 
variable, namely xz, then one could express it as, 


y = f(z) (12) 

Its inverse function is then, 
fly)=2 (13) 
8 


Example 23. Both the function and its inverse may be thought of as being 
an operator operating on an input to produce an output. The function, 


y = f(z) 


is understood diagrammatically as, 
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while its inverse function, 
«=f *(y) 


is displayed as a diagram as, 


Theorem 16 is related to the domain and range of inverse functions. 
Example 24 gives some examples of inverse functions. 


Theorem 16. An inverse function must always be a one-to-one mapping. 


Proof. Let f(-) be a function. Then f(-) can be either one-to-one or many- 
to-one, and therefore f—!(-) could turn out to be either one-to-one or one-to- 
many. But since f~! is also a function, so for each of the values in its domain 
the corresponding value in its range must be unique. This means that in 
cases where f—! turns out to be one-to-many, some constraints must be put 
on its input in order to make the output one-to-one, which then makes all the 
outputs from f—1(-) one-to-one. q 


Example 24. Table 1 gives some of the functions and their corresponding 
inverse functions which are fundamental in mathematics. 


f(:) £7) 
addition subtraction 
multiplication division 
power root 
exponential logarithm 


Table 1 Some of the functions and their corresponding inverse functions. 
Then, letting a be a constant, Table 1 becomes Table 2 


fe) eet) 
rt+a «£-a 
n-a 2 
x ft 
a” log, « 


Table 2 The notational forms of functions and their inverses. 


in which division and logarithm are both undefined for a = 0. 
Definition 37 gives some of the basic building blocks of mathematics. 
Example 25 then shows how these are built one on top of another. 


Definition 37. The inverse of the addition, 


y=at+a 
is the subtraction, 
y-a=2 
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The inverse of the multiplication, 


y = ax 
is the division, 
Log 
a 
The inverse of the power, 
y=ax° 


is the root, 


The inverse of the exponential, 


is the logarithm, 
log, y =% 
8 


Example 25. Figure 10 is a diagram which shows how addition makes 
multiplication makes power function. 


—— => a? 
b 
t 
g-x=2? 
t 
Bes SE 
See z-a 
a 
t 
g+u2= 224 


Figure 10 The building blocks of mathematics, from addition to multi- 
plication to power function. 


In Figure 10 we start from considering a variable at the base. If in- 
stead of doing this we begin by considering addition of some constant a, then 
eventually a becomes a parameter in our more complicated functions. 


a---a 
Sey - a 
x 
t 
at+...+a 
en 2 a-n 
n 
t 
a+t+az= 2a 
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Figure 11 Starting from a constant to obtain in the end the exponential 
function. 


Our derivation in Figure 10 gives us x’ when 0 is an integer, and similarly 
that in Figure 11 gives a” when z is an integer, but both the power- and the 
exponential functions can be extended to cover cases where the powers are 
noninteger, that is to say, when they are real or complex numbers. In these 
cases, however, the output may no longer be real. 

Next, we look at the exponential function, which is defined in Definition 
38. Example 26 discusses this further, and the exponential to the power of 
zero is looked at in Theorem 17. 


Definition 38. An ezponential function is defined as y = a”, where a > 0 
anda #1. 
8 


Example 26. The domain of the exponential function y = a® is the set 
of all real numbers, while its range the set of all positive real numbers. The 
function is convex and increasing when a > 1, and convex and decreasing 
when 0<a< 1. At 2 =0, the value of the function is y = 1 for any a > 0. 


Theorem 17. For any a #0, 


lim a* =1 
z—0 
8 
Problem 16. Try prove Theorem 17. 
8 


Theorem 18 gives some of the rules of exponential function and Example 
26 looks at some exponential function. 


Theorem 18. Three basic rules of the exponential function are, 


a™a”™ = qmtn 
m 

a =qrn 
qa” 


Proof. Write, say, a™ as, 


—— 
m 


and similarly for a”. Then all three equations above become obvious. q 
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Example 27. Figure 12 gives a graph of the exponential function when 
a> 1. Figure 13 gives a graph of the exponential function y = a” when 
0<a<l. 


ae 


y= 2? 
f 


95 2 x 0 1 2 3 
Figure 12 Example of the graph of the exponential function when a > 1. 
Here the graph is that of y = 2”. 


y = (0.5)* 
8 T T T 


0 L Ll L Ll 
-3 -2 -1 0 1 2 3 


Figure 13 An example of graph of the exponential function y = a” when 
0<a<1. Herea=0.5. 
From Figure’s 12 and 13, one may see that the graph of y = a*, where 
0 <a< 1, is the same as the graph of y = b *, where 4 =b> 1. This is 
obvious since by putting a = ; into y = a” one arrives at y = b ”, and if 
0<a<1,thend>1. 


One of the place we find use of an exponential function is the growth- 
and decay curves mentioned in Definition 39. Example 28 gives several basic 
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growth functions. 


Definition 39. Let a > 1. Then the graph of y = a” is called a growth 
curve, while that of y = a~* is called a decay curve. 


§ 


Example 28. There are basically three laws of growth, namely unlimited, 
limited and logistic growth, all of which involve an exponential function. The 
model is for unlimited growth, 


y(t) = ae" 


for limited growth, 
y(t) =m(1-e-”) 


and for logistic growth, 
m 


t) = ——_.] 
y( ) 1+ae—T™t 


where a, m and r are constants. 


Interests, and future- and present values are discussed in Example’s 29 
and 30. 


Example 29. The value of a principal p compounded annually at an interest 
rate 2 for t years is, 
s=p(1+i)® 


where 7 is expressed in decimal points. For compounding m times a year, 


then, 
j mt 
=p(1l4+— 
ica 


If the compounding is continuous, at 100 per cent interest for one year, then, 


mM—- Co 


1 m 
s=p lim (1+=) = pe 
m 
where e is the Euler’s constant, e = 2.71828... 


Example 30. For multiple compounding, 
-\ mt 
-\t 7 
p(1 + ie) =(1+ =) 
m 
the effective annual rate of interest is, 


j m 
i= (1+) -1 


The effective annual rate of interest for continuous compounding is, 
ie =e" -1 
Definition 40 and Example 31 are about discounting. 
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Definition 40. Discounting is the process of finding the present value p of 
a future sum of money s. 
8 


Example 31. Discounting when under annual compounding is, 
s=p(1+i)' 


when under multiple compounding, 


j —mt 
p=s(1+5) 
m 


and when under continuous compounding, 
p= se 


When discounting, the interest rate i is called the rate of discount. 


Example 32. A discrete growth s = p(1+i/m)"™ can be converted to a 
continuous growth s = pe” thus, 
mt 
) = pert 


p(1+ 
mt 
In (1+ ) =Ine”™ 
a 
r=mln (1+ =) 
m 


i mt F 
s=p (1 7 =) = pent x)t 


Example 33. Reversing the sign of x, that is replacing x by —2, has the 
effect of reflection of the original graph with respect to the y-axis. Reversing 
the sign of y, that is replacing y by —y, gives a reflection of the same with 
respect to the x-axis. The graphs of y = a+* remain always above the x-axis, 
in other words the function y = a** maps —oo < 4 < oo toy > 0. The two 
functions y = a® and y = a ® are the reflection of each other with respect to 
the y-axis. It can be easily seen that the functions y = —a** are the reflection 


with respect to the z-axis respectively of y = a**. 


3|-- 


3|-- 


Therefore, 


Definition 41 introduces the logarithmic function, and Example 34 gives 
some elaboration regarding this. Some examples of natural logarithm are 
given in Example 35. Theorem 19 gives rules for logarithm. 
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Definition 41. The logarithmic function with base a is defined to be the 
inverse of the exponential function, and is written y = log, x, where a > 0 and 
a #1. The logarithmic function of base 10 is called the common logarithmic 
function, and one of base e, where e = limn_yoo (1 + 4)" is called the natural 
logarithmic function. By the notation y = log, « we mean that the logarithm 
base a of x is the power to which a must be raised to get z. 


§ 


Example 34. The domain of the logarithmic function y = log, 2 is the 
set of all positive real numbers, its range the set of all real numbers. The 
function is concave and increasing for a > 1, and is convex and decreasing for 
0<a<1. Note also that log, x is the power which a must be raised to get 
x. 


Example 35. Note that e™* = a = Ine® where a > 0, e™® = & = Ine® 
where x > 0, and e!™) = f(x) = Ine**) where f(a) > 0. 


Theorem 19. Four basic rules for logarithm function are listed in the 
following. 
log, m + log, n = log, mn 


log, m — log, n = log, ~ 


log, m* = zlog, m 


log, n 
] See 
nen log, 6 
§ 
Problem 17. Prove Theorem 19, the theorem for rules of logarithm. 
§ 


Definition 42 sets out the meaning of the elasticity of substitution. Ex- 
ample 36 discuss the values of the elasticity of substitution. Definition 43 is 
on the constant elasticity of substitution production function. 


Definition 42. The elasticity of substitution o is defined as, 


az) = (4) 
oe ae) 
Coes: 

“7 a 


where f is called the least-cost input ratio, and eg the input-price ratio. 


§ 


Example 36. The value o = 0 means there is no substitutability, that is the 
two inputs are complements of each other and both must be used together in a 


30 14 January, 2007 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


fixed proportion. The value o = oo means that the two goods may substitute 
each other perfectly. Ultimately, 0 < 0 < co. 


Definition 43. <A constant elasticity of substitution production function 
is a production function where, unlike the Cobb-Douglas function, has an 
elasticity of substitution whose value is constant but not necessarily 1. In its 
typical form, it is, 


mex [ey 


q=a(ak? +(1-a)l-*)~ 

where a is called the efficiency parameter, a the distribution parameter, 6 the 

substitution parameter. Furthermore, § determines 0, anda >0,0<a< 1, 
and § > -1. 

8 


Example 37 discusses logarithmic transformation of nonlinear functions. 


Example 37. Some nonlinear functions can be converted to linear functions 
using logarithmic transformation, for example the Cobb-Douglas production 
function, 
q = ak*l? 
which becomes 
Ing=lIna+alnk+lnl 


Other nonlinear functions can not be converted, for example the constant 
elasticity of substitution production function, 


ai 
q=alak-° +(1-a)l-*] * 
which becomes just another nonlinear function, 
1 
Ing=Ina- as [ak-? + (1-a)l-F] 


others 


Example’s 38 and 39 give some examples of the use of nonlinear function 
in economics, namely the nonlinear total revenue and the nonlinear total cost. 


Example 38. Let the total revenue be r; = pq, and the demand function 
p =a-—bgq, where q is the quantity sold. Then r; expressed as a function of 
q is nonlinear, for r; = (a — bq)q = aq — bq’. 


Example 39. A more realistic equation for the total cost instead of c, = 
a+bgq is the nonlinear function c = aq? — bq? +cq+d in which the production 
cost increases with quantity in at a decreasing rate (c¢, < 0) up to the inflection 
point at q = a after which it increases at an increasing rate (c/ > 0). During 
the first stage the cost per unit decreases once the initial investment has been 
spent. During the second stage the cost per unit increases since more capital 
needs to be invested in order to allow more production capacity. 


Definition 44 talks about polynomial, and Example 40 about quadratic 
equation. 
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Definition 44. A polynomial is an expression in the form )>; , age” *, 
Here n is called the order of the polynomial. If n = 2 the polynomial is 
known as a quadratic polynomial, if n = 3 a cubic polynomial, ifn =4a 
quartic, n = 5 a quintic and n = 6 a seztic. If we let p(x) be a polynomial, 
then a polynomial equation is the equation p(x) = 0. A polynomial function 
is a function of the form y = f(x) = p(z). 


8 
Example 40. The quadratic equation ax? + bx + c = 0 has the solutions, 
—bivb? — 4ac 
p= (14) 
a 


These solutions are called the roots of the quadratic equation. Equation 14 
is called the ‘minus-b formula’. The values of x obtained from the minus-b 
formula give the intersections of the graph of the quadratic function 

f(x) = p(x) =az? +br +e 


on the x-axis. The value b? —4ac determines how the graph of f(x) lies relative 
to the z-axis, that is, 


> 0, there are two x-intersections 
b? — 4ac< =0, the graph touches the x-axis at one point 
<0, the graph never touches the z-axis 


Furthermore, the graph reverses its direction with respect to the y-axis at 
the critical point where f'(x) = 0, that is when x = —p. Consequently the 


critical point is 
b b? ie 
2a’ 4a 


The graph of f(x) is symmetric with respect to the vertical line which passes 
through the turning point, that is to say, the line x = -- The y-intercept is 
at the point (0, c). 


Example 41 is about another form of nonlinear function, the hyperbolic 
function. 


Example 41. A hyperbolic relation is an expression of the form, 


(px +q)(ry+s)=t 


From this we obtain, 


Ge Cae 
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where p, q, r, s and ¢ are constants, hence so are a a b ai c=. In 
economics we sometimes find hyperbolic functions of the form, 
a 
= 15 
Uae a (15) 


For example, a demand function of a good may be given by, 
m 
qra=— 
Pp 


which leads to, 
m 


ear 


Pp 


where p and gq are respectively price and quantity demanded of a good, while 
m and a are constants. 


The graph of Equation 15 has the z-axis, that is the line y = 0, as its 
horizontal asymptote, and has the line « = —c/b as its vertical asymptote. If 
all the parameters are positive, then the curve in the first quadrant decreases 
with a decreasing rate. 
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Matrix 
15%” November 2005 
Definition 45. Let A= {aj}, B = {bij} and C = {cj} be three matrices. 
Then C = A+B is called the addition of the matrices A and B if ej; = aj; +bij 


for all i and j. 


§ 


Definition 46. Let A = (a;;) be an m x n matrix and B = (by) an n x p 
matrix. Then the product AB is an m x p matrix C = (cy) where, 


n 
ca = Sandu 
k=1 


where 1<i<mand1<Ii<p. 


§ 


Definition 47. The expression obtained by eliminating the n variables 
%1,...,2%, from n equations, 


411%, +...+ a1n%pn = 0 
(16) 
AniX1 +... + Anny = 0 


is called the determinant of this system of equations, Equation 16. The de- 
terminant of matrix A denoted by various different notations, for example 
det(A), |Al, >> (tai bacg ad ‘), D(a1b2c3 ons ‘), or |a1 bacg “ay I. 


§ 


Example 42. For a linear system of three variables, Equation 16 can be 
written as, 

a,x + agy + agz =0 

byx boy b3z =0 (17) 


t+ eyt+cz =0 


Eliminating x, y and z from Equation 17 gives us, 


aybec3 — a, b3cq + agb1c2 — agbic3 + a2b3c, — agbec, = 0 


Definition 48. <A minor Mj; of any matrix A is the determinant of a 
reduced matrix obtained by omitting the i** row and the j** column of A. 


§ 


Theorem 20. Determinant can be determined by, 


k 
| A] = > Aig Cig 
i=l 
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where Cj; is called the cofactor of ajj. The cofactor C;; can also be denoted 
as a, and, at 

Ci = (-1)' Mi; 
where Mj; is a minor of A. 


§ 


Problem 18. Prove Theorem 20,the theorem for finding the determinant 
of a matrix by Laplacian expansion. 


§ 


Definition 49. Any pairwisely ordered pair in a permutation p is called a 
permutation inversion in p ifi > j and p; < pj. 


§ 


Theorem 21. Determination of the determinant can also be determined 


by, 
| A | = Sop Il ix (i) 
w i=1 


where 7 is a permutation which ranges over all permutations of {1,...,n}, 
and I(7) is called the inversion number of 1. 


§ 


Problem 19. Prove the theorem for the determination of determinant by 
permutation, Theorem 21. 


8 
Theorem 22. Let a be a constant and A an n x n matrix. Then, 
|aA|=a"|A| 
[=Al=(—D" A] 
|AB|=|A||B| 
|] =|AA“|=|A||[A*|=1 
1 
Al= 
| | | Aq? | 
8 
Problem 20. Prove Theorem 22, the theorem on properties of determinant. 
8 


Definition 50. <A function in two or more variables is said to be multilinear 
if it is linear in each variable separately. 


§ 


Theorem 23. Determinants of matrix are multilinear in rows and columns. 


§ 
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Example 43. Consider an 3 x 3 matrix, 


a, a2 a3 
A=l/a4 a5 46 
a7 ag ag 


What Theorem 23 says about multilinearity of determinants amounts to say- 
ing that, 


at 0 0 0 a2 0 0 0 a3 
|A|=]a4 a5 a6|+]a@s a5 ae ]t+|as a5 a6 
a7 ag ag a7 ag ag a7 ag ag 
and 
a, a2 ag 0 a2 a3 0 a2 a3 
|A|=) 0 as ag|+}as a5 ag ]/+| 0 as ag 
0 ag ag 0 ag ag a7 ag ag 


Problem 21. Prove the theorem on the multilinearity of determinants, 
Theorem 23. 


§ 


Definition 51. A conformal mapping is a transformation that preserves 
local angle. The terms function, map and transformation are synonyms. 


§ 


Definition 52. A similarity transformation is a conformal mapping the 
transformation matrix of which is, 


A' = BAB! 


Here A and A’ are similar matrices. 


§ 


Theorem 24. Similarity transformation does not change the determinant. 


Proof. The proof for this is simply, 


1 
BAB Nabe tee aT a 


Example 44. 


|B-1AB —I| =| BAB - B-1IB| 


=|B-M(A-ANB| 
=|B||A-Ar|| BI 
=|A—AI| 
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Definition 53. Let A be a square, n x n matrix. Then the trace of A is, 


Tr(A) = ay 
i=1 


Definition 54. The transpose of a matrix A = {a;;} is AT = {a;;}. 


Definition 55. The complex conjugate of a matrix A = {a;;} is A = {a;;}, 
where G=p—qiifa=p+qi. 
8 


Definition 56. Let ¢(n) or $(x) be a positive function, and let f(n) or 
f(z) be any function. Then f = O(@) if |f| < A@ for some constant A and 
all values of n and x. Here O is called the big-O notation which denotes 
asymptoticity. The notation f = O(¢) is read, ‘f is of order ¢’. 


8 
Theorem 25. Some other properties of the determinant are, 
|A| =|A?| 
|A| =| A] 
|I+¢A| =1+Tr(A) + O(e’), for € small 
8 


Example 45. For a square matrix A, 

a. switching rows changes the sign of the determinant 

b. factoring out scalars from rows and columns leaves the value of the de- 
terminant unchanged 

c. adding rows and columns together leaves the determinant’s value un- 
changed 

d. multiplying a row by a constant c gives the same determinant multiplied 
by c 

e. if a row or a column is zero, then the determinant is zero 

f. if any two rows or columns are equal, then the determinant is zero 


Problem 22. Prove the properties of determinant given in Theorem 25. 


§ 


Theorem 26. Some properties of matrix trace are, 
Tr(A) = Tr (A™) 
Tr(A + B) = Tr(A) + Tr(B) 
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Tr(@A) = aTr(A) 


8 
Problem 23. Prove that, 
Ci eee oo 
8 
Theorem 27. 
(AB)? = BT AT 
Proof. Sa ‘i z 
(B A }y = (b ie (a ej 
= bpiajr 
= Ajk OK: = (AB) ji = (AB)}, 
q 


Definition 57. Let A be a square matrix. Then the inverse of A, if it 
exists, is 4—! such that, 
AAT=I 


Furthermore, A is said to be nonsingular or invertible if its inverse exists, 
otherwise it is said to be singular. 


§ 


Example 46. For a 2 x 2 matrix, 


the inverse of A is, 


= 1 d —b 
At= 
ad — be & | 


If A is a 3 x 3 matrix, then the inverse of A is, 
A“! = {det (mij)} 
|A| 


where mj; is a minor of A. 


If A is an n x n matrix, then A~! can be found by numerical methods, for 
example Gauss-Jordan elimination, Gaussian elimination, and LU decompo- 
sition. 
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Example 47. The Gaussian elimination procedure solves the matrix equa- 
tion Ax = b by first forming an augmented matrix equation [Ab] and then 
transform this into an upper triangular matrix [{a/,}b’], where aj, are all 
zero except when i < j. Then, 


1 k 
= ! rf, : 
vi a! b; ai; vj 
uw j=l! 


The Gauss-Jordan elimination procedure finds matrix inverse by first forming 
a matrix [A J], and then use the Gaussian elimination to transform this matrix 
into [I B]. The result matrix B is in fact A+. 

The LU decomposition forms from the matrix A a product LU of two matrices, 
one lower- while the other upper triangular. This gives us three types of 
equation to solve, namely when 7 < j,i = j andi > Jj, where i and j are the 
indices of row and respectively column of the matrix product. Then, 


Ax = (LU)x = L(Ux) =b 


Letting y = Ux we have Ly = b, and therefore, 


by 
__— 
aoa Fe 
y 4-1 
Y= le b; ad Do hiius 
j=l 
where 7 = 2,...,n. Then since Ux =y, 
Yn 
Ly = — 
Unn 
1 nm 
T= — | Yi S- UjgXj 
Nei a 
j=i+1 


wherei =n-—-1,...,1. 


Theorem 28. Let A and B be two square matrices of equal size. Then, 
(AB tSp 44% 
Proof. Let C= AB. Then B= A!C and A=CB™|, therefore, 
C= AB=(CB"')(A'C) =CB'A'C 
Hence CB! A~! = TJ, and thus B-1A~! = (AB)"!. q 
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Definition 58. The Einstein’s summation is the simplification of notation 
by omitting a summation sign, keeping in mind that repeated indices are 
implicitly summed over, for example 57, ainajj becomes ajgaij, and >; aiai 
becomes a;4a;. 


§ 


Definition 59. The multiplication of two matrices A = {a;;} and B = {bj;} 
is the matrix C = AB such that cj, = Aig DjK- 


8 
Theorem 29. The matrix multiplication is associative. 
Proof. 
[(ab)c];; = (ab) incr; = (Gubie) Cag = Gu (Die Cej) = Ga(be)ig = [a(be)];; 
q 


Example 48. From Theorem 29, which shows us the associativity of matrix 
multiplication, we could write the multiplication of three matrices as [abc];;, 
which is the same as writing ajbj,c,;. And this applies in a similar manner 
to the multiplication of four or more matrices. 


Theorem 30. If A and B are two square and diagonal matrices, then 
AB = BA. But in general matrix multiplication is not commutative. 


§ 


Problem 24. Prove Theorem 30, which is a theorem about non-commutati- 
vity of matrix multiplication. 


§ 


Definition 60. A block matriz is a matrix which is is made up of small 
matrices put together, for example, 


le > 


where A, B, C and D are matrices. 


§ 


Theorem 31. Block matrices may be multiplied together in the usual 
manner, for example, 


Cy Di Cy Do} | CyAg+DiCyp C1 Bo+ D,De2 


provided that all the products involved are possible. 


Problem 25. Prove Theorem 31. 
8 
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Definition 61. Let A = {aj} be an n x n matrix. Then A is called a 
diagonal matrix if aj4j = 0 when i A 7. Here 1 < i,j <n. In other words, 
a diagonal matrix has its components in the form aij; = cjdj;, where c; is a 
constant and 4;; is the Kronecker delta, 


§ 


Theorem 32. A square matrix A can be diagonalised by the transformation 
A = PDP, where P is made up of the eigenvectors of A and D is the 
diagonal matrix desired. 

8 
Problem 26. Prove Theorem 32, the theorem on matrix diagonalisation. 


§ 


Example 49. Matrix diagonalisation can greatly help reducing the number 
of parameters in a system of equations. For instance, the systems Ax = y 
when diagonalised becomes PDP-'x = y, that is Dx' = y’, where x! = P~!x 
and y! = P~'y. In this case, if A is an n x n matrix, we say that our new 
system obtained through the process of diagonalisation has canonicalised from 
n Xn to n parameters. 


Definition 62. A symmetric matrix is a square matrix A which satisfies 
AT=A. 
8 


Example 50. If A is a symmetric matrix, then A~! A? = J. 


Definition 63. Let A be a square matrix. Then A is said to be orthogonal 
if AA? = I. 


8 
Example 51. Definition 63 is the same as saying that A~' = A’. 
Theorem 33. A matrix A is symmetric if it can be expressed as A = 
QDQ’, where Q is an orthogonal matrix and D is a diagonal matrix. 

8 
Problem 27. Prove Theorem 33, the problem on symmetric matrix. 

8 


Example 52. Any square matrix A may be decomposed into two terms 
added together, that is A, + A, where A, is a symmetric matrix and A, an 
antisymmetric matrix, called respectively a symmetric part and an antisym- 
metric part of A. Furthermore, 


A= 5 (4447) 
and, 
A, = 5 (4-47) 
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Linear algebra 
22"4 November 2005 


Definition 64. Let A be a square, nonsingular matrix. Then the inverse 
matriz A~' of A is a unique matrix for which, 


AA Sat SA TA 


8 
Example 53. An inverse matrix may be found using the formula, 
Art ods AdjA 
| A| 


Example 54. Matrix equations of the form Ax = b can be solved with the 
help of the inverse matrix A~! as x = A~'b, where A is an n X n matrix, x 
a vector of size n whose components are variables, and b a vector of size n 
containing constants. 


Theorem 34. Let A be the coefficient matrix and A; a matrix formed 
from A by replacing the column of coefficients of x; with the column vector 
of constants. Cramer’s rule solves a system of linear equations through the 
use of determinants as follows. 


oe, = Me 
‘Al 
8 
Problem 28. Prove Cramer’s rule, Theorem 34. 
8 


Definition 65. Let a system of n functions not necessarily linear be 


yi = fi (@1,---, En) 


Yn = fn (@1,---, En) 


Then a Jacobian determinant comprises all the first-order partial derivatives 
of the system arranged in ordered sequence, that is 


Oy. |, Om 
0x1 Orn 
— | Oyty.+:0Yyn — a < 
| J | ~~ | 021,.-.,0an | = a 
Oyn  .,, OYn 
0x1 Orn 


§ 
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Theorem 35. Let a system of n equations be y; = fj (@1,.--,%n), 1 = 
1,...,n. If|J| = 0, then y; are functionally dependent. On the other hand 
if | J| #0, then y; are functionally independent. 


§ 


Problem 29. Prove Theorem 35. 
8 


Definition 66. A determinant | H | composed of all the second-order partial 
derivatives, with the direct partials on the principal diagonal and the cross 
partials off the same, is called a Hessian. In other words, let a multivariable 
function be z = f(x,y). Then the Hessian of z is 


Zz Zz 
| H| — | ve zy 
Zyn Zyy 
where Zzy = Zyz- Moreover, the first principal minor is |H1| = zzz and the 
second principal minor is 
_ |%ex Zay| _ 2 
| H2| = ~ = Zrv%yy — (Zey) 
ry yy 


§ 


Theorem 36. Let a multivariable function be z = f(x,y), and let the 
first-order conditions zz = zy = 0 are met. Then a sufficient condition for z 
to be at optimum is ZzgZyy > (ey)? together with zzz, Zyy < 0 in case of a 
maximum and Zz, Zyy > 0 in case of a minimum. 


§ 


Problem 30. Prove Theorem 36, the theorem on the optimality of a mul- 
tivariable function. 


§ 


Definition 67. From Definition 66, if | H, | > 0 and | H2| > 0 the Hessian 
| H| is said to be positive definite, and the second-order conditions for the 
minimum are met. If | Hi | < 0 and | H2| > 0 it is said to be negative definite, 
and the second-order conditions for the maximum are met. 


§ 


Algorithm 1 Procedure to test for the optimality of multivariable func- 
tions of two variables. 
z= f(x,y) 
find z, and zy 
if z, = 0 and z, = 0 then 
find 222, Zzy and Zyy 
find H; and H2 
if | H, | > 0 and | H2| >0 then 
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| H| is positive definite 
elseif | H,| <0 and | H2 > 0| then 
| H| is negative definite 
endif 
endif 


Definition 68. Let y = f(#1,...,2n) be function of n variables. Then the 
n*b_order Hessian for this function is 


Yr oct Yin 
|H| = : 
Yn1 **" Unn 


Then the first principal minor | Hy | is simply x11, and the i** principal minor 
is 
Yur ctr Yui 
| Hi| = 
Ya oct Yu 
8 


Theorem 37. Let y = f(21,.-..,2,) be function of n variables. Let the 
Hessian of y be represented by | H|. Then if all the principal minors of | Hi | 
are positive, then | H | is positive definite and the second-order conditions for 
a relative minimum are met. If the sign of the principal minors alternates 
between negagive and positive, then | H | is negative definite and the second- 
order conditions for a relative maximum are met. 


§ 


Example 55. For y = f(21, 22,23) the third-order Hessian is 


Yl =Y12=«—«Y/13 
|H|=|yo1 yoo Yaz 
Y¥31 Y32 = Y33 
where 
0? 0? 


=> = ———., andsoon 
Yi1 dx? >» Yi12 022021 ’ 
The first-, second- and third-order Hessian’s are respectively 


Yi1 Y12 «18 
and |H3|=|yo1 yoo Yy23 
Y¥31 Y32 = Y33 


Yi1l —Y12 


A => 5 H. 4 
|fil=yu, | A2| valoda 


If |H,| > 0, |H2| > 0 and | H3| > 0, then H is positive definite and the 
second-order condition for minimum is fulfilled. If | H,| < 0, | H2| > 0 and 
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| H3| < 0, then | H| is negative definite and the second-order condition for 
maximum is satisfied. 


Definition 69. A discriminant is a determinant of a quadratic form. Let 
the quadratic form be z = ax? + bay + cy”, which is in matrix form 


rate lly] [ir 


Then the discriminant is 


b 
a pn 
ipi=| 2 
2 
The first principal minor of the discriminant is |D ,| = a, and the second 
principal minor 
b 2 
as b 
Dz| = 2)=ac-— — 
| P2| | g c 4 


Theorem 38. Let a quadratic form be 
z= ax" + bry + cy” 


and let the discriminant of z be | D |. If | Di | > 0 and | D2| > 0, then | D| is 
positive definite and z > 0 for all x,y 4 0. If | Di| <0 and |De2| > 0, then 
| D | is negative definite and z < 0 for all z,y 4 0. 


§ 


Theorem 39. Let f(x,y) be a function subject to a constraint g(x,y) = k, 
where k is a constant. Then the optimisation of f can be done by first 
transforming f together with g into a new function 


and then solve the following equations, 


Fi, (,y, A) = 0 
Fy(a,y,) =0 
F,(z,y,A) =0 


to obtain the critical values x0, yo and Ao at which F and hence f are opti- 
mised. 


§ 


Problem 31. Prove Theorem 38, the definiteness of a function by the 
discriminant. Prove the theorem for constrained optimisation with Lagrange 
multipliers, Theorem 39 above. 


§ 
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Definition 70. In the constrained optimisation with Lagrange multipliers 
in Theorem 39 above, f is called an objective or origin function and F the 
Lagrangian function. 


§ 


Definition 71. Let f(21,...,2n) be a function of n variables subject to 
constraints g(@1,..-,%n). Let 


F(a1,-..,%n,A) = f(&,-.-,2n) +A(k — g(a1,---, Bn) 


Then the bordered Hessian | H | is defined as either 


Fy Fi: Fin 
Foy 92 
|H|=| : 
Fra Fan gn 
fn 8692 gn 0 
or 
0 91 gn 
! n Fu Fin 
|H| =]. a te 
This is simply the Hessian 
Fyy sss Fin 
Fra wo Finn 


bordered by the first derivatives of the constraint with zero on the principal 
diagonal. The order of a bordered principal minor being determined by the 
order of the principal minor being bordered, | H| = | H,| since in this case 
an n Xx n principal minor is being bordered. 


8 
Theorem 40. Let f(z1,...,%,) be a function of n variables subject to con- 
straints g(21,...,2n). Let | H| be the bordered Hessian defined in Definition 
71. Then if | H2|,...,|H»| < 0, then the bordered Hessian |H | is positive 


definite, and therefore is a sufficient condition for a minimum. If | H2| > 0 
| H3| < 0, | H4| > 0, and so alternatingly on, then | H'| is negative definite, 
which is a sufficient condition for a maximum. 


§ 


Problem 32. Prove Theorem 40. 
8 
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Example 56. Let f(x,y) be a function to be optimised subject to a con- 
straint g(z,y) = k, where k is a constant. Then the Lagrangian function 
becomes 


The first-order conditions for optimisation are F, = Fy = Fy = 0. The 
second-order conditions for optimisation can be expressed together as a bor- 


dered Hessian 


2 Frye Fry Ix 
|H|=|Fyce Fyy Gy 
Gz Jy O 
or 
= 0 Gz Iy 
9y  Fyo Fy 
Note 1. Theorem 39 gives the first-order conditions for optimising a func- 


tion subject to some constraints. Theorem 40 gives the second-order condi- 
tions for optimising a function subject to some constraints. 


§ 


Definition 72. A Marshallian demand function gives an expression of the 
amount of a good that a consumer will buy as a function of commodity prices 
and income available. It is derived by maximising the utility subjected to a 
budgetary constraint. 


§ 


Example 57. Let a utility be u = qiq2 which is subject to a constraint 
Pigi + peg. = 6, where b is the amount of income available, that is to say, our 
budget. Then the Lagrangian function is 


U=Hna+rA(b— pin — p2q2) 


The first partial derivatives are then 


U1 = q2—- ADi =0 (18) 
U2 = dq — Ap2 =0 (19) 
Ur = b—pigi — prge = 0 (20) 


where ui, u2 are respectively ug, and ug,. Simultaneously solving Equation’s 


18, 19 and 20 leads us to 
q2 


Pr 


a 
P2 


=X 


Hence q2 = qipi/p2 and qi = q2p2/pi and from Equation 20 we have, 


141 P2q2 
b=pig shepard = peg pipieee 
P2 Pl 
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which yield us, for q, and q2, the Marshallian demand functions which max- 
imise satisfaction of the consumer subject to income and prices. 


Next, we test the second-order conditions by firstly finding ui, = 0, u22 = 0, 
u12 = ua1 = 1, gi = pr and gz = po, which give us 


a 0 1 P1 
|H|=|1 0 pe 
Pi pe O 


which gives | H2| = 2pip2 > 0 Hence | H| is negative definite and thus u is 
maximised. 


Definition 73. The production process of producing one good usually 
requires the input of many other intermediate goods. Let x; be the total 
demand for product 2, and let b be the final demand for the product from the 
ultimate users. Then, 


j= 4%, +...4+ Aint, + 0; 


fori =1,...,n, where aj; is a technical coefficient which represents the value 
of input 7 required to produce one monetary unit’s worth of product j. If we 
consider the total demand for every one of the products, then 


x=Ax+b 


where 
Ly a1 -* Gin by 


Zn Qnt ""* ann bn 


It follows from this that x = (I — A)~'b. The matrix A is known as the 
matrix of technical coefficients. It is also known as the input-output table, the 
rows being the inputs and the columns the outputs.. The matrix J — A is 
known as the Leontief matriz. 


§ 


Example 58. In a complete input-output table, labour and capital would 
also be included as inputs. These give the value added by the firm. They 
are normally put as an extra row at the bottom of the matrix of technical 
coefficients A. The vertical summation of each column of the table is then 
equal to 1. 


Definition 74. Let A be a square matrix. Then a scalar \ such that the 
equation 
Av =v (21) 
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holds for some vector v # 0 is called an eigenvalue } of A, and the vector v is 
called an eigenvector of A corresponding to the eigenvalue 4. The eigenvalue 
is also known as the characteristic root, or the latent root, while the eigenvector 
is also known as the characteristic vector, or the latent vector. 


8 
Note 2. From Equation 21 it follows directly that 


(A—Al)v =0 (22) 


Then A—AI is called the characteristic matriz of A. Since v assumes a unique 
value and by assumption v # 0, it follows that A—AJ must be singular, which 
means that its rows must be a multiple of one another. Now A — XI is zero 
if and only if the characteristic determinant | A —I| of A is zero. In other 
words 


|A-AI] =0 (23) 


which is called the characteristic equation of A. With Equation 23 there will 
be an infinite solution for v in Equation 22. In particular, if v is a solution, 
that is if it is an eigenvector, so is kv for any k # 0. We force a unique 
solution by using the normalisation 


Soup =1 


Then the sign-definiteness of A can be determined from the characteristic 
roots »’s. Thus if all X’s are positive, then A is positive definite; and if 
negative, negative definite. Let at least one \ be zero, which is neither positive 
nor negative, if all the remaining ’s are nonnegative, then A is positive 
semidefinite; and if they are nonpositive, negative semidefinite. Lastly, if 
some of the 4’s are positive while others are negative, then A is indefinite. 


8 
Problem 33. Prove all the necessary details need to be proved in Note 2. 
8 
Note 3. We have seen in Note 2 how, having found A;, where i = 1,...,n, 


we find through normalisation the corresponding, unique v;. On the other 
hand if we have found first the v,’s, their corresponding ;’s may be found 
by first forming a transformation matriz 


T =[v1 .-- Val 
+ The word eigenvalue is a half-translation of the German word Eigen- 
wert. The latter means ‘appropriate value’ since Wert means ‘value’ and 


eigen means ‘proper or appropriate’. 
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and then the corresponding eigenvalues or the characteristic roots are ob- 
tained from 


fe Oh “ae. 0 
rrap=|° ” 
ee 
8 
Problem 34. Prove all the necessary details in Note 3. 
8 


Definition 75. The vector equation, Equation 21, has as its solutions the 
zero vector v = 0 together with all the corresponding eigenvalue-eigenvector 
pairs. The set of all the eigenvalues of A is called the spectrum of A. The spec- 
tral radius of A is then the largest of all the absolute values of the eigenvalues 
of A, that is to say, 

max |Aa| 


The set of all eigenvectors v;;, together with 0, forms a vector space called 
the eigenspace of A corresponding to j;. 


§ 


Bibliography 
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Examples for linear algebra 
14" January, 2007 


1. Solve using the Gaussian elimination, 


321 + 8x5 = 53 
621 + 22 = 50 


Solution. Write the equations in an augmented matrix, 


[4 B]=|5 8 Aa 


6 2 50 


Row 1 times $, 


1 & 358 
3 3 

le 2 EA 
1, 


Row 2 subtracted by 6 times of Row 


1 8 53 
3 3 
F —14 | 


: I 
Row 2 times —z, 


128 
01 4 
Row 1 subtracted by 8 times of Row 2, 
10 7 
01 4 
Therefore £; = 7 and %2 = 4. 
# 
2. Let 
1 -3 3 
A=[3 -5 3 
6 -6 4 


Find all the eigenvalues and a basis of each eigenspace. Can A be diago- 
nalised? Why? 


Solution. Form the characteristic matrix tJ — A and find its determinant to 
obtain the characteristic polynomial A(t) of A; 


t-1 3 -3 
A(t) =|tI-A|=| -3 t+5 —-3 | =(t+2)?(t-4) 
-6 6 t-4 
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The roots of A(t) are —2 and 4, hence the eigenvalues of A are —2 and 4. 


Next, find a basis of the eigenspace of the eigenvalue —2. Substitute 


¢ = —2 into the characteristic matrix tJ — A, thus obtain the homogeneous 
system 

-3 3 -3 x 0 

-3 3 -3 y |=} 0 

—6 6 -6 z 0 


In other words, 
—3z + 3y — 3z = 0 
—32 + 3y — 3z =0 
—6z + 6y —6z =0 


that is x -—y+z=0. The system has two independent solutions, for example 
x=1,y=1,z2=0and2=1,y=0, z=-1. Therefore u = (1,1,0) and 
v = (1,0, —-1) are independent eigenvectors which generate the eigenspace of 
the eigenvalue -2. In other words, u and v form a basis of the eigenspace of 
-2. This means that every eigenvector belonging to -2 is a linear combination 
of u and v. 


Similarly, find a basis of the eigenspace of the eigenvalue 4. Substitute 
t = 4 into the characteristic matrix tJ — A to obtain the homogeneous system 


a es ae a 0 
291 OS | ge Soho 
-6 6 0 z 0 


This is the same as writing 


32 + 3y — 3z = 0 
—3z + 9y — 3z =0 


—6x2 + 6y =0 
which can be reduced to 
zt+y-—z=0 
2y—z=0 


Since the number of independent equations is less than the number of variables 
by one, this system has only one free variable. Therefore any particular non- 
zero solution, for example x = 1, y = 1, z = 2 generates its solution space. 
Hence w = (1,1, 2) is an eigenvector which generates, and so form a basis of 
the eigenspace of the eigenvalue 4. 


Since A has three linearly independent eigenvectors, A is diagonalisable. 


# 
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Let P be the matrix the columns of which are the three independent 
eigenvectors, 


1 1 1 
P={1 0 1 
0 2 


-1 
Then 
—2 0 0 
P-'AP=[ 0 -2 0 
0 0 4 


The diagonal elements of P~! AP are the eigenvalues of A corresponding to 
the columns of P. 


Bibliography 
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Exercises for linear algebra 
14%” January, 2007 


3. Find values of the variables «, y and z for each of the following systems of 
equations, using Gaussian elimination. 


2a—2y-—z = 3 4a+3yt+2 = 5 

(i) ex-ytz = 2 (it) eae ee = 4 

( o+y+2z2 = 3 ( ety-z = 3 
2x+2y—7z = 10 3a+2y+z2 = 16 

(iit) —“%-y = 5 (iv) etyte = 5 
3a+2y+z2 = —-1 2a-y-z = 7 
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Linear programming 
29" November 2005 


Definition 76. <A problem of optimisation is one in which one tries to 
maximise or minimise a certain quantity called the objective, which depends 
on a finite number of variables. These may be either independent or related 
to one another through some constraints. 


§ 


Definition 77. A mathematical programme is an optimisation problem in 
which the objective and the constraints are given as functions or mathematical 
relationship. In other words, 


optimise: z= f(@1,..-,%n) 


subject to: gi(%1,.--,2n) 


IV ILIA 
= 
I 
co 
3 


Some constraints are explicitly stated as requirements, others are hidden con- 
ditions. These latter need to be pin-pointed through the study and under- 
standing of the model and its inputs. 


§ 


Definition 78. A linear programme is a mathematical programme all the 
functions involved of which are linear. This means that, 


f(®1,---,%n) = 101 +--+ + Cnt 
Gi(@1,---; mn) = At. +--+ + inEn 
where 7 = 1,...,m and c; and aj, 7 = 1,...,n, are constants. If there is 


an additional restriction on the input variables that they be all integers, then 
the optimisation problem is called an integer programme. A mathematical 
programme which is not a linear programme is said to be nonlinear. 


§ 


Definition 79. A quadratic programme is a mathematical programme in 
which all the constraints are linear and the objective function is in quadratic 
form, which is in general, 


n n n 
f(®1,---,%n) = So > cigaeay + So dia; 
i=1 


i=1 j=l 


where cj and d; are constants. 
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Definition 80. A linear programme is said to be in standard form if all 
the constraints are equalities and if one feasible solution is known. In other 
words, our problem is now 
optimise: z= cx 
subject to: Ax=b 
with: x>0 


§ 


Definition 81. One may change any linear programme into the standard 
form by adding a slack variable to the left-hand side of a constraint of the 
form }> aijxj < by to obtain 


n 
) Ayr; + Lp, _ b; 
j=l 


where pp, > n and k = 1,2,.... Similarly one may add a surplus variable 
to the right-hand side of a constraint of the form Taya; > b; to obtain 


D1 ijt; = be + hq 
nm 
> ijt — &q, = bi 
j=l 


where q > n and / =1,2,.... Next, all the slack and surplus variables are 
added to the objective function with zero coefficients. Then if we add an 
artificial variable to the left-hand side of each constraint where there is no 
slack variable, then the initial feasible solution is x9 = b, where x is the 
vector of slack and artificial variables. The artificial variables are added to 
the objective function with a large negative coefficient —M. 


8 

Definition 82. A set of n vectors of m dimensions {v1,...,Vn} is said to 

be linearly dependent if there exist some constants ay,,...,Q,, not all of which 
are zero, such that 

ayVv, +++: +anVn = 0 (24) 


It is said to be linearly independent if the condition in Equation 24 implies 
Oy = Ag =-''=a,=0. 


§ 


Theorem 41. Consider a set of n vectors of m dimensions. If n > m, then 
the set is linearly dependent. 


§ 


Problem 35. Prove Theorem 41. 
8 
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Definition 83. <A vector v is called a conver combination of vectors v1,..-, 
v, if there exist some nonnegative constants (1,...,2n, where 


Brt--+Bu=1 


such that 
v=fAivi +:+++ Bavn 


§ 


Definition 84. A set of m-dimensional vectors is said to be convex if for any 
two vectors belonging to the set the line segment between them also belongs 
to the set. 


§ 


Theorem 42. All points on the line segment joining any two vectors may 
be expressed as a convex combination of the two vectors. 


§ 


Problem 36. Prove Theorem 42. 
8 


Definition 85. A vector v is called an extreme point of a convex set if it 
can not be expressed as a convex combination of two other vectors in the set. 


In other words, an extreme point of a convex set K is a point x in K that 
cannot be written as « = dy + (1—9)z with 0 < 6 < 1, y and z in K, and 
y # z, that is to say, an extreme point is a point which is not an interior 
point of any line segment belonging to K. 

An equivalent definition of an extreme point is that x is an extreme point of 
a convex set K if K\ {x} is convex. 


§ 


Definition 86. <A metric space is a non-empty set X for which is defined a 
concept of distance. The distance d is called a metric on X, having such 
properties that, for any points x and y in X, we have d(z,y) > 0, and 
d(x, y) =0 implies x = y; d(x, y) = d(y, x); and d(a,y) < d(a,z) + d(z,y). 


Let X be a metric space with metric d, let A be a subset of X and let x be 
any point of X. Then the distance from x to A is defined as 


d(x, A) = inf {d(z,a) : a € A} 
whereas the diameter of A is defined as 
d(A) = sup {d(a1, a2) : a1 andaz € A} 
Then a set is said to be bounded if its diameter is finite. 
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Further, let zo be a point in X and r a positive real number. Then the open 
sphere S,(xo) with centre xo and radius r is the subset of X defined by 


S, (ao) = {x : d(x, x9) < r} 


A point x in X is called a limit point of A if each open sphere centred on x 
contains at least one point of A different from x. A subset F' of X is said to 
be closed if it contains all its limit points. 


§ 


Definition 87. A linear space, aka a vector space, is a non-empty set L on 
which is defined two binary processes, say addition and scalar multiplication. 
Addition is defined such that for any x, y and z in L, then x+y is again in L; 
ety=yta;c+(ytz) = (e@+y) +2; there exists a unique identity element 
0, aka zero element or the origin, such that x + 0 = x for every x; and there 
exists a unique inverse element —x for every x, such that x+(—2x) = 0. Scalar 
multiplication is defined with regard to scalars, some instances of which are 
real and complex numbers, such that for any scalar a and any z and y in L, 
ax is again in L; a(a+y) = axr+ay; (a+ B)t = ant Ba; (a8)x = a(Bx); 
and lx = x, where 1 is the identity for scalar multiplication. A normed linear 
space is a linear space on which is defined a norm, that is a function which 
maps each element z in the space to a real number ||z|| in such a manner 
that ||z|| > 0, and ||z|| = 0 if and only if 2 = 0; ||z +y|| < |lz|] + |ly||; and 
law] = Ja [le 


8 
Theorem 43. A normed linear space is a metric space. 

8 
Problem 37. Prove Theorem 43. 

8 


Theorem 44. Any vector in a closed and bounded convex set with a finite 
number of extreme points can be expressed as a convex combinations of the 
extreme points. 


§ 


Problem 38. Prove Theorem 44. 
8 


Definition 88. For two vectors, that is points, x and y in R”, we write 
x > y if and only if 7; > y; for all 1 < i <n. A system of m weak 
linear inequalilties in n variables can be written as Ax > b, where A is an 
m Xn matrix. A fundamental question concerning such system is whether it 
is consistent, that is to say, whether there exists some x such that Ax = b. 
A system may be inconsistent, or it may have a set of solutions which is 
unbounded. If we sketch our problem on a graph, we may see that it’s solution 
set is convez. 
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§ 


Theorem 45. The solution space of a set of simultaneous linear equations 
is a convex set the number of extreme points of which is finite. 


§ 


Problem 39. Prove Theorem 45. 
8 


Theorem 46. Let S be the set of all feasible solutions to the linear pro- 
gramme in standard form in Definition 80, in other words, S is the set of all 
vectors x that satisfy Ax = b and x > 0, where A is an m x n matrix. Then 
S is a convex set, and the number of its extreme points is finite. The ob- 
jective function attains its optimum, provided that one exists, at an extreme 
point of S. If m <n, then the extreme points of S have at least n — m zero 
components. 


8 
Problem 40. Prove Theorem 46. 


Algorithm 2 Procedure for finding basic feasible solutions. 


Input: Ax = b, Ais an m x n matrix, m <n, rank A =m 
[a, --- a,|/¢A 
(ja, +---+2,a, = b) ¢ (Ax =b) 
for i=m+1tondo 
2; <0 
endfor 
(@1,---,%n) © solve za, +--+ +2,a, =b 


Definition 89. The simplex method is a matrix procedure which solves 
linear programmes of the standard form as described in Definition 80 where 
b > 0. Starting from a basic feasible solution x9 we locate successively other 
basic feasible solutions giving better values for our objective. For minimi- 
sation programmes the method uses Table 3, for maximisation programmes 
the same table is also used but with the sign of entries in the bottom row 
reversed. 


8 
Table 3 Table used for minimisation programming in simplex method. 
Xo co b 
—cjb 
8 


Table 4 Description of the simplex method. 


God’s Ayudhya’s Defence 14 January, 2007 59 


Business Mathematics, notes and projections Kit Tyabandha, PhD 


while negative number exists in d do 


Locate the most negative number in the bottom row of the simplex 
table, excluding the last column. The column in which we find 
this number is called the work column. If more than one such 
column exist, choose one of them. 

Find the smallest of the ratios between the elements in the last column 
and the elements in the work column of the same row, if these 
latter are positive. The element in the work column that yields 
this smallest ratio is called the pivot element. If there are more 
than one of these, choose one. If none of the elements in the work 
column is positive, the programme has no solution. 

Using elementary row operations, convert the pivot element to 1 and 
reduce all other elements in the work column to 0. 

Replace the x-variable in the pivot row and first column by the z- 
variable in the first row and pivot column. This new first column 
then becomes the current set of basic variables. 


endwhile 

The optimal solution is one in which all the basic variables assume the cor- 
responding values in the last column, while the remaining variables are zero. 
The optimal value of the objective function is then the value of the last row 
and last column for a maximisation programme, and the negative of this value 
if the programme is one of minimisation. 


§ 


Algorithm 3 gives a procedure for the simplex method. Here M repre- 
sents a large positive integer, p a ratio, c a column, r a row, and xo contains 
all the basic variables. Also the last row in Table 3 is represented here by 
d=c! —cf'A ande=—cfb. 


Algorithm 3 Algorithm for the simplex procedure. 
j <0 
while there exists a negative number in d do 
jegjtl 
for i=1 ton do 
{c’s} < (column number of the most negative number in the bot- 
tom row) 
(work column) ¢ choose one of the {c}’s 
k + (work column) 
endfor 
Ppivot *— M 
c+ 0 
soln + 0 
for i=1 tom do 
if (aj), > 0 then 
soln + 1 
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(b;); 


ae 35) ik 
if p < Ppivoe then 
ret 
endif 
endif 
endfor 
if soln = 0 then 
no solutions exist 
endif 
convert} A, such that (a;),, = 1 and (aj),, =0, l<i<m,ifr 
(Xo), = Le 
endwhile 
for 1 =1 to m do 
(xo); + (bj); 
endfor 
for i=m+1ton do 
x30 
endfor 
a+ 6; 
if the programme is one of minimisation then 
ze —2* 
endif 


Definition 90. The two-phase method is a procedure modified from the 
simplex method to cope with cases when artificial variables exist in the initial 
solution x9, in order to minimise the round-off errors that occur in the cal- 
culation. The last row in Table 3 in this case is d= ec? —cf A = di + Md, 
and consequently we have Table 5 which is used here. Algorithm 3 is then 
firstly applied to the last row, and then again to those elements directly above 
the zeros in that row. When an artificial variable is removed from the first 
column of the table, it ceases to be basic and may be removed from the top 
row of the table together with the entire column under it. When the last row 
contains only zeros, it may be deleted from the table. The programme has 
no solution if non-zero artificial variables are present in the final basic set. 


§ 


Table 3 Table used for minimisation programming using the two-phase 
method. 


+ With the use of elementary row operations. 
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§ 


Definition 91. Given a linear programme in the variables 21,...,2n, there 
exists another linear programme associated with it, called its dual, which is in 
the variables w1,...,Wm-. The original programme is called the primal. The 
primal completely determines the form of its dual. The symmetric dual of a 
primal linear programme in the matrix form 
minimise: z=c’x 
subject to: Ax >b 
with: x>0 


is the linear programme 


maximise: z=b’w 


subject to: A’w<c 


with: w>0 
The dual variables w1,...,Wm are known as shadow costs. The unsymmetric 
dual of the primal 
T 


minimise: z=c x 
subject to: Ax=b 
with: x>0 
is 
maximise: z= b’w 
subject to: Alw<c 


The unsymmetric dual of the primal 


maximise: z=c’ x 
subject to: Ax=b 
with: x>0 
is 
minimise: z= b’w 
subject to: A’w>c 


§ 


Note 4. We may see from Definition 91 that the dual of a programme 
in standard form is not itself in standard form. These duals are said to be 
unsymmetric. 

8 


Theorem 47. If an optimal solution exists for either the primal or the dual 
programme, then the other programme also has an optimal solution. If the 
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duality is symmetric, then the two functions have the same optimal value. If 
the duality if unsymmetric, then the optimal value of each function can be 
derived from that of the other. 


§ 
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Examples for linear programming 
14% January, 2007 


A. Use simplex method. 


maximise: z= 2, +9%24+ 23 
subject to: 21 +222 + 323 < 9 
321 + 2x2 +223 < 15 
with: all variables non-negative 


Solution. 
Ly 1 
a 2 1 2k 0 9 z 
re _ —_ i = 4 
| Mart ae; 2920 ee as] | 
v4 0 
x5 0 
cf. 
T 


optimise: z=c x 
subject to: Ax=b 
with: x>0 


Note: (ce? — cj A) and —c¢b in case of a minimisation problem, whereas 
—(c? —cf'A) and cf b in case of a maximisation problem. 


Tableau 1; 
v4 9 
£5 15 
0 


The most negative number in the last row is —9. Therefore x2-column 
becomes the work column. And then, 


2 
ty 2 + positive > $=45 
ts 2 — positive 4 2=7.5 
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Since min(4.5,7.5) = 4.5, the value of z2 on the row corresponding to 24 
becomes our pivot element. Then carry out a series of elementary row op- 
erations, namely in order (I)2 ¢ (Z)i/2; (D2 «— (Di — 22; TIDe2 — 


Now the last row is all non-negative. Therefore x3 = 3, 2$ = 6, xj 
r3=23=O0and2*=% 
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Integer programming 
6°" December 2005 


Definition 92. Algorithms which change the boundary of the solution 
region in order to find the optimal solution of an integer programme are 
called cut algorithms. The branch-and-bound algorithm does this by splitting 
the solution region into two and then discard the one which does not contain 
the optimal solution. The Gomory algorithm, on the other hand, reduces the 
feasible region with the help of a new constraint without the region being 
splitted. 

8 


Definition 93. We call branching a process by which a programme whose 
solution contains a non-integral j < x; < k is made into two separate pro- 
grammes having the additional constraint 7; < j in one, and x; > k in the 
other, the objective together with all the constraints of the original problem 
of which remain the same. Here j and k are positive integers and j < k. 


§ 


Definition 94. In the branch-and-bound algorithm, if the objective is max- 
imisation, the value of the objective obtained when the first integral approx- 
imation occurs is said to be the lower bound for the problem, and if the 
objective is minimisation it is said to be the upper bound of the same. 


§ 


Algorithm 4 Branch-and-bound algorithm for integer programming. 


find first approximation 
while approximations not all integers do 

choose 2; from all non-integral variables such that 

min (|x; — |w;||, |e; — [z;]|) is maximised 

branch 

choose the branch whose value of the objective is maximum 
endwhile 
solution¢last approximation 


Example 59. (Problem 6.9; Bronson, 1982) 


maximise: z= 2, +2%24+ 23 
subject to: 2%, + 3%. + 323 < 11 


with: all variables non-negative and integral 
Solve by branch-and-bound algorithm. 


Solution. Draw a simplex table of Programme 1. 
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w= + = 3.6, cf = «3 = 2] = 0, 2” 2 Since 3 < 23 < 4, branch 


into two programmes, namely Programme 1 where x2 < 3, and Programme 2 
where x2 > 4. Consider first Programme 2. 


maximise: z= 2, +2%.4+ 23 
subject to: 2%, + 322+ 3a3 < 11 
%<3 
with: all variables non-negative and integral 


Use the simplex method in a tabulated form. 


maximise: z= 2, +2%24+ 23 
subject to: 2%; + 3%2 + 323 < 11 
t2>A4 
with: all variables non-negative and integral 
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Draw a table for the two-phase method. 


11 
4 
0 
as 
11 
1 3 
1 
0 3 
1 2 22 
a 0 1 2 0 0 a 
2 1 3. 
2 0 1 i 1 si — 34 


The coefficient parts of the row corresponding to the basic variable x¢ 


and the last row cancel each other. The optimal result is 23 = #, xj = 23 = 
x4 =o, 9, =U and 27 ==: 


Therefore the solution is a] = 1, #3 = 3, #3 = a3 =a} =—0, and z* = 


Algorithm 5 Gomory algorithm for integer programming. 


while solution not wholly all integers do 
choose one non-integral optimal approximation 
write a relation from the row where that variable is basic 
rewrite the relation to make all fractional coefficients 
some integer plus a proper fraction 
move all the fractions to LHS, and all the non-fractions to RHS 
write a new constraint as LHS > 0 
find the solution for the original problem together with 
the new constraint 
endwhile 


Example 60. (Problem 7.1; Bronson, 1982) 
maximise: z = 27, + x2 
subject to: 22; + 5%. <17 

3x21 + 2%2 < 10 


with: 21,2%2non-negative and integral 
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Use cut algorithm. 


Solve 


Solution. Find the first approximation of Programme 1 normally using the 
simplex method. 


L3 0 17 
X4 0 10 
0 


Since 42 < 41, we know that 3 is the pivot element, and therefore we 


replace the basic variable x4 with 2. 


31 
3 3 
10 
Ly 3 
20 
3 
x _ 10 x _ 381 SRP ae ee ps *x _ 20 : 
We have af = 3, 3 = 3, 73 = %] = 0 and z* = =. Since both 


xj and «#3 are non-integers, arbitrarily choose the former to generate a new 
constraint. Then our Programme 2 becomes, 


Fee tlle NOG ak 
v4 ta7 344 = 3 = T 
2 i lone 
=t2+5-%4->= x 
a : 
2 1 dio 
—x + -XYa- 
Ce ae 
le a 
geo gee ae 
2%2 +44 >1 


and our new programme becomes 


maximise: z = 2%, +%.+0%3+4+ 024 
11 2 31 


bject to: —2_ — =-%4 = — 
Subject to 3 2 374 3 
ke 1 10 
£ =—%2 —- -%#4 = — 
cae eee aa 
2%2 +244 >1 


with: all variables non-negative and integral 
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Then our first approximation of Programme 2 is xj 
ai = ef = 0, and z* = 43. Arbitrarily choose 23 
constraint. 


bal Hg 
Ho) gu gs = 2 

1 1 I as 
ieee ae 
1 1 Pine 
pt gre 


3 
17 
2 
L 
2 
L 
2 
0 

3 

17 

2 

L 

2 

13 

2 

1 17 


* * 
3, aD) 37 U3 oe) 
to generate the new 
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Then our Programme 3 becomes, 


maximise: 


subject to: 


with: 


2 = 24, + 224+ 0x3 + 0x4 + O25 


1 
fit teS 3 
Pe re ae 
3 Q°4 6 nD) 
42 | 
LQ g%4 g75 5 


G4—-%5>1 


all variables non-negative and integral 


We draw our table for this programme. 


r2 


r3 


x7 


This becomes 
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X1 
x2 


2&3 


3 
0 
8 
4 1 
6 


0 0 0 0 1 L 


The optimum point for Programme 3 is then, xf = 3, #3 = 8, #3 = 1, 
v5 = xf = x§ = 0 and z* = 6. Therefore the solution to the original problem 
Programme 1 is a] = 3, x3 = 0 at the objective value z* = 6. 


# 


Definition 95. A transportation problem involves m sources each of which 
supplies a;,i =1,...,m, units of a certain product, and n destinations each 
of which requires 6;, 1 = 1,...,n, units of the same. The problem may be 
stated as following. 


maximise: z=) Cig Lig 


subject to: ye = &: a=1,...,m 


pee, PH 1; N 


with: all x;; non-negative and integral 


The total supply and the total demand are assumed to be equal. Were this 
not so, a fictitious destination or a fictitious source is added. 


§ 


Definition 96. The north-west corner rule finds an initial basic solution for 
the transportation algorithm of the integer programming. It begins with the 
(1,1) cell in the m x n table, and allocates as many units as possible to 211 
violating neither the constraints of supply, that is the summation along each 
row, nor those of demand, that is the summation along each column. Then 
carry on moving for each step either right or downwards, until we reach the 
lower-right corner, 2mn- 


§ 


Definition 97. <A loop, which is a sequence of cells in the table used for 
finding the solution in the transportation problem, has the following proper- 
ties. 
a. each pair of consecutive cells is on either the same row or the same column 
b. no three, or in fact any odd-numbered, consecutive cells lie in the same 
row or column 
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c. the first and the last cells are on the same row or column 
d. the path along the loop is self-avoiding, that is no cells appear more than 
once in the sequence 


8 
Algorithm 6 Transportation algorithm. 


while optimal solution not attained do 
find an initial, basic feasible solution using, for instance, the North- 
west corner rule 
let either u; = 0 or vj =0 depending on whether the i**-row or the 
j*_ column 
has the maximum number of basic solutions 
find all u; and vj,i=1,...,mandj =1,...,n from uj;+v; = cj; for 
basic variables, 
and from cj; — uj; — vj; for non-basic variables 
improve the solution 
endwhile 


Note 5. In a transportation problem, optimal solution is achieved when 
Cij — Us — uy > O for all transportation costs per unit ci; of all non-basic 
variables. 


§ 
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Financial mathematics 
13%" December 2005 


Definition 98. The present value po or the principal is the amount initially 
borrowed or invested. The future value ~; is the principal after a period of 
time ft. 


§ 


Definition 99. Interest rates expressed per annum are called nominal rates, 
i. The annual percentage rate or effective annual rate ig is the equivalent 
annual rate of different interest rates variously compounded. 


§ 


Definition 100. A sequence is a list of numbers which follows a definite 
pattern. It is called an arithmetic sequence if each of its terms is obtained 
from the term immediately preceding it by an addition of a constant d, which 
is called the common difference. It is called a geometric sequence if each of 
its terms is obtained from the previous term by a multiplication of a constant 
r, the common ratio. 


§ 


Definition 101. A series is the sum of the terms of sequence. It is called 
a finite series is one whose number of terms is finite, otherwise it is called an 
infinite series. An arithmetic series or arithmetic progression is the sum of 
the terms of an arithmetic sequence. Likewise a geometric series or geometric 
progression is the sum of the terms of geometric sequence. 


8 
Theorem 48. The value of the n‘* term of an arithmetic series is 
Tn =a+(n—1)d 
The sum of its first n terms is 
Sn = = (2a +(n—-1)d) 
8 
Problem 41. Prove Theorem 48. 
8 


Theorem 49. The n‘® term of a geometric series is 
Ears 
The sum of the first n terms of it is 


Sn =atar+--»+ar"1 


alr) 
~  J=r 
_ a(r” —1) 
— p=) 
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When the number of terms approaches infinity, the summation in cases where 


r <1 becomes é 


Oa ay 


§ 


Definition 102. A simple interest is a fixed percentage of the principal 
paid to an investor each year. A compound interest is an interest paid on the 
principal plus any interest accumulated in previous years. 


8 
Theorem 50. The present value for simple interest is 
Pt = po(1 + it) 
where 7 is the interest rate and ¢ the time in years. 
8 
Problem 42. Prove Theorem 50 
8 
Theorem 51. The present value in the case of compound interest is 
pi = po(l + i)* 
8 


Note 6. The interest may be compounded more than once a year, for 
example biannually, quarterly, monthly, weekly, daily, or continuously. Each 
time period is called a conversion period or interest period. The interest rate 
applied at each conversion is i/m, where m is the number of conversion periods 
per year. The number of conversion periods over ¢ years is then n = mt. 


8 
Theorem 52. The present value at the end of n conversion periods is 
-\n -\ mt 
a a 
Pt = Po (1+ =) = Po (1+ =) 
m m 
where all the variables and parameters are as previously defined. 
8 
Problem 43. Prove Theorem 52. 
8 


Theorem 53. When the number of compoundings per year becomes very 
large, the present value becomes 


Pt = Po et 


Proof. Since ~ = po (1+ ae and lim oo (1 + a)" = e', we have the 
proof. q 
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Theorem 54. The annual percentage rate when compounding occurs m 


times per year is 
i m™m 
ne (ee ie 
m 


Problem 44. Prove Theorem 54. 


Bibliography 
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Examples for financial mathematics 
14% January, 2007 


5. Given the demand function for a good as p = 1800—3g. Find the coefficient 
of point elasticity of demand when p is 300, 900 and 1200. Describe in words 
each of the results. Find the percentage change if the price of good increases 
by 10 per cent. 


Solution. The coefficient of point elasticity of demand is eg = = Se and 
also €g = 5, where Agqg is percentage change in quantity demanded and Ap 
percentage change in price. 
Po 300 900 1200 
qo Gsor"ro) 500-300 200 
Ed =e -3 =f @ 
eal ae <1 1 al 
demand elastic unit elastic inelastic 
Aqa eqAp -2% -10% -20% 
# 
Bibliography 


Teresa Bradley and Paul Patton. Essential mathematics for economics and 
business. 2°¢ ed. 2002(1998) 


God’s Ayudhya’s Defence 14 January, 2007 717 


Business Mathematics, notes and projections Kit Tyabandha, PhD 


Integral calculus 
10" January 2006 


Definition 103. Let f(x) be a function, and let f'(x) be its derivative. The 
reverse process of differentiation is called antidifferentiation or integration. It 
gives us the original function, which is called the antiderivative or integral of 


f(#)- 


8 
Theorem 55. Let c, n and k be constants. Then 
a. 
[ke =krt+c 
b. 
/ dx=az+ec 

C. i 

n = n+1 =] 

fe dx mane +o, n# 

d. 

[otaramnete, z>0 
e. 

[etae=inle|+e 0O#Ax<0 
f. j 
[eae — a +¢ 
g. 
[ri@ac=s f faz 
h. 
[u@sto@) = [tears f oeae 

i. 

[-t@ar=- | fae 

8 


Definition 104. The approximation 57}, (f(x;) Ax’) of the area under a 
continuous curve A is called a Riemann sum. That area under the curve is 


n 
A= Jim, > f(aj)Ax; 
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§ 


Theorem 56. Let F(x) be the integral of f(x). We call the fundamental 
theorem of calculus the expression. 


b 
i) f(a)de = F(a)? =F) — F@) 


8 
Theorem 57. 
| “fois2= eve 
i J 
b. 
[ fe@ar=F@—F@ =0 
; [tear= [pears f° Hear, a<b<c 
d. 
ik f(a)da + [sae = if (f(x) + g(a))dax 

[ si@az=e [sean 

8 


Property d of Theorem 57 is used to find the area between two curves. 


Theorem 58. The process of integration by parts is 


} (f(a) -g'(@))de = f(«) - g(a) - - (g(2) - f'(@)) de 
Proof. From 


< (f(a) - g(x) = f(a) - g(x) + g(x) - f'(@) 


we have 


Ha) - g(x) = i (F(a) -g!(a)) de + i (g(2)- f'(w)) de 
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Examples for integral calculus 
14% January, 2007 


6. Find ic x? dx, Hae (4 — 2) dx and tee ade. 


Solution. 


# 

Consider f°, 4dx = F(1) — F(-1). Here F(1) = 1 and F(-1) = -1 

But coo = F(0t+) #4 F(0~) = —oo, which means that f(z) = 1/2? is not 

continuous at 0, which is a point in [—1,1]. Therefore the definite integral 
given does not exist. 


# 
7. Find f 37 dz 
Solution. (1) We know that if u is a differentiable function of x, then 
nt+l1 
pou 22 +C 
n+1 
provided that n 4 —1. Let u = 27 +1, then du = 2zdz. And then, 
2z 1 
eS fd 
/ YP +1 / we 
= frau 
3 2 
_ 375 3 
a5 (2741)? +C 
# 
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Solution. (2) Let u= %/z? +1, then 3u?du = 2zdz. 


2z 3u? 
| yre- [7 


8. Evaluate i 322/23 + Ida. 


Solution. (1) Let u = 2? +1. Then du = 32?dz, u(—1) = 0 and u(1) = 2. 


And then, 
1 2 
/ 3a° Va + de = | Judu 
eal 0 


Solution. (2) Again, substitute the u and du as above. Then find the 
indefinite integral, 


[seve rias = f vadu 


2 


3 
2 


= 5 (+1) +C 


And then consider the corresponding definite integral, 


1 1 
i 327/23 + 1dz = = (0° +1) 
1 


= -1 


_4y2 
~ 3 
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9. Find the average value of f(x) = V4 — x? dz on the interval [—2, 2]. 


Solution. The average value required is, 


1 2 7 
a — 72 =~ 
Ee haaaar 
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Exercises for Integral calculus 


14% January, 2007 


10. Find the definite integrals, 
J (202% + 324 — 62%) dx, f Hyde, f 8Vx—7dz, 
f 4c?" dt, and f 3x73 dx. 

11. Find the definite integral 


f (2e* _ 3e**) dt 


given an initial, or a boundary condition F(0) = 3. 


12. Find the values of the definite integrals, ie 527 dz, q. 4e3 dt, L 6x3, 
and Soy (Aen dt. 


13. Find a firm’s total revenue r; function, given the marginal revenue func- 
tion rm = —.2a@? — 1.32 +500. 


14. Let the present value be p = a~™ of the sum of money a to be received 


in the future when the interest is compounded continuously. Find the present 
value p,, of a stream of future income, that is to say, the money to be received 
each year for n years. 


15. Find the following integrals and definite integrals. 
f 803 dx JP 803 dx ive4? dx 


SP Ve + 2dx J xrat > 3dt 
f (10-202) de =f? (10-22?) de = ff Aw — 3) Fda 


fr a - on f 10et% dt ie 10e1-?¢ dt 

J fas d i —— dx Sa 7a dx 

ST sae J (Jet h+e) de i. setae 
f c(6x + 3)dz ii x(6x + 3)dx f (823 + 927) dx 

JP (803 +927) de f (x3 + 1) dx tie (x3 + 7) dx 
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Kit Tyabandha. Integral calculus practice. Practices for Business Mathematics. 10 


84 


fG+2)dx 
i. (308 tet e*) dz 
f(l-e-*) dt 
fr (EP +e-%# +1) at 
J (e* + 30? + $) dt 
ie e® dt 
f-5dz 
f zo dt 
f5e~7*dt 
2 zx 
1 aT ae 
fesdt 
ie (50% ++ c*) dz 
I ove+3) at 
fie za + 7d 
J (2ve— 3) 
f — Jada 
f (4 + 2c?) dx 
fP 30-1 dx 
f2(@-1)-7dx 


outline series, 1993 


Jan 2006, Bangkok, 2006 


f7@ +2)*dz 

f 2x(a + 3)? dz 
SP 1-e-**) dt 
f tdx 


Sp (ec? + 382 + 4) dt 


3 furdu 

iF —5dz 

f ve® dx 

a, be dt 

Sf (6va+ 3) dz 
Jp esdt 

f 378 dt 

St (2v@ + 4) dt 
[3Vedx 

SP (t2Ve-4) dx 
f (2 + e”) dx 
i (5 + 2e*) dx 
i wry dz 

is 2(2 —1)-7dz 


Reference 
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f (308 +¢ +e") dz 
i 2x(x + 3)? dx 
Sle +e *+ 4) dt 
Jp eae 

f edt 

nae udu 

| dee 

iN ze" dz 

J trae 

ir (6Va+ 2) da 

fe (5a ++ +e") dz 
fp 3e7! dt 

J avJx +7dz 

aie 3./rdx 

[Gus 

i (2 +e") da 
f3a-tde 


2 29 
J wry dt 
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Integral calculus 
7h February 2006 


Definition 105. 
a. f° f(x)dx =0 
b. fo f(a)da = — f" f(a) dx 


Theorem 59. 
a. fo kf(a)de =k f° f(x)de. And if k = —1, then 


[ —f(«)dz = -[ f(x)dx 


b. fv (f(a) £ g(x) dx = J? f(a)dx + f° g(a)dx 

c. If f(x) > g(x) on [a, b], then ibs f(a)dx > g(x)dx Let g(x) = 0. Then, 
f(x) > 0 on [a, | implies aM f(z)dz >0 

d. If max f and min f are the maximum and minimum values of f on [a, b], 
then 


b 
min f-(b—a) < [ f(x)da < max f - (b-a) 


e. If f if integrable on the intervals between a, b and c, then 
b c c 
/ f(x)dx +f f(x)dz = / f(x)dx 
a b a 


Theorem 60 is called 60. 


Theorem 60. If f is continuous on the closed innterval [a, b], then at some 
point c in the interval [a, }], 


b 
fle) = == fleas 


Theorem 61. The average- or mean value of an integrable function f on 
a, b] is 
1 b 


Theorem 62. If f has a constant value c on [a, }], then 


[ feae= [ctv =c-0 


a 
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Theorem 63. If f is continuous on [a,}], then the function F(x) = 
J° f (t)dt has a derivative at every point on [a,b] and 


dF d f* 
ere | 1e=f@ 


§ 


Theorem 64. If f is continuous at every point of [a,b] and F is an an- 
tiderivative of f on [a, bj, then 


b 
[ tear =F -F@) 
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Simultaneous equations 
24" January 2005 


Definition 106. We call system of equations equations which together 
describe a mathematical model. A system of equations of n equations and v 
variables is called an n x v system or a system with n x v dimensions. If n = v 
the system of equations is called an exactly constrained system, if n < v an 
under-constrained system, and if n > v an over-constrained system. 


§ 


Note 7. A unique solution to a system exists only if there are as many 
equations as variables, that is to say, if n > v. An under-constrained system 
may have an unlimited number of solutions or no solutions, but it may never 
have a unique solution. An exactly constrained or over-constrained system 
may have a unique solution, an infinite number of solutions, or no solution. 


§ 


Definition 107. The graph of a (2 x 2) linear system of equations comprise 
two straight lines. If the two lines intersect, then the point of intersection 
(%1, 41) satisfies both equations and therefore represents a unique solution of 
the system. If they do not intersect, then there are no solutions and the two 
corresponding equations are said to be inconsistent with each other. If the 
two equations have identical graph, then the system has an infinite number 
of solutions. Such equations are called dependent or equivalent equations. 


§ 


Note 8. Consider a (2 x 2) system of linear equations in the slope-intercept 
form, 
y= mx + by 


y = Moz + beg 


if M4, x ms then 
system has a unique solution 
else 
if b, # be then 
equations are inconsistent and the system has no solution 
else 
equations are equivalent and the system has infinitely many solu- 
tions 
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Differential equation 
31%* January 2005 


Definition 108. A differential equation (DE) is an equation which involves 
derivatives. An ordinary differential equation (ODE) is a differential equation 
in which there is exactly one independent variable. A partial differential 
equation (PDE) is one where there are at least two independent variables. 
The derivatives of an ODE are ordinary-, whereas those of a PDE are partial 
derivatives. 


§ 


Definition 109. Consider a differential equation. The order of it is the 
order of the highest derivative appearing in it. Its degree is the degree of 
the highest ordered derivative therein. A primitive is a relation between the 
variables that involves n essential arbitrary constants, which gives rise to a 
differential equation of order n. The n constants are called essential if they 
cannot be replaced by a smaller number of constants. 


§ 


Example 61. The differential equation y’” + 3(y”)* + 2y’ = sina is an 
ordinary differential equation of order 3 and degree one. The differential 
equation (y")* + (y')® + y = 2x is an ODE which has an order 2 and degree 
2: 


Problem 45. The problem of finding solutions of differential equations is 
essentially that of finding the primitive which gave rise to the equation. 

§ 
Example 62. The differential equation y'” = 0 has a primitive y = Ax? + 
Ba+C, y!" —6y" +11y’ —6y = 0 has y = Cye2* + Coe?” + Cae®, y? (y")? +y? = 
7 has @=CP 4 yt =r. 


§ 


Definition 110. Existence theorems give conditions by which one could 
determine whether a differential equation is solvable. A particular solution 
of a differential equation is one obtained from the primitive by assigning 
definite values to the parameters, that is to say, the arbitrary constants. A 
singular solution is a solution which cannot be obtained from the primitive by 
any manipulation of the arbitrary constants. The primitive of a differential 
equation is usually called the general solution of the equation. 


§ 


Definition 111. A differential equation is said to be variable separable if 
an integrating factor can be readily found. Such equation has the form 


fo(a) - go(y)da + fo(x)- a y)dy =0 
Through the use of the integrating factor 
dL 
f2(@) - g2(y) 
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the primitive of this is then 


fi(z) gy) i 
fo(x) a / g2(y) se 


§ 


Definition 112. A differential equation of the first order and first degree 
may be written in the form 


M(a,y)dz + N(x,y)dy =0 


If this such equation admits a solution f(z, y,C) = 0 where C is an arbitrary 
constant, then there exist infinitely many integrating factors xi(x, y) such that 
€(z,y)[M(a,y)dx + N(x, y) dy] = 0 is exact, and there exist transformations 
of the variables which render the latter separated. But since no general rules 
exist for doing this, the use in practice is still somewhat limited. 


8 
Definition 113. A function f(x,y) is said to be homogeneous of degree n 
if 

f (Aw, Ay) = A” f (x,y) 

8 

Note 9. The equation 
(aye + by + c1) dx + (agx + boy + co) dy = 0 
where a,b — agb; = 0, is reduced through the transformation 
dt—ad 
ax+bjy=t and dy= — 
1 
to the form 
P(a,t)dx + Q(x,t)dt =0 
8 


Note 10. The equation 
(aya + by + c1) dx + (agx + boy + co) dy = 0 
where a,b — agb; 4 0, is reduced through the transformation 
g=a'th and y=y't+k 


in which x = h and y = k are the solutions of the equations a,x +biy+c, = 0 
and agxz + boy + co = 0 into the homogeneous form 


(ayz' + by’) dx’ + (agz’ + bey’) dy’ = 0 
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8 
Note 11. The equation of the form 
y- f(xy)dz + x - g(xy)dy = 0 
through the transformation 
z udz—zdx 
ryY=2%, Y=, dy = 2 
x x 
into the form 
P(a,y)dx + Q(z, z)dz =0 
which is variable separable. 
8 
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Examples for differential equations 
14% January, 2007 


16. Solve 2?(y + 1)dx + y?(a —1)dy = 0. 


Solution. The integrating factor is ae by the multiplicative applica- 
tion of which the given equation becomes 


2 


y? 
d. —dy =0 
a—1 ea] “ 


And then from this, 


1 1 
14+ —— -14+—— = 
(2+ +z) det (u +o) ay 0 


1 1 
on Se = 
[ («+ toy )act f(y +5) a 0 


2 


2 
9 +2+In(2-1)+ 2 —-y+h@ +I = 


x? + Qe — 2y +In(@ —1)(y+1)=Ce 
(@ +1)? + (y—1)? + 2In(# — 1)(y+1) =C 
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Difference equation 
20% February 2005 


Definition 114. A difference equation gives the relationship between an 
independent variable and a dependent variable, which changes at fixed, equally 
spaced intervals in time. The order of a difference equation is the number of 
time intervals spanned within the equation. 


§ 


Example 63. In the difference equation y+, = 1-ly;, an independent 
variable is the time t while the dependent variable is the income y. The order 
is the span of ¢ intervals within the equation, which in this case is (¢(+1)—t = 1. 


Theorem 65. The general solution of a homogeneous first-order differ- 
ence equation is of the form y, = Aa’, where ¢ and y are respectively the 
independent and dependent variables. 


§ 


Note 12. Given a homogeneous first-order difference equation of the form 
Yyt+1 — by, = 0 and some conditions, we may find the parameters a and A of 
our general solution y; = Aa* by first finding a by substitution of this general 
solution for ¢ and t+ 1 in the difference equation, and then find A from the 
conditions given. 


§ 


Theorem 66. The stability of the solution to a difference equation in 
general form y; = Aa’ is, 


range of a time path of y, solution time path 
-wo<a<-l ab s++0 unstable alternates 
-l1<a<0 a’ +0 stable alternates 
0<a<l a’ +0 stable tends to zero 
l<a<o a’ > unstable tends to infinity 


§ 


Theorem 67. The solution of a non-homogeneous difference equation is the 
sum of a complementary function and a particular integral, that is yx = y-+Yyp.- 
The complementary function is the solution of the homogeneous part of the 
difference equation. The particular integral is a function which satisfies the 
full difference equation. 


§ 


Note 13. The general form of the particular integral is deduced from the 
right-hand side of the difference equation. In particular, if c and b are con- 
stants, 
right-hand side general form of particular integral 
€ Yp =k 
cbt Yp = kb! 
8 
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Examples for difference equations 
14% January, 2007 


17. Solve the difference equation y,41 — 1.1y; = 0 by iteration for year 2, 3, 
4 and 5, given the income in year 1 is 20,000 Bahts.. 


Solution. We have yi41 = 1.1y, therefore 


yi = 20000 

yo = 1.1(20000) = 22000 
y3 = 1.1(22000) = 24200 
ya = 1.1(24200) = 26620 
ys = 1.1(26620) = 29282 


# 


18. Write out the solution of the difference equation y+; — 1.ly; = 0 for t = 
1, 2, 3, 4 and 5 in terms of y;. Deduce the general expression for y; in terms 
of yi. Evaluate yso given yz = 20000. 


Solution. 
yo = lly 
yg = Lly2 = 1.1 (11m) = (1-1)? 
ya = 1.ly3 = 1.1(1.1)?y, = (1.1)3y1 
ys = Leys = 1.1(1.1)?y, = (1.1)4y1 
# 
In general, 
yt = (1-1) "ye 
# 
t = 40, y, = 20000; 
yao = (1.1)?920000 = 822895.56 
# 


19. Find a general solution of the difference equation y+; — 0.9y; = 0. If 
y2 = 100, find the particular solution. Evaluate y1, yo, y3, yeo and yso- 


Solution. The general form of the solution is y; = Aa’. Therefore y+. = 
Aa‘t!, Then the difference equation becomes 
Yiti — 0.94, = 0 
Aa‘*? — 0.9a° = 0 
Aa‘ (a — 0.9) = 0 
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Since A 4 0 and a’ # 0, therefore a — 0.9 = 0, that is a = 0.9. Thus the 
general solution is 


UL = A(0.9)' 
# 
From y2 = 100; 
Yy2 = A(0.9)? 
100 = A(0.9) 
A = 123.46 
The particular solution is then 
Yp = 123.46(0.9)! 
# 
Then we tabulate the required calculation, 
tye 
1 123.46(0.9) = 111.11 
2  123.46(0.9)? = 100.00 
3: 123.46(0.9)? = 90.00 
20 123.46(0.9)?° = 15.01 
50  123.46(0.9)°° = 0.64 
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Exercises for difference equation 
14” January, 2007 


20. For each of the following difference equations, state 
i The order of the equation. 
ii Whether the equation is homogeneous or not. 
(a) Pet1 —0-8p,=0 — (b) Yrs = 8 — Yet (c) Yr2 =80+% 
21. Solve each of the following difference equations for the indicated variable 
by the iteration method. 
a. Yeri — 0-8y; = 10, given y; = 1. Then find ys. 
b. pepo = 4pe41 — 8pz, given py = 20, pop = 18. Then find ps. 
c. pe = 0.6p4_1 + 80, given p; = 100. Find ps. 
22. An amount of money A is invested at r% compounded annually. 
a. Show that the value of the investment at the end of any year ¢ is given 
by the difference equation pi4i1 = (1+1r)p; 
b. Show that the general solution of the equation is p, = po(1 +r). 
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Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Function 


A function of one independent variable is a relation in the form y = f(x) such 
that there exists one and only one value of y in the range of f for each real 
number x in the domain of f. The variable y is called the dependent variable. 
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An implicit function is a function in which both dependent- and independent 
variables appear on the same side. 


An explicit function is one where the dependent variable is on the left hand 
side-, and the independent variable on the right hand side of the equation. 
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Polynomial 
Anx” +... + a,x + a9 


Polynomial equation 


Anat” +... +a,2 +a9 = 0 


Quadratic equation 


ax? + bea +c=0 
ae —b+ Vb? — dac 
ae 2a 
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Definition of limits 


If f(x) is a function which draws closer to a unique finite real number / for all 
values of x as the latter draws closer to a, but x # a, then J is called the limit 
of f(x) as x approaches a. In notation this is, 


lim f(x) =1 


na 
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Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition of limits 


For a function f(x), limp, f(z) = 1 if and only if for every « > 0, there exists 
5 > 0 such that |f(z) —1| < « whenever 0 < |x —a| <0. 
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Rules of limits 


limk=k 
za 
lim 2” = a” 
za 


jim kia) =k iim f(x) 


lim [f(z) + g(x)] = lim f(x) + lim g(x) 


za za 
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lim [£(¢) - g(¢)] = lim f(¢) - lim g(c) 


hn [Sty ma 


n 
woe 
lim [f(z)] = [im f(2)] , provided that n > 0 


, provided that lim g(x) £0 
za 


za 
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Definition of a derivative 
Let y = f(z). Then, the derivative of y with respect to z is, 
dy . fa+h) — f(x) 
dea h 
The various notations for the derivative include 


df(«) df 


dz ’ dz’ f'(«), y’, Dy and D(f(z)) 
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Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


General rules of differentiation 


Let u, v and w are functions of x, and c is a constant. 


d Gain ) du, du 
dz 4 dx dz 
Business mathematics, Graph and derivative, 25" October 2005 -9- From 21% 
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d du 
—(cu) = c— 
da ce) dx 
d bay ue au 
da d da 
i (uvw) = uu— + uw 2 +uw = 
da da da da 
a ae v (a) — 4 (3) 
dx \v vy 
d,» n-1 du 
—(u") = nu” — 
da Ci da 
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dy _ dy du 
dx dudz 

du 1 
dr da 
dx s 
dy _ an 
= id 
dx = 
de* = e” 
dx 
dlnz 1 
dx  & 


Business mathematics, Graph and derivative, 25°" October 2005 -11- From 21% 
October 2005 , as of 14*” January, 2007 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Laws of exponents 


Let p and q be real numbers, a and b positive numbers, and m and n positive 
integers. Then, 


a? - qt = q?t4 
aP fas 
Se a? q 
a’ 

(a?)4 = qP4 


a° = 1, provided that a 4 0 
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Kit Tyabandha, PhD Department of Mathematics, Mahidol University 
aP= as 
a? 
(ab)? = a? bP 
Ya = at 
Van = an 
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Laws of logarithms 
log, mn = log, m + log, n 
m 
log, — = log, m — log, n 
n 
log, m? = plog,m 
log, n 
log, n = —— 
Sa log, a 
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global maximum 


| 


SEE, | 
y : Ea \ | | 


\ \ 
\ \ 


eigliter 


global minimum 


(a) (b) 


Figure 1 (a) global mazimum-, and (b) global minimum points of a 
function. 
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we 
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Figure 2 Local minimum- and local maximum points of a function. 
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concave up 


curves down 


y’> O 


y’<0 


concave down 


curves up 


(a) (b) 


Figure 3 The two curvature types, namely (a) concave up (y" > 0) and, 
(b) concave down (y" <0). 
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y’>0 
y’<0 y>0 y'>0 ie 
y’<0 
(a) (b) 


Figure 4 The four possible curvatures in two dimensions, considering 
both y’ andy", (a) y" > 0 and, (b) y"” <0. 
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fe aa 
f y’<0 ae / 
[ j yee 
ag ees 
(a) (b) 


Figure 5 Inflection points, where y" = 0, (a) with y" increasing and, 
(b) with y"” decreasing. 
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y'=0 | y’=0 
ea y’>0 
— —e 
y’=0 y’ =0 
y’<0 y’<0 
(a) (b) 


Figure 6 Stationary inflection points, where both y’ = 0 andy" = 0, (a) 
with y"” increasing and, (b) with y" decreasing. 
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Figure 7 shows a plot of the case where a = 1, b= 2, c=3 andd=4. 


y =ax* + be? + cx +d, where a=1,b=2,¢=3and d=4 


Figure 7 The cubic function y = x? + 227 +3244. 
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Figure 8 shows the case wherea = 1,b=5,c=4andd=3 
ra ao vi 
see 3 ere S 


Figure 8 The cubic function y = x? + 527 +4a 4+ 3. 
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Figure 9 The hyperbolic function (x — 1)(y — 2) =3. 
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The general demand function is of the form 


da = f(D, Y, Ps; De, ta; @, ae .) 


where qq is the quantity demand of good 2, p the price of x, y the income of 
the consumer, ps the price of substitute goods, p, the price of complementary 
goods, t, the taste or fashion of the consumer, and a the advertisement level. 
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In its simplest case where all other factors are constant, the demand equation 
takes the form 
p=C1 — C244 


where p is the price-, while gg the quantity demanded of the good z, and c, 
and cz are positive constants. 
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The general supply function is of the form 


Is = f(p, c, po, te, N, 0, , na) 


where q, is the quantity supplied of good z, p the price of x, c the cost of 
production, po the price of other goods, t, the available technology, n the 
number of producers in the market, and o other factors, for example tax and 
subsidies. 
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The simplified relation for the supply is 


p=C1 + C2qs5 


where q, is the quantity of x supplied, and c, and c2 are positive constants. 
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Definition 15 Total cost and revenue 
The total cost c; comprises a fixed cost and a variable cost, that is 
Ce = Ce tly 
The total revenue r; is price times output, 
T, = pq 
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Definition 16 Marginal cost and revenue 
The marginal cost cm is the change in total cost caused by the production of 
an additional unit. 


The marginal revenue rm is the change in total revenue coming from the sale 
of an extra good. In other words, 


tm = d 
pte 
me "dq 


where q is the output. 
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Definition 17 Average cost and revenue 


The average cost c, is the total cost per unit output, 


Ct 
=> 


q 


The average revenue rq is the total revenue per unit output, 


Tt 
Tfo=— 


q 


Business mathematics, Graph and derivative, 25" October 2005 -30-— From 21% 
October 2005 , as of 14°” January, 2007 


God’s Ayudhya’s Defence 14 January, 2007 111 


Business Mathematics, notes and projections Kit Tyabandha, PhD 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


From Definition 15, r; = pq, and from Definition 17, r; = raq, therefore 
P—-Ta 
From Definition 17, cg = aS and from Definition 15, cg = cz + cy. Therefore 
Ca = Caf + Cav 


where the average fixed cost cas = = and the average variable cost cgy = me 
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Procedure 1 Total-, average-, and marginal graphs 


Given: c,(x) 
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for each f € {cz,ca,Cm} do 
find the critical values x, for f’ = 0 
find f” 
n+ [xe| 
for 1=1 to n do 
if f(2§) > 0 then 
f(x§) is convex and is the relative minimum of f 
elseif f’ (xf) < 0 then 
f(x) is concave and is the relative maximum of f 
endif 
endfor 
find inflection points xy from f” = 0 
endfor 
plot c:(x#), then ca(x) and cm(x) 
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The relationship between input and output is called a production function, 


q={f(l,k,r, te, s,e,.-..) 


where | is labour, k phical capital such as buildings and machines, r raw 
materials, t. technology, s land, and e enterprise. 
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Assuming a short period of time, then / becomes the only independent variable 
and the other remaining factors are parameters, that is fixed, and therefore 
q=f(). 


Then the marginal product of labour is 


_ Wg 
Pim = dl 
and the average product of labour is 
Pla = = 
la jj 
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de 


The marginal propensity to consume is pem = ag? 


The marginal propensity to save is pgm = s.. 


The average propensity to consume is Deg = c 
The average propensity to save is Psa = Pu 
Here y is the income, c the consumption, and s the saving. 


Business mathematics, Graph and derivative, 25" October 2005 -36— From 21% 
October 2005 , as of 14*” January, 2007 


114 14 January, 2007 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


The profit is 
T=Tt-—Ce 


At the break-even point 


that is 
Te = Ch 
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Calculus of multivariable functions 


Definition 21 functions of n independent variables 


A function y = f(@1,...,%n) is called a function of n independent variables 
if there exists one and only one value of y in the range of f for each tuple of 
real number (21,...,2%n) in the domain of f. Here y is called the dependent 
variable while 7;,i=1,...,n, the independent variables. 
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Definition 22 partial derivatives 
Let a multivariable function be 


y = £(21,.-.,2n) 
The partial derivative of y with respect to 7;, where 1 <i <n, is a measure 


of the instantaneous rate of change of y with respect to 2; while x; is held 
constant for all 7 4%, where 1 < j <n. This partial derivative is defined as 


Oy a f(..., 0; + Ag;,...) —f(a1,...,2n) 


Ox; Axio Aa; 


and can be written in either one of the following forms. 


Oy Of 
> a fe; (¥1,---,Xn), fe;, OF Ya; 
Ox,’ Ox; i (1 n)> fai Ya: 
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Theorem 9 product rule 
Let z = g(a, y) -h(z,y). Then, 


dz _ hy a 
dx OO Ox 


and 
Oz oh og 


dy ° dy ay 
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Theorem 11 generalised power function rule 
Let z = [g(a,y)]”. Then, 


OF» geet sO 
Ox me Ox 
and 5 5 
z _ ygnt 88 
dy "® Oy 
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Definition 23 second-order partial derivatives 


Let z = f(z, y). Then, the second-order direct partial derivatives are 
O (Oz ree 0 (Oz 
Ox \ Ox Oy \ dy 


072 ; 022 
foo, (fe)a; a2 and respectively fyy, (fy)y, dy? 


These are also written 
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The cross partial derivatives are 


8 (2) nq 2 (2 
Oy \ Ox or Ox \ Oy 


These are also written as 


o? 6? 
fry, (fe) y, Bee and respectively fyz, (fy)z, sao 
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Theorem 12 critical points 


For a multivariable function z = f(z,y) to be a relative maximum at (a,b) 
necessarily f,,f, = 0, and frz,fyy <0 and fro -fyy > (fey)? at that point. For 
the same at the same to be a relative minimum, necessarily f,,fy = 0, and 
for, fyy > O and fra-fyy > (fey)? there. Moreover, an inflection point is a point 
(a,b) at which fyz -fyy < (fey)?, and both fy, and fy, have the same sign. On 
the other hand, a saddle point is a point (a,b) at which firs -fyy < (fey)?, but 
foe and fy, are of different signs. 
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Procedure 2 Procedure for determining a critical point of a function 
with two independent variables 


Given z = f(x,y) and a point (a, 6), at this point, 
if f, = 0 and f, = 0 then 
(a, 6) is a critical point 
if fro -fyy > (fey)? then 
if fez <0 and fy, <0 then 
(a,b) is a relative maximum of z 
elseif fz, > 0 and fy, > 0 then 
(a,b) is a relative minimum of z 
else { 


endif 
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elseif fy. -fyy < (fey)? then 
if fo2 -fyy > 0 then 
(a, 6) is an inflection point 
elseif frz - fy, <0 then 
(a,b) is a saddle point 
else { 
endif 
else 
test inconclusive 
endif 
else 
(a, 6) is no critical point 
endif 
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Problem 15 details in the critical point procedure 


There are two dead ends in Procedure 2. The first one (+) is the case where 
fou: Sdyy. > (ig) end either (fg — 0,14 — 0), (ae = 0g <0). (ee = 
0,fyy > 0), fee < 0,fyy = 0), (fax > 0, fyy = 0), (fer < 0,fyy > 0), or 
(fre > 0,fyy < 0). The second one ({) is where fre - fy, = 0. Find out 
what happen in these cases, and thus complete the missing lines of logic in 
Procedure 2. 
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Definition 24 derivative and differential 


By derivative sit we mean an infinitesimally small change in y with respect 
to an infinitesimally small change in x. By differential dy and dz we mean 


an infinitesimally small change in the values of y and respectively 2. 
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Example 13 derivative and differential of functions of one variable and 
two variables 

For a function of one variable y = f(x), the total derivative is 

dy 
da 


_ (vy 
dy = (5) dx 


For a function of two variables z = f(x,y) partial derivatives are, the first- 
order partial derivatives 


and the differential of y is 
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and the second-order partial derivatives 
Oz Oz o Oz 


de? = 2e25 dy? = yy, dyda =eay ond a = 2yz 
The total differential of z is 
Of of 
dz={|—]d — jd 
«= (35) e+ (5) 


and for small changes which are not infinitesimal, dx becomes Ax and the 
incremental change formula is 


of of 
Aze (=) Ag + (5) Ay 
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Definition 25 general production function 
The general production function is q = f(l,k), where q is output of the pro- 
duction, | labour and k capital. The Cobb-Douglas production function in its 
general form is 
q=al*kF (1) 


where a is a constant and0O<a<10<6<1,!]>0andk>0. 
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Example 14 Cobb-Douglas production function 
With the Cobb-Douglas production function, the marginal product of labour 
is, 


0 
= st Saal (2) 


and the marginal product of capital 


0 es 
ime = = xn Saal ere (3) 


From this we see that pj, > 0 and ppm > 0. 
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Theorem 13 law of diminishing returns to labour 


From the Cobb-Douglas production function we have the law of diminishing 
returns to labour, which states that q; < 0. 
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Proof. From Equation 1 in Definition 25, 


d’>q_ 9 (3) OB: — 


Uae aan) Ok 


Since 0 < a <1, therefore q, < 0. 
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Example 15 changes in marginal product values 
Using the Cobb-Douglas production function, 


Get = Qk = aaPl?—1Ke-1 


Therefore, giz > 0 and qx > 0. In other words, pj, increases as capital input 
k increases, and respectively pz, increases as labour input / increases. 
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Example 16 average functions of labour and capital 


For the Cobb-Douglas production function in Equation 1 the average product 
of labour is q 
Pla = 7 = TKO (4) 


and the average product of capital is 


tie Salk (5) 


Business mathematics, Calculus of multivariable functions, 1°* November 2005 —19- 
From 23” October 2005 , as of 14°" January, 2007 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Example 17 marginal functions of labour and capital 


Again using the Cobb-Douglas production function of Equation 1, the margin- 
al product of labour is 


0 
Pim = si = aal?~*k? (6) 


and the marginal product of capital is 


0 = 
Phim = 5p = apleKe! (7) 
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Example 18 comparison between marginal and average functions 
From the APL equation, Equation 4, and the MPL equation, Equation 6, 
and since 0 < a < 1, therefore py < pia. Similarly from the APK equation, 
Equation 5, and the MPK equation, Equation 7, since 0 < 8 < 1, we have 
Pkm < Pka- 
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Example 19 conditions for using labour and capital 


A producer likes to have a positive marginal function, which means that the 
productivity increases as the input increases. But the second derivative is 
negative, which means that this rate of increase slows down as time goes by. 
In practice, the conditions for using labour are, 


_ Og dpim _ dq 
Pim = aI > 0, =a = Or < 0, and pim < Pla (8) 


The conditions for using capital are similarly, 


_ Og dpim _ 0q 
= Ok > VU, AL = Ok? < 0, and prem < Pka (9) 


Pkm 
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Definition 26 production function graphs 


An isoquant is a graph in two dimensions, k = f(1), plotted to represent a 
production function q = f(I,k). The slope 


dk 
dl 
is called the marginal rate of technical substitution. The value of this slope at 


(lo, ko) is denoted by 
dk 


A isk 


Business mathematics, Calculus of multivariable functions, 1°* November 2005 —23- 
From 234 October 2005 , as of 14°" January, 2007 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Theorem 15 slope of an isoquant 
The slope of an isoquant is the ratio of the marginal products. 
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Proof. The total differential of g = f(1,k) is 


_ (% oq 
dq = (3) di + (#) dk 


Along any isoquant, dg = 0, therefore, 


a= (20) a+ (2) a ao) 
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This directly yield, after some manipulation, 


dk ee 
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Definition 27 returns to scale 


In the Cobb-Douglas production function equation, Equation 1, let both in- 
puts / and k change by the same proportion, and let A be the constant of this 
proportionality. Then gz = a(Al)°(Ak)°, which leads to q = A°+8q,. When 
a+ =1, the case is described as constant returns to scale, when a+6 <1 
as decreasing returns to scale, and when a+ 6 > 1 as increasing returns to 
scale. 
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Definition 28 homogeneous Cobb-Douglas production function 


The homogeneous Cobb-Douglas production function of order r is, 
f(Al, Ak) = AE(L, k) 


where r = (a, f). 
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Definition 29 utility function 


A utility function expresses utility as a function of goods consumed. In its 
general form this is, 
u = f(a, y) 


where z and y are the quantities of goods X and respectively Y consumed. 
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Definition 30 Cobb-Douglas utility function 
The Cobb-Douglas utility function is in its general form, 


u=arry® 


where a is a constant, and0O<a<1,0<8<1,2>0andy>0. 
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Definition 31 marginal utility 
The marginal utility for a utility function with one variable, u = f(x), is 
S* = uy = Ugm- The marginal utility for a utility function with two variables, 
u = f(x,y), is gu = Uz = Urm and ou = Uy = Uym- 
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Definition 32 indifference curves 
The indifference curve is a graph y = f(x) drawn to represent a utility function 


u = f(x,y). Its slope gu is called the marginal rate of substitution. Setting 


the total differential equal to zero, 


Ou Ou 
0=du= ($2) dx + ($) dy 


we find 
OU 22 te 
dz Uy 
and 
dy _ —_ Uam 
dz Uym 
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Definition 33 partial elasticities of demand 


Let a demand function be 
da = f(pa, Y, Po) (11) 


where q, is the quantity demanded of good a, pg the price of a, y consumer’s 
income, and py the price of another good b. Then, the price elasticity of 
demand is, 
= 04a Pa 

OPa Wa 


Ed 
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The income elasticity of demand is, 


= Oda y 
OY qa 


Ey 
And the cross-price elasticity of demand is, 


— 24a Pb 
° Opp da 
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Example 20 partial elasticity with respect to labour 
With the demand function as in Equation 11, the partial elasticity with re- 
spect to labour is, 
_ Oql 


ca Ol q 
And from Equation’s 6 and 4, this leads to, 


Egi = 
q 
Pla 


For the Cobb-Douglas production function, Equation 1, then eq = a. 
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Example 21 partial elasticity with respect to capital 
Again, with the demand function as in Equation 11, the partial elasticity with 
respect to capital is, 


Oqk 
Eqk = Ok q 
Then, from Equation’s 7 and 5, 
Pkm 
E — 
uf Pka 


For the Cobb-Douglas production function, Equation 1, we have éq, = f. 
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Definition 34 function 


A function is an operator or a procedure which accepts a permissible input 
and transforms it into a unique output. The input is some nonempty set. If 
a function is defined to be 


y = f(x) 
then x is the input vector, y the output, and f(-) the function itself. 
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Example 22 function 


If f(-) is the function of dressing, then its input is possibly a person and its 
output a dressed person. 


If f(-) is the function of making up, then the input is perhaps a girl and the 
output a made-up girl. 
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Definition 35 variables and parameters of a function 
Let 
y = f(a, x) 
be a function, where a is a set of all its parameters, and x a set of all its 


variables. Then y is its dependent variable and 2;, for all i € x, are its 
independent variables. 


In other words, x; vary, y follows, and a; could assume any value within the 
range of its permissible ones, but its value must be constant. 
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Definition 36 inverse function 
An inverse function is an expression of the independent variable in terms of 
the dependent variables. The inverse of the function f(-) is denoted by f~1(-). 
If f(-) is a function which admits one independent variable, namely x, then 
one could express it as, 


y = f(a) (12) 


Its inverse function is then, 
fl_y)=2 (13) 
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Example 23 operator and inverse operator 


Both the function and its inverse may be thought of as being an operator 
operating on an input to produce an output. The function, 


y ={(a) 
is understood diagrammatically as, 


fe.) 


while its inverse function, 


is displayed as a diagram as, 
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Theorem 16 inverse function 


An inverse function must always be a one-to-one mapping. 
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Proof. Let f(-) be a function. Then f(-) can be either one-to-one or many- 
to-one, and therefore f~1(-) could turn out to be either one-to-one or one-to- 
many. But since f~* is also a function, so for each of the values in its domain 
the corresponding value in its range must be unique. This means that in 
cases where f—! turns out to be one-to-many, some constraints must be put 
on its input in order to make the output one-to-one, which then makes all the 
outputs from f~1(-) one-to-one. q 
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f(:) f*() 
addition subtraction 
multiplication division 
power root 
exponential logarithm 


Table 1 Some of the functions and their corresponding inverse functions. 
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xv-:-a ai 
x Vt 
a® log, x 


Table 2 The notational forms of functions and their inverses. 


in which division and logarithm are both undefined for a = 0. 
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Definition 37 some inverse functions 


The inverse of the addition, 
y=art+a 


is the subtraction, 
y-a=2 


The inverse of the multiplication, 
y = ax 


is the division, 
y 
— — pi 
a 
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The inverse of the power, 


is the root, 


The inverse of the exponential, 


is the logarithm, 


loggy = & 
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Example 25 building-blocks of mathematics 


Figure 10 addition makes multiplication makes power function 


—— => x 
b 
t 
L-x=2? 
t 
eo+...+2 
ee nd zr-a 
a 
ak 
gta = 24 
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Figure 11 Starting from a constant to obtain in the end the exponential 
function. 
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Definition 38 exponential function 


An exponential function is defined as 


zx 


y=a 


where a > O anda. 
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Example 26 exponential function 


The domain of the exponential function y = a” is the set of all real numbers, 
while its range the set of all positive real numbers. The function is convex 
and increasing when a > 1, and convex and decreasing when 0 < a < 1. At 
x = 0, the value of the function is y = 1 for any a > 0. 
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Theorem 17 exponential to the power of zero 


For any a £ 0, 


lim a® = 1 
z—-0 
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Theorem 18 rules of exponential function 


Three basic rules of the exponential function are, 


a™a” = qm 
q™ 
pa, m—-n 
an a 
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Proof. Write, say, a™ as, 
a ee a 


—— 
m 


and similarly for a”. Then all three equations above become obvious. q 
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Example 27 graphs of the exponential function 


y= 
8 


0 
S3 ~2 a 0 1 2 3 


Figure 12 Example of the graph of the exponential function when a > 1. 
Here the graph is that of y = 2°. 
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y = (0.5)* 
\ 
\ 
a \ 
a ‘ 
. er 


Ss 2 a 0 1 2 3 
x 


Figure 13 An example of graph of the exponential function y = a” when 
O0<a<1. Herea=0.5. 
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Definition 39 growth- and decay curves 


Let a > 1. Then the graph of 
y=a" 


is called a growth curve, while that of 


Tz 


ya 


is called a decay curve. 
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Example 28 growth functions 


There are basically three laws of growth, namely unlimited, limited and lo- 
gistic growth, all of which involve an exponential function. The model is for 
unlimited growth, 

y(t) = ae” 


for limited growth, 
y(t) =m(1—e7") 
and for logistic growth, 


y() = 1l+ae—7™m 


where a, m and r are constants. 
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Example 29 interest compounding 


The value of a principal p compounded annually at an interest rate i for t 
years is, 
s=p(1+i)' 


where 7 is expressed in decimal points. For compounding m times a year, 


then, 
j mt 
S=p (1 + =) 
m 


If the compounding is continuous, at 100 per cent interest for one year, then, 


1 m™m 

s=p lim (1+=) = pe 
m—-> oo m 

where e is the Euler’s constant, e = 2.71828... 
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Example 30 multiple compounding 
For multiple compounding, 


j mt 
p(l+ie)' =p (2 + =) 
m 
the effective annual rate of interest is, 
4 m™m 
= (1 ie =) = 
m 
The effective annual rate of interest for continuous compounding is, 


de =e" —-1 
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Definition 40 discounting 


Discounting is the process of finding the present value p of a future sum of 
money 8s. 
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Example 31 discounting 
Discounting when under annual compounding is, 
s=p(1+i)' 


when under multiple compounding, 


j —mt 
p=s(1+=) 
m 


and when under continuous compounding, 
p=se™ 


When discounting, the interest rate i is called the rate of discount. 
Business mathematics, Exponential, log and nonlinear functions, 8°” November 2005 
~26- From 28° October 2005 , as of 14*” January, 2007 


146 14 January, 2007 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Example 32 converting exponential- to natural exponential functions 


A discrete growth s = p(1+i/m)™ can be converted to a continuous growth 


s = pe™ thus, 
j mt 
et) ow 
m 


‘ mt 
In (1 + =) =Ine”™ 
m 


r=min (1+) 
m 


j mt ; 
s=p (1 % ) Speer en) 
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Example 33 reflections of graphs 


Reversing the sign of x, that is replacing x by —z, has the effect of reflection 
of the original graph with respect to the y-axis. Reversing the sign of y, that 
is replacing y by —y, gives a reflection of the same with respect to the z-axis. 
The graphs of y = a*” remain always above the z-axis, in other words the 
function y = a** maps —co < x < co toy > 0. The two functions y = a® 
and y = a ® are the reflection of each other with respect to the y-axis. It can 
be easily seen that the functions y = —a** are the reflection with respect to 
the z-axis respectively of y = a**. 
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Definition 41 logarithmic function 


The logarithmic function with base a is defined to be the inverse of the expo- 
nential function, and is written 


y = log, x 


where a > 0 anda #1. The logarithmic function of base 10 is called the 
common logarithmic function, and one of base e, where 


1 n 
e= lim (1 + *) 
n—->Co n 
is called the natural logarithmic function. By the notation y = log, « we mean 


that the logarithm base a of x is the power to which a must be raised to get 


az. 
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Example 34 logarithmic function 
The domain of the logarithmic function 


y = log, z 


is the set of all positive real numbers, its range the set of all real numbers. 


The function is concave and increasing for a > 1, and is convex and decreasing 
for 0 <a <1. Note also that log, x is the power which a must be raised to 
get x. 


Business mathematics, Exponential, log and nonlinear functions, 8° November 2005 
~30- From 28% October 2005 , as of 14" January, 2007 


148 14 January, 2007 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Example 35 examples of natural logarithm 
Note that 


el — g=Ine® 


where a > 0, 


e™® — x = Ine 


where x > 0, and 
el) — f(x) = In ef) 


where f(x) > 0. 
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Theorem 19 rules of logarithm 


Four basic rules for logarithm function are listed in the following. 
log, m + log, n = log, mn 


log, m — log, n = log, = 
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Definition 42 elasticity of substitution 
The elasticity of substitution o is defined as, 


d(# a(#) 
, _ ae) 
Oe, a(Zzt) 
Pik tr 
PL Pe 
Pk 


where & is called the least-cost input ratio, and = the input-price ratio. 
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Example 36 values of the elasticity of substitution 
The value o = O means there is no substitutability, that is the two inputs 
are complements of each other and both must be used together in a fixed 
proportion. 


The value o = co means that the two goods may substitute each other per- 
fectly. Ultimately, 0 <a < o. 
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Definition 43 constant elasticity of substitution production function 


A constant elasticity of substitution production function is a production func- 
tion where, unlike the Cobb-Douglas function, has an elasticity of substitution 
whose value is constant but not necessarily 1. In its typical form, it is, 

2a 

q=a(ak-? +(1—a)l-*) 

where a is called the efficiency parameter, a the distribution parameter, 6 the 
substitution parameter. Furthermore, § determines 0, anda >0,0<a<1, 
and § > —1. 
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Example 37 logarithmic transformation of nonlinear functions 


Some nonlinear functions can be converted to linear functions using logarith- 
mic transformation, for example the Cobb-Douglas production function, 


q = ak@l8 


which becomes 
Ing=Ina+alnk+/lnl 
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Other nonlinear functions can not be converted, for example the constant 
elasticity of substitution production function, 


q=alak*+(1- al-8]* 


which becomes just another nonlinear function, 
: -p 8 
Ing=lIna— an [ak ° +(1—a)l*] 
others 
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Example 38 nonlinear total revenue 
Let the total revenue be 
Tt = pq 
and the demand function 
p=a—lbq 
where q is the quantity sold. Then r; expressed as a function of q is nonlinear, 
for 
re = (a — bg)g = ag — bq? 
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Example 39 nonlinear total cost 
A more realistic equation for the total cost instead of 

c =a+bq 
is the nonlinear function 
co = aq? — be? +ceqt+d 
in which the production cost increases with quantity in at a decreasing rate 
(c, < 0) up to the inflection point at 
b 


~ 6a 
after which it increases at an increasing rate (c// > 0). 
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Definition 44 polynomial 


A polynomial is an expression in the form 


n 
) aye” * 
i=0 


Here n is called the order of the polynomial. If n = 2 the polynomial is known 
as a quadratic polynomial, if n = 3 a cubic polynomial, if n = 4 a quartic, 
n = 5 a quintic and n = 6 a seztic. If we let p(x) be a polynomial, then 
a polynomial equation is the equation p(x) = 0. A polynomial function is a 
function of the form 


y = f(x) = p(x) 
Business mathematics, Exponential, log and nonlinear functions, 8°” November 2005 
-40- From 28°" October 2005 , as of 14°” January, 2007 


God’s Ayudhya’s Defence 14% January, 2007 158 


Business Mathematics, notes and projections Kit Tyabandha, PhD 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Example 40 quadratic equation 
The quadratic equation ax? + bx + c = 0 has the solutions, 
= —b+ Vb? — dac 


ala eee (14) 


These solutions are called the roots of the quadratic equation. Equation 14 
is called the ‘minus-b formula’. The values of x obtained from the minus-b 
formula give the intersections of the graph of the quadratic function 


f(z) = p(a) = ax? + bx +e 
on the z-axis. 
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The value b? — 4ac determines how the graph of f(z) lies relative to the z-axis, 
that is, 


> 0, there are two z-intersections 
b? —4ac< =0, the graph touches the z-axis at one point 
<0, the graph never touches the z-axis 


Furthermore, the graph reverses its direction with respect to the y-axis at the 
critical point where f(z) = 0, that is when 
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Consequently the critical point is 


b 6b ie 
2a’ 4a 
The graph of f(z) is symmetric with respect to the vertical line which passes 


through the turning point, that is to say, the line 


at 
2a 


The y-intercept is at the point (0, c). 
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Example 41 hyperbolic function 
A hyperbolic relation has the form, 
(pc +q)(ry+s)=t 
From this we obtain, 


ene). 


~o as ‘ 
where p, qg, r, s and ¢ are constants, hence so are a = a b= ai c= 8, 
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In economics we sometimes find hyperbolic functions of the form, 


_ 4a 
~ be +e 


y (15) 


For example, a demand function of a good may be given by, 
m 
qtra=— 
p 


which leads to, 
m 


 qta 


where p and q are respectively price and quantity demanded of a good, while 


m and a are constants. 
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The graph of Equation 15 has the z-axis, that is the line y = 0, as its hori- 
zontal asymptote, and has the line 


as its vertical asymptote. If all the parameters are positive, then the curve in 
the first quadrant decreases with a decreasing rate. 
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Definition 45 matriz addition 
Let A= {aij}, B = {bij} and C = {cj} be three matrices. Then 


C=A+B 
is called the addition of the matrices A and B if 
Cig = Aig + bi; 


for all i and j. 
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Definition 46 matrix multiplication 


Let A = (aj;) be an m X n matrix and B = (b,;) an n x p matrix. Then the 
product AB is an m x p matrix C = (cz) where, 


n 
cu = y Gindel 
k=1 


where 1 <i<mand1<Il<p. 
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Definition 47 determinant 


The expression obtained by eliminating the n variables 7,...,2, from n 
equations, 
44121 +... + Ginn = 0 


(16) 
Qnit1 +... + Anntn = 0 


is called the determinant of this system of equations, Equation 16. The de- 
terminant of matrix A denoted by various different notations, for example 
det (A), |Al, >> (tay becg me -), D(ayib2c3 ae ), or |a1 becg vee |. 
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Example 42 determinant of systems of three variables 


For a linear system of three variables, Equation 16 can be written as, 


aye + agy +a3z =0 
bia + bey + b3z = 0 (17) 
ct + coy +c3z =0 J 


Eliminating x, y and z from Equation 17 gives us, 


a1 bece3 — a, bgcg + agbice — agbic3 + agb3c1 — agbec, = 0 
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Definition 48 minor 


A minor Mj; of any matrix A is the determinant of a reduced matrix obtained 
by omitting the i* row and the j*® column of A. 
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Theorem 20 Laplacian expansion 


Determinant can be determined by, 


k 
[Al =o ayy 
i=1 


where Ci is called the cofactor of aj;. The cofactor C;; can also be denoted 
as a4, and, 7 
Cig = (-1)" Mi 


where M,; is a minor of A. 
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Definition 49 permutation inversion 


Any pairwisely ordered pair in a permutation p is called a permutation inver- 
sion in p ifi > Jj and p; < pj. 
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Theorem 21 determination of determinant by permutation 


Determination of the determinant can also be determined by, 


n 
|A| = yen Il Qi, ni) 
Tv t=1 


where 7 is a permutation which ranges over all permutations of {1,...,n}, 
and I(z) is called the inversion number of 7. 
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Theorem 22 properties of determinant 


Let a be a constant and A an n x n matrix. Then, 


|aA| =a"|A| 
|-A| =(-))"| 4] 
|AB|=|A||B| 
|Z] =|AA*|=|A|[A“| =1 
1 
|A| = |A~*| 
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Definition 50 multilinearity 


A function in two or more variables is said to be multilinear if it is linear in 
each variable separately. 
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Theorem 23 multilinearity of determinants 


Determinants of matrix are multilinear in rows and columns. 
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Example 43 multilinearity of determinants 
Consider an 3 x 3 matrix, 


a7 ag ag 
What Theorem 23 says about multilinearity of determinants is the same as 
saying that, 


ay, 0 0 0 a2 0 0 0 a3 
|A|=|]a@, as ag |/+]as as ag|/t+]a, as ag 
a7 ag ag a7 ag ag a7 ag ag 
and 
a, a2 a3 0 a2 a3 0 a2 a3 
|A|=| 0 as ag|/+]as as ag |/+] 0 as ag 
0 ag ag 0 ag ag a7 ag ag 
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Definition 51 conformal mapping 


A conformal mapping is a transformation that preserves local angle. The 
terms function, map and transformation are synonyms. 
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Definition 52 similarity transformation 


A similarity transformation is a conformal mapping the transformation matrix 
of which is, 
A'= BAB"! 


Here A and A’ are similar matrices. 
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Theorem 24 similarity transformation and determinant 


Similarity transformation does not change the determinant. 


Proof. The proof for this is simply, 


1 
Peete ae Peel alee =e 
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Example 44 similarity transformation 


| B-!AB—I| =|B-1AB — B-10IB| 


=|B-1(A—I)B| 
=|B-1||A—AI||B| 
=|A-Al| 
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Definition 53 matriz trace 


Let A be a square, n x n matrix. Then the trace of A is, 


Tr(A) = S ii 
i=l 
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Definition 54 matrix transpose 
The transpose of a matrix 
A= {aij} 
is 
AT = {aji} 
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Definition 55 complex conjugate 


The complex conjugate of a matrix 
A= {aij} 


is a 
A= {aij} 
whereG=p-—qiifa=p+q. 
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Definition 56 big-O notaton 
Let ¢(n) or d(x) be a positive function, and let f(n) or f(x) be any function. 
Then f = O(@) if |f| < Ad for some constant A and all values of n and 
xz. Here O is called the big-O notation which denotes asymptoticity. The 
notation f = O(¢) is read, ‘f is of order ¢’. 
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Theorem 25 properties of determinant 


Some other properties of the determinant are, 
|A| =| A? 
|A| =] 4] 
|I+eA|=1+ Tr(A) + O(e’), for € small 
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Example 45 notes on determinants 


For a square matrix A, 

a. switching rows changes the sign of the determinant 

b. factoring out scalars from rows and columns leaves the value of the 
determinant unchanged 

c. adding rows and columns together leaves the determinant’s value 
unchanged 

d. multiplying a row by a constant c gives the same determinant mul- 
tiplied by c 

e. if a row or a column is zero, then the determinant is zero 

f. if any two rows or columns are equal, then the determinant is zero 
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Theorem 26 matriz trace 


Some properties of matrix trace are, 
Tr(A) = Tr (A*) 


Tr(A+ B) = Tr(A) + Tr(B) 
Tr(a@A) = a Tr(A) 
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Problem 23 matrix transpose 


Prove that, 
(ay Sas)" 
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Theorem 27 property of matrix transpose 


(AB)? = BT AT 


Proof. ae . e 
(BrA hes = (0") x (@ Des 


= dpiajr 
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Definition 57 matriz inverse 


Let A be a square matrix. Then the inverse of A, if it exists, is A~+ such 
that, 
AA'=f 


Furthermore, A is said to be nonsingular or invertible if its inverse exists, 
otherwise it is said to be singular. 
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Example 46 matrix inverse 


For a 2 x 2 matrix, 


the inverse of A is, 
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If Ais a3 x 3 matrix, then the inverse of A is, 


eee 1 
A = [4] {det (maj) } 


where mj; is a minor of A. 

If A is an n xX n matrix, then A~+ can be found by numerical methods, for 
example Gauss-Jordan elimination, Gaussian elimination, and LU decompo- 
sition. 
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Example 47 finding matrix inverse 
The Gaussian elimination procedure solves the matrix equation Ax = b by 
first forming an augmented matrix equation [Ab] and then transform this 
into an upper triangular matrix [{aj,} b’], where aj, are all zero except when 
i<j. Then, 
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The Gauss-Jordan elimination procedure finds matrix inverse by first forming 
a matrix [A J], and then use the Gaussian elimination to transform this matrix 
into [I B]. The result matrix B is in fact A~'. 
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The LU decomposition forms from the matrix A a product LU of two matrices, 
one lower- while the other upper triangular. This gives us three types of 
equation to solve, namely when 7 < j,i = Jj andi > Jj, where i and j are the 
indices of row and respectively column of the matrix product. Then, 


Ax = (LU)x = L(Ux) =b 
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Letting y = Ux we have Ly = b, and therefore, 


= by 
ae tii 


i-l1 
y= [bi So ligy; 
Lig = 


where 2 = 2,...,n. 
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Then since Ux =y, 


wherez1 =n-—1,...,1. 
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— Yn 
Lyn = — 
Unn 
n 
1 
ta Ya 5 Usjz j 
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Theorem 28 matriz 


14*” January, 2007 
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Hence CB-1A~! = J, and thus B-1A! = (AB). 


Department of Mathematics, Mahidol University 


inverse 


Let A and B be two square matrices of equal size. Then, 


(AB) =B7A7 


Proof. Let C= AB. Then B= A~!C and A=CB |, therefore, 


CHAR = (CB WAC) =]CB AG. 


q 
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Definition 58 Einstein’s summation 


The Finstein’s summation is the simplification of notation by omitting a sum- 
mation sign, keeping in mind that repeated indices are implicitly summed 
over, for example 7, ajnaij becomes 


Qin dij 
and 57, a;a; becomes 


ajiaj 
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Definition 59 matrix multiplication 


The multiplication of two matrices A = {aj} and B = {bij} is the matrix 
C = AB such that 
Cik = ijdjr 
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Theorem 29 associativity of matrix multiplication 


The matrix multiplication is associative. 
Proof. 
[(ab)c];; = (ab)inces = (andre) Cag 


dit (biKCKj) = Au (bc)ag = [a(be)];; 
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Example 48 matrix multiplication 


From Theorem 29, which shows us the associativity of matrix multiplication, 
we could write the multiplication of three matrices as [abc];;, which is the same 
as writing ajbj,cnj- And this applies in a similar manner to the multiplication 
of four or more matrices. 
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Theorem 30 non-commutativity of matrix multiplication 


If A and B are two square and diagonal matrices, then 
AB=BA 


But in general matrix multiplication is not commutative. 
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Definition 60 block matrix 
A block matrix is a matrix which is is made up of small matrices put together, 


for example, 
A B 
C D 


where A, B, C and D are matrices. 
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Theorem 31 block matrix multiplication 
Block matrices may be multiplied together in the usual manner, for example, 


A; By} | A2 Bo] _ | A142+BiC2. Ai Bo + BiD2 
Cy Dy} [Co De Ci42+ DiC. C,Bo.+ Di D2 


provided that all the products involved are possible. 
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Definition 61 diagonal matrix 


Let A = {aj} be an n x n matrix. Then A is called a diagonal matrix if 
aij = 0 when i A j. Here 1 < i,j <n. In other words, a diagonal matrix 
has its components in the form aij = cjd4j, where c; is a constant and 4;; is 
the Kronecker delta, 
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Theorem 32 matriz diagonalisation 


A square matrix A can be diagonalised by the transformation 
Aa PoP! 


where P is made up of the eigenvectors of A and D is the diagonal matrix 
desired. 
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Example 49 matrix diagonalisation 


Matrix diagonalisation can greatly help reducing the number of parameters in 
a system of equations. For instance, the systems Ax = y when diagonalised 
becomes 

PDP-'x=y 


that is Dx’ = y’, where x! = P~!x and y’ = P~'y. In this case, if A is an 
n Xn matrix, we say that our new system obtained through the process of 
diagonalisation has canonicalised from n x n to n parameters. 


Business mathematics, Matrix, 15¢” November 2005 —44-From 5%” Nov 2005 , as of 
14*” January, 2007 


178 14 January, 2007 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Definition 62 symmetric matrix 
A symmetric matrix is a square matrix A which satisfies 
AMSA 
Example 50 symmetric matriz 
If A is a symmetric matrix, then 
AAP = 7 
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Definition 63 orthogonal matrix 
Let A be a square matrix. Then A is said to be orthogonal if 
AAT =I 
Example 51 orthogonal matriz 
Definition 63 is the same as saying that 
Av} = At 
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Theorem 33 symmetric matrix 


A matrix A is symmetric if it can be expressed as 
A= QDQ* 


where @ is an orthogonal matrix and D is a diagonal matrix. 
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Example 52 symmetric matriz 


Any square matrix A may be decomposed into two terms added together, 
that is A, + Ag where A, is a symmetric matrix and A, an antisymmetric 
matrix, called respectively a symmetric part and an antisymmetric part of A. 
Furthermore, 


As =5 (44.47) 


and, 
A, = 5 (A — A’) 
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Definition 64 inverse matriz 


Let A be a square, nonsingular matrix. Then the inverse matric A~+ of A is 
a unique matrix for which, 


AA =T=A CA 
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Example 53 finding the inverse of a matriz 


An inverse matrix may be found using the formula, 


1 
Alt=—AdjA 
| A| 
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Example 54 solving linear equations with the inverse 


Matrix equations of the form 
Ax =b 


can be solved with the help of the inverse matrix A~! as 
x=A'b 
where A is an n X n matrix, x a vector of size n whose components are 


variables, and b a vector of size n containing constants. 
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Theorem 34 Cramer’s rule 


Let A be the coefficient matrix and A; a matrix formed from A by replacing 
the column of coefficients of 7; with the column vector of constants. Cramer’s 
rule solves a system of linear equations through the use of determinants as 
follows. 


o, = 4 
<= 
| A| 
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Definition 65 Jacobian determinant 


Let a system of n functions not necessarily linear be 


Y1 = fi (@1,-.-,2n) 


Yn = fn (L1,---,2n) 
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Then a Jacobian determinant comprises all the first-order partial derivatives 
of the system arranged in ordered sequence, that is 


Oy... Ou 
0x4 Orn 
| J| = | OY1 yy OYn = 
021,...,0tn : 7 2 
Oyn 6 .,  DYn 
On41 Orn 
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Theorem 35 functional dependence from Jacobian determinant 
Let a system of n equations be 
Yi= fi (@1,---,2n) 


a@=1,...,n. 
If | J| =0, then y; are functionally dependent. 
On the other hand if |. J| 40, then y; are functionally independent. 
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Definition 66 Hessian 


A determinant | H | composed of all the second-order partial derivatives, with 
the direct partials on the principal diagonal and the cross partials off the 
same, is called a Hessian. In other words, let a multivariable function be 


z= f(x,y) 
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Then the Hessian of z is 
Zye Zyy 


where Rey = Xyu- Moreover, the first principal minor is 
y y ? 


and the second principal minor is 


z z 
| Ho), =|" O88) = zeatyy — (Zey)” 
Zay @yy 
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Theorem 36 optimality of a multivariable function 
Let a multivariable function be 
z= f(x,y) 
and let the first-order conditions 
Zq = 2y = 0 
are met. Then a sufficient condition for z to be at optimum is 
Zanzyy > (ay) 
together with 
Zen, Zyy <0 
in case of a maximum and 
Zen, yy > 0 


in case of a minimum. 
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Definition 67 definiteness 
From Definition 66, 


if | H, | > 0 and | H2| > 0 the Hessian | H | is said to be positive definite, and 
the second-order conditions for the minimum are met. 


If | Hi | < 0 and | H2| > 0 it is said to be negative definite, and the second- 
order conditions for the maximum are met. 
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Algorithm 1 Procedure to test for the optimality of multivariable functions of 
two variables. 


z= f(x,y) 
find z, and zy 
if z, =0 and z, =0 then 
find zzz, Z2y and Zyy 
find Hi and He 
if | Hi| > 0 and | H2| > 0 then 
| H| is positive definite 
elseif | H, | < 0 and | H2 > 0| then 
| H| is negative definite 
endif 
endif 
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Definition 68 general Hessian 


Let y = f(21,...,2n) be function of n variables. Then the n‘*-order Hessian 
for this function is 
Yur 7") Yin 


|H|= 
Ynl 77° Ynn 


Then the first principal minor | H, | is simply 211, and the i‘” principal minor 
is 


Yur occ Yui 
| Hi| = : 
Ya cc: Yi 
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Theorem 37 positive definiteness through Hessian 


Let y = f(x1,.-.,2%n) be function of n variables. Let the Hessian of y be 
represented by | H |. 

Then if all the principal minors of || are positive, then | H| is positive 
definite and the second-order conditions for a relative minimum are met. 

If the sign of the principal minors alternates between negagive and positive, 
then | H | is negative definite and the second-order conditions for a relative 
maximum are met. 
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Example 55 positive definiteness for three dimensions 
For 
y = f (21, £2, 23) 


the third-order Hessian is 


Yur 1218 
| H | = 1/421 Y22 Y23 

Y31 =Y32 = Y33 

where ; ; 
a oy = ms and so on 
we Ox?’ via 022021’ 
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The first-, second- and third-order Hessian’s are respectively 


[Fi |= yn, |He| =| oe 
Y21 —Y22 
and 
Yi Y1l2——s13 
| H3| = Y21 Y22 Y28 
Y31 =Y32 = Y33 
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If |Hi| > 0, |H2| > 0 and | H3| > 0, then # is positive definite and the 
second-order condition for minimum is fulfilled. 


If | Hi| < 0, | H2| > 0 and | H3| < 0, then | H| is negative definite and the 
second-order condition for maximum is satisfied. 
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Definition 69 discriminant 


A discriminant is a determinant of a quadratic form. Let the quadratic form 
be 
= 2 
z=axn~+bry+cy 


which is in matrix form 


a 2) [a 
cele ale ali] 
Then the discriminant is ‘ 
ely! 
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The first principal minor of the discriminant is 


|Dil|=a 
and the second principal minor 
b 2 
a> b 
D2|= 2} =ac-— — 
| D2 8 c | 4 
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Theorem 38 definiteness of a function by the discriminant 


Let a quadratic form be 
z= ax” + bry + cy? 


and let the discriminant of z be | D|. 


If |D,| > 0 and | D2| > 0, then | D| is positive definite and z > 0 for all 
x,y #0. 
If |D,| < 0 and | D2| > 0, then | D| is negative definite and z < 0 for all 
x,y #0. 
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Theorem 39 constrained optimisation with Lagrange multipliers 
Let 
f(x,y) 


be a function subject to a constraint 


g(a,y) =k 


where k is a constant. Then the optimisation of f can be done by first 
transforming f together with g into a new function 


F(a,y, ) = f(x,y) + A(k — g(a, y)) 
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and then solve the following equations, 


F,(a, y, A) =0 
Fy (x,y, X) = 0 
Fy(z,y, A) =0 


to obtain the critical values 29, yo and Ao at which F and hence f are opti- 
mised. 
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Definition 70 constrained optimisation 


In the constrained optimisation with Lagrange multipliers in Theorem 39 
above, f is called an 


objective or origin function 
and F the 


Lagrangian function 
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Definition 71 bordered Hessian 


Let 
f(v1,---,2n) 
be a function of n variables subject to constraints 
g(@1,---,Xn) 
Let 
F(a,.--,2n,r) = f(@,.--,2n) +A (K— g(a1,---, En) 
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Then the bordered Hessian | H | is defined as either 
Fu Fin s+ Fin 
Foy 92 
|H|=| : : : 
Frnt Fan Qn 
fn 92 gn 9 
or 
O got Gn 
= gn Fu Fin 
Le Sy3 ; 
Gn Fut ean Fun 
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This is simply the Hessian 


bordered by the first derivatives of the constraint with zero on the principal 
diagonal. The order of a bordered principal minor being determined by the 
order of the principal minor being bordered, 


|H| =|An| 


since in this case an n Xx n principal minor is being bordered. 
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Theorem 40 bordered Hessian 
Let f(x1,.-.,%n) be a function of n variables subject to constraints 


G(@1,---,%n) 


Let | H | be the bordered Hessian defined in Definition 71. 


Then if | H2|,...,|_Hn| <0, then the bordered Hessian | H | is positive defi- 
nite, and therefore is a sufficient condition for a minimum. 


If | H2| >0|H3| <0, | H4| > 0, and so alternatingly on, then | A | is negative 


definite, which is a sufficient condition for a maximum. 
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Example 56 constrained optimisation using bordered Hessian 


Let f(x,y) be a function to be optimised subject to a constraint g(x,y) = k, 
where & is a constant. Then the Lagrangian function becomes 


The first-order conditions for optimisation are 


F, =F, =F, =0 
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The second-order conditions for optimisation can be expressed together as a 
bordered Hessian 
= Fox Fryy 9x 
|H |=) Fyz Fyy Gy 
Ga Dy 0 
or 
0 Ga Dy 
Iy Frye Fyy 
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Note 1 first- and second-order conditions for optimisation of a function 
subject to some constraints 
Theorem 39 gives the first-order conditions for optimising a function sub- 
ject to some constraints. Theorem 40 gives the second-order conditions for 
optimising a function subject to some constraints. 
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Definition 72 Marshallian demand function 


A Marshallian demand function gives an expression of the amount of a good 
that a consumer will buy as a function of commodity prices and income avail- 
able. It is derived by maximising the utility subjected to a budgetary con- 
straint. 
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Example 57 derivation of a Marshallian demand function 


Let a utility be 
U= 142 


which is subject to a constraint 


Pig + pag = 6 


where b is the amount of income available, that is to say, our budget. Then 
the Lagrangian function is 


U =q1q2 + A(b— pit — p2q) 
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The first partial derivatives are then 


U1 = q2—- AP =0 (18) 
u2 = qm — Ape =0 (19) 
ux =b—pigi — prg2 = 0 (20) 


where u1, U2 are respectively ug, and Ug- 
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Simultaneously solving Equation’s 18, 19 and 20 leads us to 


Beery a: Mh 
Pr Pp2 


Hence q2 = qipi/p2 and q = g2p2/pi and from Equation 20 we have, 


Mh P2q2 
= pega + pi 
2 Pi 


b=pig + pe 


which yield us, for g, and qz, the Marshallian demand functions which max- 
imise satisfaction of the consumer subject to income and prices. 
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Next, we test the second-order conditions by firstly finding ui, = 0, uaz = 0, 
U12 = Ua = 1, Nn=Pi and 92 = Pp2, which give us 


- 0 1 P1 
|H|=|1 0 pp 
Pi po O 


which gives | H2| = 2pip2 > 0 Hence | H | is negative definite and thus uw is 
maximised. 
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Definition 73 input-output analysis 


The production process of producing one good usually requires the input of 
many other intermediate goods. Let x; be the total demand for product 34, 
and let b be the final demand for the product from the ultimate users. Then, 


ye = OX +...+Gintn +0; 
for? =1,...,n, where a,j is a technical coefficient which represents the value 


of input 7 required to produce one monetary unit’s worth of product j. 
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If we consider the total demand for every one of the products, then 
x= Ax+b 

where 

Z1 Qi -"° Gin by 

In Qn1 ""* Ann bn 
It follows from this that 

x=(F=A)}""b 
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The matrix A is known as the matrix of technical coefficients. It is also 
known as the input-output table, the rows being the inputs and the columns 
the outputs.. The matrix I — A is known as the Leontief matriz. 
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Example 58 complete input-output table 


In a complete input-output table, labour and capital would also be included 
as inputs. These give the value added by the firm. They are normally put 
as an extra row at the bottom of the matrix of technical coefficients A. The 
vertical summation of each column of the table is then equal to 1. 
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Definition 74 eigenvalue and eigenvector 


Let A be a square matrix. Then a scalar \ such that the equation 
Av = Av (21) 


holds for some vector v # 0 is called an eigenvalue } of A, and the vector v is 
called an eigenvector of A corresponding to the eigenvalue 4. The eigenvalue 
is also known as the characteristic root, or the latent root, while the eigenvector 
is also known as the characteristic vector, or the latent vector. 
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Note 2 eigenvalue 


From Equation 21 it follows directly that 
(A-ADv =0 (22) 


Then A—AI is called the characteristic matriz of A. Since v assumes a unique 
value and by assumption v # 0, it follows that A—AJ must be singular, which 
means that its rows must be a multiple of one another. Now A — XI is zero 
if and only if the characteristic determinant | A — AI | of A is zero. 
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In other words 
| A—AI| =0 (23) 


which is called the characteristic equation of A. With Equation 23 there will 
be an infinite solution for v in Equation 22. In particular, if v is a solution, 
that is if it is an eigenvector, so is kv for any k # 0. We force a unique 
solution by using the normalisation 


Soup =1 
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Then the sign-definiteness of A can be determined from the characteristic 
roots X’s. 


Thus if all \’s are positive, then A is positive definite; and if negative, negative 
definite. 


Let at least one be zero, which is neither positive nor negative, if all the 
remaining \’s are nonnegative, then A is positive semidefinite; and if they are 
nonpositive, negative semidefinite. 


Lastly, if some of the \’s are positive while others are negative, then A is 
indefinite. 
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Note 3 
We have seen in Note 2 how, having found »;, where i = 1,...,n, we find 


through normalisation the corresponding, unique v;. On the other hand if 
we have found first the v,’s, their corresponding \;’s may be found by first 
forming a transformation matrix 


T =[v1 ..- Val 


and then the corresponding eigenvalues or the characteristic roots are ob- 
tained from 


1 0 es (0 
TAL = 3: ee < < 
A. dee 
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Definition 75 
The vector equation, Equation 21, has as its solutions the zero vector v = 0 
together with all the corresponding eigenvalue-eigenvector pairs. The set of 
all the eigenvalues of A is called the spectrum of A. The spectral radius of A 
is then the largest of all the absolute values of the eigenvalues of A, that is to 
Say, 


The set of all eigenvectors v;;, together with 0, forms a vector space called 


the eigenspace of A corresponding to j. 
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Definition 76 optimisation 


A problem of optimisation is one in which one tries to maximise or minimise 
a certain quantity called the objective, which depends on a finite number of 
variables. These may be either independent or related to one another through 
some constraints. 
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Definition 77 mathematical programme 


A mathematical programme is an optimisation problem in which the objective 
and the constraints are given as functions or mathematical relationship. In 
other words, 


optimise: z= f(r1,...,2n) 


b;, i=1,...,m 


IVIL IA 


subject to: 9;(X1,.--,2n) J 


Some constraints are explicitly stated as requirements, others are hidden con- 
ditions. These latter need to be pin-pointed through the study and under- 
standing of the model and its inputs. 
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Definition 78 linear programme 


A linear programme is a mathematical programme all the functions involved 
of which are linear. This means that, 


f(@1,---,8n) =av1 t++- +n 2n 
Gi(@1,-+-52n) = a1 +++ + indy 
where 7 = 1,...,m and c; and aj;, j7 = 1,...,n, are constants. If there is 


an additional restriction on the input variables that they be all integers, then 
the optimisation problem is called an integer programme. A mathematical 
programme which is not a linear programme is said to be nonlinear. 
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Definition 79 quadratic programme 


A quadratic programme is a mathematical programme in which all the con- 
straints are linear and the objective function is in quadratic form, which is in 


general, 
n n n 
fiysaj tn) = SOY. cig rea + So dix 
i=1 


i=1 j=l 


where cj; and d; are constants. 
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Definition 80 standard form 
A linear programme is said to be in standard form if all the constraints are 
equalities and if one feasible solution is known. In other words, our problem 


is now T 
optimise: z=c’x 


subject to: Ax =b 
with: x >0 
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Definition 81 initial feasible solution 


One may change any linear programme into the standard form by adding a 
slack variable to the left-hand side of a constraint of the form )7 ajjx; < b; 
to obtain 


n 
) Aiyjlj +Lp, = b; 
j=1 


where p, >n and k =1,2,.... 
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Similarly one may add a surplus variable to the right-hand side of a constraint 
of the form 7 aya; > b; to obtain aia; = bi + oq, 


n 
Daiji — &q = Oi 
j=l 


where q > n and! =1,2,.... Next, all the slack and surplus variables are 
added to the objective function with zero coefficients. 
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Then if we add an artificial variable to the left-hand side of each constraint 
where there is no slack variable, then the initial feasible solution is xo = b, 
where x is the vector of slack and artificial variables. The artificial variables 
are added to the objective function with a large negative coefficient —M. 
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Definition 82 linear dependence 


A set of n vectors of m dimensions {vj,...,Vv,} is said to be linearly dependent 
if there exist some constants a1,...,@, not all of which are zero, such that 
Q1V, +--+: +QnVn =O (24) 


It is said to be linearly independent if the condition in Equation 24 implies 
ay =ag=-:-=a, = 0. 
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Theorem 41 


Consider a set of n vectors of m dimensions. If n > m, then the set is linearly 
dependent. 
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Definition 83 convex combination 


A vector v is called a conver combination of vectors vi,...,Vn if there exist 
some nonnegative constants (1,...,8n, where 


Bite + Bn =1 


such that 
v= Pivit---+Bnvn 
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Definition 84 convex set 


A set of m-dimensional vectors is said to be convez if for any two vectors 
belonging to the set the line segment between them also belongs to the set. 
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Theorem 42 


All points on the line segment joining any two vectors may be expressed as a 
convex combination of the two vectors. 
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Definition 85 extreme point 


A vector v is called an extreme point of a convex set if it can not be expressed 
as a convex combination of two other vectors in the set. 
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In other words, an extreme point of a convex set K is a point x in K that 
cannot be written as « = dy + (1—9)z withO < 6 <1, y and z in K, and 
y # z, that is to say, an extreme point is a point which is not an interior 
point of any line segment belonging to K. 
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An equivalent definition of an extreme point is that x is an extreme point of 
a convex set K if K\ {x} is convex. 
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Definition 86 bounded set 


A metric space is a non-empty set X for which is defined a concept of distance. 
The distance d is called a metric on X, having such properties that, for any 
points x and y in X, we have d(x, y) > 0, and d(x, y) = 0 implies 


r=y 


Furthermore, 
d(x,y) = dy, x) 
and 


d(z,y) < d(x,z) +d(z,y) 
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Let X be a metric space with metric d, let A be a subset of X and let x be 
any point of X. Then the distance from x to A is defined as 


d(x, A) = inf {d(x,a) : a € A} 
whereas the diameter of A is defined as 


d(A) = sup {d(a1, a2) : a1 andag € A} 
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Then a set is said to be bounded if its diameter is finite. 
Further, let zo be a point in X and r a positive real number. Then the open 
sphere S,(xo) with centre xo and radius r is the subset of X defined by 


S,(ao) = {x : d(a, 20) <r} 
A point x in X is called a limit point of A if each open sphere centred on x 


contains at least one point of A different from x. A subset F of X is said to 
be closed if it contains all its limit points. 
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Definition 87 linear space 


A linear space, aka a vector space, is anon-empty set L on which is defined two 
binary processes, say addition and scalar multiplication. Addition is defined 
such that for any z, y and z in L, then x + y is again in L; 


Lry=yrs 


st+ytz)=(a@t+y) +2 


there exists a unique identity element 0, aka zero element or the origin, such 
that x +0 =~2 for every x; and there exists a unique inverse element —2 for 
every x, such that x + (—2x) =0. 
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Scalar multiplication is defined with regard to scalars, some instances of which 
are real and complex numbers, such that for any scalar a and any z and y in 
L, ax is again in L 

a(z+y)=art+ay 


(a+ B)z = axr+ Bu 
(a8)x = a(x) 


and lz = z, where 1 is the identity for scalar multiplication. 
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A normed linear space is a linear space on which is defined a norm, that is 
a function which maps each element x in the space to a real number ||z|| in 
such a manner that ||z|| > 0, and ||«|| = 0 if and only if z = 0 


Ilz + yll < llell + llyll 


and 
llax|] = |a| ||| 
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Theorem 43 


A normed linear space is a metric space. 


Theorem 44 


Any vector in a closed and bounded convex set with a finite number of extreme 
points can be expressed as a convex combinations of the extreme points. 
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Definition 88 definition of the problem 


For two vectors, that is points, s and y in R”, we write x > y if and only 
if 2; > y for all 1 <i <n. A system of m weak linear inequalilties in n 
variables can be written as Ax > b, where A is an m x n matrix. 


A fundamental question concerning such system is whether it is consistent, 
that is to say, whether there exists some x such that Ax = b. 


A system may be inconsistent, or it may have a set of solutions which is 
unbounded. If we sketch our problem on a graph, we may see that it’s solution 
set is conver. 
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Theorem 45 solution space 


The solution space of a set of simultaneous linear equations is a convex set 
the number of extreme points of which is finite. 
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Theorem 46 eztreme-point solution 


Let S be the set of all feasible solutions to the linear programme in standard 
form in Definition 80, in other words, S is the set of all vectors x that satisfy 
Ax =b and x > 0, where A is an m X n matrix. Then S is a convex set, and 
the number of its extreme points is finite. The objective function attains its 
optimum, provided that one exists, at an extreme point of S. If m <n, then 
the extreme points of S have at least n — m zero components. 


Business mathematics, Linear programming, 29°" November 2005 -26— From 5‘ 
November 2005 , as of 14” January, 2007 


216 14 January, 2007 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Algorithm 2 
Procedure for finding basic feasible solutions.: 
Input: Ax = b, Ais an m x n matrix, m <n, rank A =m 
[a, --- a,|]+A 
(aja; +---+2,a, =b) + (Ax =b) 
for i=m+1tondo 
endfor 
(@1,...,%) © solve 1a, +---+2,a, =b 
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Definition 89 simplex method 


The simplex method is a matrix procedure which solves linear programmes of 
the standard form as described in Definition 80 where b > 0. Starting from a 
basic feasible solution xg we locate successively other basic feasible solutions 
giving better values for our objective. For minimisation programmes the 
method uses Table 3, for maximisation programmes the same table is also 
used but with the sign of entries in the bottom row reversed. 
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Table 3 


Table used for minimisation programming in simplex method.: 
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Table 4 Description of the simplex method. 
while negative number exists in d do 


Locate the most negative number in the bottom row of the simplex 
table, excluding the last column. The column in which we find 
this number is called the work column. If more than one such 
column exist, choose one of them. 
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Find the smallest of the ratios between the elements in the last column 
and the elements in the work column of the same row, if these 
latter are positive. The element in the work column that yields 
this smallest ratio is called the pivot element. If there are more 
than one of these, choose one. If none of the elements in the work 
column is positive, the programme has no solution. 
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Using elementary row operations, convert the pivot element to 1 and 
reduce all other elements in the work column to 0. 

Replace the x-variable in the pivot row and first column by the x- 
variable in the first row and pivot column. This new first column 
then becomes the current set of basic variables. 


endwhile 
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The optimal solution is one in which all the basic variables assume the cor- 
responding values in the last column, while the remaining variables are zero. 
The optimal value of the objective function is then the value of the last row 
and last column for a maximisation programme, and the negative of this value 
if the programme is one of minimisation. 
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Algorithm 3 Algorithm for the simplex procedure. 
3-0 
while there exists a negative number in d do 
jejtl 
for i1=1 tondo 
{c’s} < (column no. of the most negative no. in the bottom row) 
(work column) «+ choose one of the {c}’s 
k © (work column) 
endfor 
Ppivot — M 
cH0 
soln + 0 
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for i=1tom do 
if (a;),, > 0 then 


if p< pee then 
rest 
endif 
endif 
endfor 
if soln = 0 then 
no solutions exist 


endif 
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convert} A, such that (a;),,, =1 and (a;),, =0,1<i<m,i€r 
(xo), — Lk 

endwhile 

for i=1to mdo 
(x0); <= (b;); 

endfor 

for i=m-+1tondo 
ri+0 

endfor 

Zee 

if the programme is one of minimisation then 
ze —2* 

endif 
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Definition 90 two-phase method 


The two-phase method is a procedure modified from the simplex method to 
cope with cases when artificial variables exist in the initial solution xo, in 
order to minimise the round-off errors that occur in the calculation. The last 
row in Table 3 in this case is 


d=c’ —cj A=d,+ Mdp 
and consequently we have Table 5 which is used here. Algorithm 3 is then 
firstly applied to the last row, and then again to those elements directly above 


the zeros in that row. 
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When an artificial variable is removed from the first column of the table, it 
ceases to be basic and may be removed from the top row of the table together 
with the entire column under it. When the last row contains only zeros, it 
may be deleted from the table. The programme has no solution if non-zero 
artificial variables are present in the final basic set. 
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Table 3 Table used for minimisation programming using the two-phase 
method. 
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Definition 91 duality in linear programming 
Given a linear programme in the variables x1,...,%, there exists another 
linear programme associated with it, called its dual, which is in the vari- 
ables w1,...,Wm.- The original programme is called the primal. The primal 
completely determines the form of its dual. 
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The symmetric dual of a primal linear programme in the matrix form 


minimise: z=c!x 
subject to: Ax >b 
with: x >0 


is the linear programme 


maximise: z= b’w 
subject to: A’w<c 
with: w>0 
The dual variables w1,...,Wm are known as shadow costs. 


Business mathematics, Linear programming, 29°" November 2005 -41— From 5‘ 
November 2005 , as of 14°" January, 2007 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


The unsymmetric dual of the primal 


minimise: z=c’ x 
subject to: Ax =b 
with: x >0 
is _ 
maximise: z=b’ w 


subject to: A’w<c 
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The unsymmetric dual of the primal 


maximise: z=c x 
subject to: Ax =b 
with: x >0 
is 
minimise: z= b’ w 


subject to: A’w>c 
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Note 4 


We may see from Definition 91 that the dual of a programme in standard 
form is not itself in standard form. These duals are said to be unsymmetric. 
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Theorem 47 


If an optimal solution exists for either the primal or the dual programme, 
then the other programme also has an optimal solution. If the duality is 
symmetric, then the two functions have the same optimal value. If the duality 
if unsymmetric, then the optimal value of each function can be derived from 
that of the other. 
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Definition 92 cut algorithms 
Algorithms which change the boundary of the solution region in order to find 
the optimal solution of an integer programme are called cut algorithms. 


The branch-and-bound algorithm does this by splitting the solution region 
into two and then discard the one which does not contain the optimal solution. 


The Gomory algorithm, on the other hand, reduce the feasible region with 
the help of a new constraint without the region being splitted. 


Business mathematics, Integer programming, 6°" December 2005 -1- From 5°” 
November 2005 , as of 14” January, 2007 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 
Definition 93 branching 
We call branching a process by which a programme whose solution contains 
a non-integral 
J<aj<k 
is made into two separate programmes having the additional constraint 


ai<j 


in one, and 


in the other, the objective together with all the constraints of the original 
problem of which remain the same. Here j and k are positive integers and 
J<k. 
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Definition 94 bounding 


In the branch-and-bound algorithm, if the objective is maximisation, the value 
of the objective obtained when the first integral approximation occurs is said 
to be the lower bound for the problem, and if the objective is minimisation it 
is said to be the upper bound of the same. 
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Algorithm 4 Branch-and-bound algorithm for integer programming. 


find first approximation 
while approximations not all integers do 
choose z; from all non-integral variables such that 
min (|x; — |x4]|, |i — [vi]|) is maximised 
branch 
choose the branch whose value of the objective 
is maximum 
endwhile 
solution¢last approximation 
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Example 59 branch-and-bound example 
(Problem 6.9; Bronson, 1982) 


maximise: z= 2%, +2%24+ 23 
subject to: 22%, + 3%2 + 323 < 11 


with: all variables non-negative and integral 


Solve by branch-and-bound algorithm. 
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Draw a simplex table of Programme 1. 


4 0 
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Replace x4 for x2 as the basic variable. 


x2 
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Since 3 < x5 < 4, branch into two programmes, namely Programme 1 where 
v2 <3, and Programme 2 where x2 > 4. Consider first Programme 2. 


maximise: z= 2%, + 2%2+ 23 
subject to: 2a, + 3%2 + 343 < 11 
t2<3 
with: all variables non-negative and integral 
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Use the simplex method in a tabulated form. 
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Replace the basic variable x5 with x2. 


LA 2 
HD) 3 
6 
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Replace the basic variable x4 with 21. 
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Then consider Programme 3. 


maximise: z= 2; + 2%. +23 
subject to: 2a, + 3%. + 323 < 11 
t>A4 
with: all variables non-negative and integral 
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Draw a table for the two-phase method. 


L4 0 11 
x6 —M 4 
0 

-15 
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Change x4 for x2 in the basic variables. 


X1 Xa x3 v4 Ls x6 
Xo 2 1 1 $ 0 0 + 
x6 —3 0 -1 —; -1 1 z 
1 2 22 
7 0 1 . 0 0 > 
2 1 34 
5 0 1 3 1 -1 =F 
The coefficient parts of the row corresponding to the basic variable xg and 
the last row cancel each other. The optimal result is #3 = oan xi = 23 = 


* pk * _ 22 
v4 = 2; = 2§ =O and 2° = =. 
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(3) to 24, 2% = a (0, =) 


(1) z* = #, (0, 3) 
(2) 25 3, a= if (1,3) 


Therefore the solution is ef = 1, #3 = 3, #3 = xj = wf =0, and z* =7. 


# 
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Algorithm 5 Gomory algorithm for integer programming. 


while solution not wholly all integers do 
choose one non-integral optimal approximation 
write a relation from the row where that variable is basic 
rewrite the relation to make all fractional coefficients 
some integer plus a proper fraction 
move all the fractions to LHS, and all the non-fractions 
to RHS 
write a new constraint as LHS > 0 
find the solution for the original problem together 
with the new constraint 
endwhile 
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Example 60 
(Problem 7.1; Bronson, 1982) 


maximise: z= 2%, +22 
subject to: 2%, + 54% <17 
321 Te 2x9 < 10 


with: 2 ,22non-negative and integral 


Use cut algorithm. 
Solve 


Business mathematics, Integer programming, 6*” December 2005 -17— From 5‘ 
November 2005 , as of 14” January, 2007 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Find the first approximation of Programme 1 normally using the simplex 
method. 


8 
nn 
8 
N 
8 
ro) 
8 
> 
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Since 2 < 42, we know that 3 is the pivot element, and therefore we replace 


the basic variable x4 with x1. 


x3 


Z1 


10 


We have zt = BY, x3 = 4, 23 =x} =Oand 2*=2 
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Since both 2} and 23 are non-integers, arbitrarily choose the former to gen- 
erate a new constraint. Then our Programme 2 becomes, 


ee ee 10a 
+ 5%2+5%4 == = 
Pa ee 8 3 
1 
gree = 3 oi 
24! ithe 
3°2 3°74 3 = 
1 
= ea ee 
gr2t gm a3 
2%2 +44 >1 


and our new programme becomes 
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maximise: z = 2%, + %2 +023 + 024 


: 11 2 31 
subject to: 22 - =%4 = — 


3 3 3 
pee ipiien € 
gee ee 
2%2+%4>1 


with: all variables non-negative and integral 
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Ly v2 x3 v4 5 x6 
1 
2 
3 
1 
3 
2 
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Now 22 replaces xg in the basic variables and becomes the pivot element. 


Oo wiF[NR [Qo 
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This becomes, 


Ly 3 
17 
“3 oy 
1 
tr 3 
13 
2 
Then our first approximation of Programme 2 is af = 3, «3 = 4, x3 = 4, 


T= ef = 0; and z= Sy 
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Arbitrarily choose x4 to generate the new constraint. 


ne i 1 
i: = —z — a 
i a a a 

1 1 1 

= —a —_—_—- 1 4 i 
gu as 2 2 
1 1 1 
ee | 
ae) 
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Then our Programme 3 becomes, 


maximise: z = 2%, + 22+ 0x3 + O04 + O25 


1 
subject to: 21+ gis =3 


Sy) 11 17 


ba Goat oo = 
x Per ies a 
ae; oD 


G4—-%5>1 


with: all variables non-negative and integral 
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ge “i Or 30 0 
vz 0 Pe. (30? 0 oi 30 3 
v2 0 O ae 30 EB Sh OL 4 
zz 0 180 [Or seh le peg 1 
a -M | 0 0 0 1 24 71 1 1 
DA SO 10 Oe oO 0 
OO: 02 al. ae ah ‘et = 
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Then x4 replaces the basic x7 to become the pivot element. 


X41 3 
2 0 
23 11 
t4 1 

0 
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Next, 21 remains basic and becomes a pivot element. 


vy 1 0 0 0 $ 0 3 
v2 0 1 0 0 0 5 0 
Ls 0 0 0 0 -1 -3 8 
4 0 0 0 1 -1 -1 1 

0 -1 0 0 $ 6 
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This becomes 


wy 3 
2 0 
3 8 
L4 1 

6 
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The optimum point for Programme 3 is then, z{] = 3, 23 = 8, zj = 1, 
Ly = 2 = 25 = 0 and z* = 6. 


Therefore the solution to the original problem Programme 1 is 2] = 3, «3 = 0 
at the objective value z* = 6. 


# 
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Definition 95 transportation problem 
A transportation problem involves m sources each of which supplies a;, 1 = 
1,...,m, units of a certain product, and n destinations each of which requires 
b;, i =1,...,n, units of the same. The problem may be stated as following. 
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Mm n 
maximise: z= S- S cigaig 
j=1 j=l 


n 
subject to: S Ly =a, t=1,...,m 
j=l 


m 
yo egeo; j=l1,...,n 
— 


with: all x;; non-negative and integral 


The total supply and the total demand are assumed to be equal. Were this 
not so, a fictitious destination or a fictitious source is added. 
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Definition 96 north-west corner rule 


The north-west corner rule finds an initial basic solution for the transportation 
algorithm of the integer programming. It begins with the (1,1) cell in the 
m x n table, and allocates as many units as possible to 21; violating neither 
the constraints of supply, that is the summation along each row, nor those of 
demand, that is the summation along each column. Then carry on moving 
for each step either right or downwards, until we reach the lower-right corner, 


Zmn- 
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Definition 97 loop 


A loop, which is a sequence of cells in the table used for finding the solution 
in the transportation problem, has the following properties. 
a. each pair of consecutive cells is on either the same row or the same 
column 
b. no three, or in fact any odd-numbered, consecutive cells lie in the 
same row or column 
c. the first and the last cells are on the same row or column 
d. the path along the loop is self-avoiding, that is no cells appear more 
than once in the sequence 
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Algorithm 6 Transportation algorithm. 


while optimal solution not attained do 
find an initial, basic feasible solution using, for instance, 
the North-west corner rule 
let. either uj =0 or vj = 0 depending on whether 
the i**-row or the j**-column 
has the maximum number of basic solutions 
find all u; and vj,i=1,...,mand j =1,...,n 
from u; + vj = cj; for basic variables, and 
from cj — ug — vj; for non-basic variables 
improve the solution 
endwhile 
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Note 5 optimum for transportation problem 


In a transportation problem, optimal solution is achieved when 


Cij — Ui — Uj > 0 


for all transportation costs per unit cj; of all non-basic variables. 
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Definition 98 Present and future values 


The present value po or the principal is the amount initially borrowed or 
invested. The future value p; is the principal after a period of time t. 
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Definition 99 Annual percentage rate 


Interest rates expressed per annum are called nominal rates, i. The annual 
percentage rate or effective annual rate i, is the equivalent annual rate of 
different interest rates variously compounded. 
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Definition 100 Sequence 


A sequence is a list of numbers which follows a definite pattern. It is called an 
arithmetic sequence if each of its terms is obtained from the term immediately 
preceding it by an addition of a constant d, which is called the common 
difference. It is called a geometric sequence if each of its term is obtained 
from the previous term by a multiplication of a constant r, the common 
ratio. 
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Definition 101 Series 


A series is the sum of the terms of sequence. It is called a finite series is 
one whose number of terms is finite, otherwise it is called an infinite series. 
An arithmetic series or arithmetic progression is the sum of the terms of an 
arithmetic sequence. Likewise a geometric series or geometric progression is 
the sum of the terms of geometric sequence. 
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Theorem 48 Arithmetic series 


The value of the n*" term of an arithmetic series is 
T, =a+(n—1)d 
The sum of its first n terms is 


Sn = 5 (2a + (n—1)d) 
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Theorem 49 Geometric series 
The n** term of a geometric series is 
Tar 
The sum of the first n terms of it is 
Sy =atar+-::+ar"} 
a(l—r”) 
=a eee 
a(r” —1) 
Sree aed 
When the number of terms approaches infinity, the summation in cases where 


r <1 becomes is 


~ 1l-r 
Business mathematics, Financial mathematics, 13° December 2005 -6— From 5** 
November 2005 , as of 14°" January, 2007 
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Definition 102 Simple and compound interests 
A simple interest is a fixed percentage of the principal paid to an investor 
each year. A compound interest is an interest paid on the principal plus any 
interest accumulated in previous years. 
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Theorem 50 Present value for simple interest 


The present value for simple interest is 
Pt = po(1 + it) 


where 2 is the interest rate and ¢ the time in years. 


Business mathematics, Financial mathematics, 13° December 2005 -8— From 5‘ 
November 2005 , as of 14” January, 2007 


God’s Ayudhya’s Defence 14 January, 2007 249 


Business Mathematics, notes and projections Kit Tyabandha, PhD 


Kit Tyabandha, PhD Department of Mathematics, Mahidol University 


Theorem 51 Present value, compound interest 
The present value in the case of compound interest is 


pt = po(1 +i)' 
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Note 6 Compounding more than once a year 


The interest may be compounded more than once a year, for example bian- 
nually, quarterly, monthly, weekly, daily, or continuously. Each time period 
is called a conversion period or interest period. The interest rate applied at 
each conversion is ¢/m, where m is the number of conversion periods per year. 
The number of conversion periods over t years is then n = mt. 
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Theorem 52 Present value, compounding several times per year 


The present value at the end of n conversion periods is 


j n j mt 
p= (145) =m (1+ =) 
m m 


where all the variables and parameters are as previously defined. 
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Theorem 53 Present value, continuous compounding 
When the number of compoundings per year becomes very large, the present 
value becomes 
Pt = poe 
Proof. Since p, = po (1+ ay and limmsoo (1+ +)" = e?, we have the 
proof. q 
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Theorem 54 Annual percentage rate 


The annual percentage rate when compounding occurs m times per year is 


j m 
in= (145) —1 


Business mathematics, Financial mathematics, 13°" December 2005 —13- From 5** 
November 2005 , as of 14*” January, 2007 


252 14 January, 2007 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD 


God’s Ayudhya’s Defence 


Date 
25 Oct 2005 
1 Nov 2005 
8 Nov 2005 
15 Nov 2005 
22 Nov 2005 
29 Nov 2005 
6 Dec 2005 
20 Dec 2005 
10 Jan 2006 
17 Jan 2006 
24 Jan 2006 
31 Jan 2006 
7 Feb 2006 
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Appendix 


Course Outline 


Topic of lecture 
Graph and derivative 
Calculus of multivariable functions 
Exponential, log and nonlinear functions 
Matrix 
Linear algebra 
Linear programming 
Integer programming 
Financial mathematics 
Integral calculus 
Integral calculus 
Simultaneous equations 
Differential equation 
Trigonometric functions and power series 
Differential equation 


Materials 


14 January, 2007 
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S 
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Content Date begun Dated Last updated 
Graph and derivative 20 Oct 05 25 Oct 05 10 Dec 05 
Calculus of multivariable 

functions 20 Oct 05 1 Nov 05 18 Nov 05 


Exponential-, logarithmic 
and nonlinear functions 28 Oct 05 8 Nov 05 18 Nov 05 


Matrix 5 Nov 05 15 Nov 05 10 Dec 05 
Linear algebra 5 Nov 05 22 Nov 05 10 Dec 05 
Examples for 

linear algebra 17 Jan 06 31 Jan 06 31 Jan 06 
Exercises for 

linear algebra 7 Feb 06 7 Feb 06 7 Feb 06 
Linear programming 5 Nov 05 29 Nov 05 10 Dec 05 


Examples for 

linear programming Nov 05 6 Dec 06 10 Dec 06 
Integer programming 5 Nov 05 6 Dec 05 11 Dec 05 
Financial mathematics 5 Nov 05 13 Dec 05 20 Dec 05 
Examples for 

financial mathematics 18 Feb 06 20 Feb 06 20 Feb 06 
Simultaneous equations 5 Nov 05 24 Jan 06 26 Jan 06 


or 


Differential equation 5 Nov 05 31 Jan 06 31 Jan 06 
Integral calculus 5 Nov 05 10 Jan 06 13 Feb 06 
Integral calculus 6 Feb 06 7 Feb 06 7 Feb 06 
Examples for 

integral calculus 17 Jan 06 7 Feb 06 7 Feb 06 
Exercises for 

integral calculus 10 Jan 06 7 Feb 06 7 Feb 06 
Difference equation 5 Nov 05 20 Feb 06 19 Feb 06 
Examples for 

difference equation 19 Feb 06 20 Feb 06 20 Feb 06 

Appendix 

Midterm examination 19 Dec 05 20 Dec 06 20 Dec 06 

[solution] 5 Jan 06 
Quiz 1 30 Jan 06 31 Jan 06 13 Feb 06 

[solution] 16 Feb 06 
Quiz 2 7 Feb 06 7 Feb 06 7 Feb 06 

[solution] 23 Feb 06 
Quiz 3 13 Feb 06 14 Feb 06 13 Feb 06 

[solution] 14 Feb 06 
Final examination 19 Feb 06 =.20 Feb 06 20 Feb 06 

[solution] 22 Feb 06 


Teaching method 


There were lectures and practices in class. All practice exercises, quizzes 
and exams were done in an opened-book manner. Practice of problem and 
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exercise in class were preferred to assignment and homework, since in the 
latter many students did not do the work themselves. In some of these practice 
the questions for the students were all different in order that they may start 
to think for themselves. 


Teaching media 


A camera projector and a microphone were the hardware media used. The 
material used as a media were lecture projections in this collection and lecture 
hand-outs, which will go into another collection. Because of the distance be- 
tween the class and its lecturer, feedbacks of homework were sent to students 
by email. Textbooks on Calculus from the library were used in some of the 
exercise sessions, for the students to learn how to use books as a resource to 
help them solve problems and answer questions. 


Evaluation methods 


means per cent 
Attendance 10 
Homework 10 
Quiz 1 10 
Quiz 2 10 
Quiz 3 10 
Midterm exam 20 
Final exam 30 


Evaluation will be based on the distribution plot, as relative performance 
among students. 
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Quiz 1 
Business Mathematics 
31%* January 2006 
Time: 1 hour (10-11pm) 


1. Let 
3 141 
A=|{2 4 2 
1 1 3 


Find all the eigenvalues and a basis for each eigenspace. If possible, find 
invertible matrices P such that P~! AP is diagonal. 


Solution. Form the characteristic matrix, 


#-3 -1 —1 
tfl-A=| -2 t-4 -2 
—1 -1 ¢-2 


> a 
= a) 
Sj oh “f= 3 
(¢ — 3) ((f— 4)(f -— 3) — 2) + ( — 8)(—2) — 2) -(2+ 4-4) 
= t? — 10¢7 + 28¢ — 24 
= (t — 2) (t? —8t +12) 


& 
| 
- 
I 


= (t— 2)°(¢-6) 
Therefore the eigenvalues are 2 and 6. 
# 
For the eigenvalue 2; 
2-3 -l1 —1 xz -1 -1 -1l My 0 
= ee y|={-2 -2 -2]}/y]={0 
—1 -1 2-3 z -1 -1 -1l z 0 


which gives x+y+z= 0. The system has thus two free variables, for example 
{z =1,y=0,z =—1} and {* =1,y =—-1,z=0}. Let wu = (1,0,-1) and 
v = (1,-1,0). Then all linear combinations of u and v forms a set of all 
possible eigenvectors, which together with zero vector comprise the eigenspace 
corresponding to the eigenvalue 2. In other words, u and v form a basis of 
the eigenspace of the eigenvalue 2. 


# 
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For the eigenvalue 6; 


6-3 -1l —1 x 3 -1 -1 x 0 
—2 6-4 —-2 yJ={[—-2 2 -2 y |= {0 
—1 -1 6-83 z -1 -1 3 z 0 


Use Gaussian elimination to find the solution of this system of equation. First 
form an augmented matrix 


ce Earns 1 
BD: Oo 29 OG 
=f oi 3-0 
Qe De =i =1 tee Cr <1 3 0); 
Sh) a> 9; 
A Be 9 2G 
oe ee 
—1(J); 
i ob. <3. 
ED DO 220 
= a ee a 


(aa ( 2 2 Se DASA. TE =a) 9 
i O)= 3" a 238° <0); 


ae oe a 
0. As = 8°50 
4 28.0 

1(I1),(1 —2 0); (111) +4(D, (-4 -8 0)4+401 -2 0); 
ee 
Oo Ae 280 
O50 40.40 


Therefore rank of A is 2. The system has only one free variable, that is to 
say, one independent solution. Any particular non-zero solution generates its 
solution space, that is the eigenspace, for example x = 1, y = 2 and z = 1. 
So w = (1,2, 1) forms a basis of the eigenspace of the eigenvalue 6. 


# 
1 1. ot 

Let Pbe | 0 -—1 2] Since det P = 1(-1) —1(2+1) = —4 £0, hence 
-1 0 1 


P has an inverse. Therefore P is the required matrix such that P~1!AP is 
diagonal. 


# 


Bibliography 
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Quiz 2 
Business Mathematics 
7h February 2006 
Time: 1 hour (10-11pm) 


1. Evaluate the following integrals. [10] 


1 2 14 2 1 
/ wdx, i: (x? — 3) da, / -dz, : (x? + Jaz) dz, i ve dz 
) 0 -14 ) 0 


Solution. Z 
[ dz = n(2\g =7 
0 


2 2 
i. (x? — 3) dz = (S- 2 
0 
# 


In Ae Ldz the limits of integration pass through 0, where there is disconti- 
nuity and where 1/z is undefined. 


Further, lim,_,9- f(z) = —oo and lim,_,9+ f(z) = 00. We look at 


8 
ee a 
3 


0 


1 
1 
/ —dz = lim (Inz|; = lim (Inl —Ina) =0—-(—o) = 
0 2 aot aot 


1 0 1 
1 

[ se-/ wde+ [ Le 

1% 1% 0 % 


0 b 
1 ; 1 . b ; 
SU ——d(-z) = 1 In(- = | In(1 
de ut im. ie ai z) lim. (n( z)\_4 jim (nb n(1)) 


and 


where 


lim (1 
ee 


the latter of which is undefined. Therefore ie +dz is also undefined. 


# 
2 ee ee ee ee 
2 pee | EG egy A | ca ee DS 
[ + va ae gtae Lat ; 
# 
Let u= «x and dv = e*dz. Then du = dz and v = e”. And then, 
1 1 
[ cetaz= (eer - f e*dx = (xe* — e*|5 = 1 
0 0 
# 
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Quiz 3 
Business Mathematics 
14% February 2006 
Time: 1 hour (10-11pm) 


Choose only one problem[10], either 


1. Solve 
y+3z=9 
22+ 2y—z=8 
—£+5z=8 


by Gaussian elimination. 


Solution. Write an augmented matrix, 


0 1 3 9 
2 2 -1 8 
-1 0 5 8 


(I) + (IIT), (0 1 3 9)H(-1 0 5 8); 


-1 0 5 8 
2 2 -1 8 
0 1 3 9 


SAT aa 0 8 Sp Ho Oe 2.1. sya 20 0 a5. 8h 


1 0 -5 -8 
0 2 9 24 
01 3 9 


II 4 III, (0 2 9 24,4 (0 1 3 9); (III)- 21D), 2 9 24)- 
2(1 3 9); 

1 0 -5 -8 

0 1 8 9 

0 0 3 6 
Therefore, directly we have z = 2, y = 9 — 3(2) = 3 and = —8+ 5(2) =2. 


or 


2. Solve 
Qe +y — 22 = 10 


32+ 2y+2z=1 
5a + 4y+32=4 
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by any method. 


Solution. Form an augmented matrix, 


21 -2 
3.2 2 
; 4 3 
3(2), (1 § -1 5); (IZ)-37), (3 2 
51), (5 4 3 4)-5(1 $ -1 5); 


2(II), (1 10 —28); (III) — 3(z1), (3 
oat 
0 1 10 
0 0 -7 

Directly, z -2 —3,y 
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2. 123 <b Sa, SE) = 
5 

-14 

I 

8 -—21)-8(1 10 —28); 
5 

~28 

21 


—28 — 10(—3) = 2 and « = 5— 4(2) + (-3) =1. 


# 
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Midterm examination 
Business mathematics 


20%" December 2005 
Time: 3 hours (14pm) 


Te 
. Given 
f(e) =az? +br +e 


where a = c= 1and b = 2. Find all the x-intercepts, y-intercepts, and critical 


points. 
. Given 

Iny=1-¢2 
Find y. 


. Given 


Find b. 


. Given 
logy = 3loga + logb — loge 


Find y. 


Solution. 


f(x) =2? +2241 


x-intercept, y = 0; 
x? +2¢+1=0 


(e« +1)\(e¢+1)=0 


z=-l 
f(x) touches x-axis at ¢ = —1 
Critical point; f'(x) = 0; 
227 +2=0 
cr=-l 
y-intercept, x = 0; 
y=l 
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y= ei-@ 
# 
gi@ty—2z)n _ a? 
b=(a4+y—2z)n 
# 
loga® + log b — loge = logy 
ab 
log — =] 
0g - ogy 
a®b 
¢ 
# 
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q=lnz+lIny 


g=eter tet 2ryt ay? 


p = 1500974 
Find the first- and second-order partial derivatives. 
[10] 
Solution. 
OG ed. 
Ox 
# 
he 
dy sy 
# 
a’*q_ 1 
Oa2 sg? 
# 
a?q _ 
Oxdy 
# 
Org. A 
dy??? 
# 
2 
ce 0 
OyOxr 
# 
ic ee early eee ery 
Ox 
# 
Oz 
—=2 2 
By z+ 2xry 
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# 
Oe -. 
0x2 
# 
072 
— =2+2 
Oxdy ee 
# 
Oy? 
# 
072 
Oydx Ae ey 
# 
Op 0.74t 0.74t 
Fe = 150(0.74)e°* = 1116" 
# 
2 
vp = 111(0.74)e°7* = 82.1469 
# 


14 January, 2007 


God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


3. Find the determinant of the following matrices. 


a. 
a b 
c d 
b. 
a be 
def 
g h i 
C. 
a bed 
e f gh 
ij kil 
mn o p 
Solution. 
a. 
a bj _ 
oh) =ad—be 
b. 
a be 
de fl=alf Z]-s]é {]+e|¢ 
g h i g 


= a(ci — fh) — b(di— fg) + c(dh 


g h 


= aci — afh — bdi+ bfg+ cdh 
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a eas f gh e gh 
ae a ee sie eae ae ee 
Be ed n oOo p mo p 
mn o p 
e fh e fg 
+e)/t g l|-dl/i gj k 
mn p mn o 


= a(f(kp — lo) — g(jp — In) + h(jo — kn)) 
— b(e(kp — lo) — g(ip — lm) + h(io — km)) 
+ c(e(jp — In) — f(ip — Im) + hin — jm)) 
— d(e(jo— kn) — f(io — km) + g(in — jm)) 
= afkp — aflo — agjp+agin+ ahjo— ahkn 
— bekp + belo + bgip — bglm — bhio + bhkm 
+ cejp —celn—cfip+cflm+ chin — chjm 
— dejo + dekn + dfio — dfkm — dgin + dgjm 
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4. The relationship between the total revenue r;, the price p, and the output 


quantity q is 


Tt = pq 


The demand function is p = a — bq, where a and b are positive constants. 


Find r;, the marginal revenue r,,, and the average revenue rg. Then find q 
at the maximum r;. And then sketch the graphs of rz, rm and ra. (assume 


az > 4b) 


Solution. 


ry = pq = (a — bq)q = aq — bg? 


Pe 
t= — =p=a-by 
q 


At maximum 1; 


r, =a—2bg =0 
a 


I~ oD 
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Tt; Tm and Tq (a2 > 4b,a=5, b= 1) 


“4 0 1 2 3 4 5 6 
qd 
Figure 14 
=4, n=a(4)- (oie Ny pe pa 
oe Oph FP NO 4p?) 48 2b) 4b 
q=0; rm =0 
Te = 0; 
g(a — bq) = 0 
a 
=-0. = 
P= 0) 5 
q=0; Tm=a 
a 
Tm = 0; aT 
q=0; rm=a 
a 
a = 0; = 
r qd b 


This can be summarised as Figure 14. 
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5. 


a. Let A co B =) and C SS represent the second-order condi- 


tions of critical point of function. Suppose the graph of a function has the 
shape as shown in Figure 15. 


Figure 15 
Which of these conditions is satisfied at points I, P and T?[3] Explain.|2] 


(e) = 32° ao +a? + sa? 6a +7 
Find f’(#) and f"(x).[1] Show that 1, —1 and 2 are the critical points.[2] 


Which of these are maximum, minimum or inflection point?[2] 


Solution. 


I ¢B 
P <¢<C 


T <A 


At I, P and T, f'(-) =0. 


At I, PandT, f"(-)=0, <0, and > 0 respectively. 


# 
f(z) = 32+ — 92? + 327 + 9x * 
f'(e)  =12x3 — 2727 +62 +9 
# 
f'(x) = (32? — 6x + 3) (w + 1)(a — 2) = 3(@ — 1)?(@ +:1)(w — 2) 
Critical points, 1, -1 and 2. At these points f’(-) = 0. 
# 
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f"(-) =0 +1 inflection point 
f"(-1) <0 4-1 maximum point 
f"(2)>0 42 minimum point 
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6. Solve the following programme by the simplex method. 


maximise: z= 321 + 4%2 + 523 
subject to: 2, +2%2+23 <2 
@+%.4+323 <1 
321 + 222 +23 <4 
with: all the variables non-negative 


[10] 


Solution. Draw simplex tables. 


X11 x2 x3 XL4 x5 2X6 
2 2 1 5 
v4 3 3 0 1 ae 0 3 
1 1 1 1 
2X3 a = 1 0 3 0 3 
2X6 § 3 0 0 = 1 oh 
0 0 0 
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Final Examination 
Business Mathematics 


20% February 2006 
Time: 3 hours (9-12am) 


1. 
. Give the definition of a homogeneous function. [2] 
. What are order and degree of a differential equation? [2] Given the differential 
equation 
(iv) ny2 ny3 _ gs 
yee + 4(y") + yy” = sina 
Give its degree and order. [2] Is this equation homogeneous? [1] 


. A primitive gives rise to a differential equation. Explain how the primitive 
3 
y = Ax? + Bx + C gives rise to the differential equation £$ = 0. [3] 


Solution. 
. A function f(x,y) is said to be homogeneous of degree n if f(A, Ay) = 
A” f(a, y)- 

# 


# 
The degree of a differential equation is the degree of its highest ordered deriva- 
tive. 


. The order of a differential equation is the order of its highest derivative. 


# 
The degree is 1, and the order 4. 
# 
The equation is not homogeneous. 
# 
dy =2Ar+B 
dz 
d’y 
—> =2A 
dx? 
d®y 
dx? 
# 


2. 

. Explain what sequence and series are. When are they said to be arithmetic?, 
when geometric? [6] 

. Given an arithmetic series, 1+3+5+7+---. What is the value of the n*® 
term? Find the value of the 50‘ term. What is the sum of the first n term? 
Find the sum of the first 50 terms. [4] 
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. Given a geometric series, 1 + s + ; + z +-++. What is the value of the n*® 
term? Find the value of the 50* term. What is the sum of the first n terms? 
Find the sum of the first 50 terms. What is the sum to infinity, that is to say, 
Soo? [5] 


Solution. 
. A sequence is a list of number following a definite pattern. A series is the sum 
of the terms of sequence. They are said to be arithmetic if each of their terms 
is obtained from the term immediately preceding it by an addition of some 
constant. They are said to be geometric if each of their terms is obtained 
from the previous term by a multiplication of some constant. 

# 
. From what is given, a = 1 and d = 2. The value of the n*” term is then T,, = 
at+(n—1)d = 1+(n—1)2d, while that of the 50° term Tso = 1+(50—1)2 = 99 
The sum of the first n terms is S, = } (2a+(n—1)d) = $(2+(n—1)2), 
and the sum of the first 50 terms is S59 = 32 (2(1) + (50 — 1)2) = 2500. 


# 
. From the given geometric series, a = 1 and r = 4. Then JT, = ar™ + = 
n—1 50-1 a(1—r” 1~(3)" 
(1)"*, Tan = 2)" = gles So = aD = HY <2 (1 2), So = 
rae eee and Soo = a 
# 


3. 

. Give the formula for finding the present value of compound interest. Let the 
annual interest rate be 10 per cent. What is the present value of 1,000 Bahts 
in ten years’ time? [4] 

. What is the present value at the end of n conversion periods in ¢t years? Here 
m is the number of conversion periods per year. Find the present value of 
1,000 Bahts in five years’ time, when m = 5 and i = 10 per cent. [4] 

. Compare the present values obtained from (a) and (b), and discuss the dif- 
ference between them. [2] 


Solution. 
. The formula for the present value of compound interest is p; = po(1+i)'. From 
what is given, i = 0.1, pp = 1000 andt = 5. Therefore, ps = 1000(1+0.1)° = 
1610.51. 

# 


. The present value of compound interest at the end of n = mé# conversion 
period is py = po (1 + ce From what is given, i = 0.1, m = 5, po = 1000 
and t = 5. Therefore ps = 1000 (1 + %+) = 1485.95 Bahts. 


- Pt > Ptym 
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4. 
. What is the average or mean value of an integrable function f(x) on [a, b]? 


[2] 

. What is an integrable function? [2] Give an example of a function that is not 
integrable on some range. [2] 

. Consider a (2 x 2) system of linear equations in the slope-intercept form, 


y= myx + by 
y = Mou + be 


Give conditions for the system to have a unique solution, no solutions, and 
infinitely many solutions. [4] 


Solution. 
. The average or mean value is 


b 
# 


. An integrable function is a function whose value never becomes infinitely large 
on the range of interest. 


# 
# 


. The system has a unique solution when m, #4 mz, no solutions when b; # be, 
and an infinite number of solutions when m, = mz and b, = be. 


The function f(x) = + is not integrable over the range [—1, 1]. 


zx 


# 


5. 
. Explain revenue and elasticity mathematically and verbally. [1] Explain 
marginal and average costs both mathematically and also in words. [1] 


b. Give the condition for a° = 1 to be true. [1] 


. The following is a formal definition of limit. 
For a function f(x), lim,_,, f(x) = 1 if and only if for every € > 0, there 
exists 6 > 0 such that |f(a) —1| < ¢ whenever 0 < |z —a| <6. 


Explain in words what we mean by this. Give your answer both in English 
and in Thai. [1] 

. What is a cut algorithm? [1] In the following solved example, identify all pivot 
elements and explain in detail how the problem was solved. [5] 


maximise: z= 27%, + x2 
subject to: 2%; + 5%. < 17 
321 + 22 < 10 


with: 21,22 non-negative and integral 
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Use cut algorithm. 
Solution 


Find the first approximation of Programme 1 normally using the simplex 
method. 


L3 0 17 
X4 0 10 
0 


Since a < i, we know that 3 is the pivot element, and therefore we replace 
the basic variable x4 with 71. 


Ly r2 L3 L4 
=m | 0 # 1 -# | & 
my 1 2 0 3 ig 
a a ae = 


We have sj = 2, 3 = 4, 3 = aj = 0 and z* = 2. Since both xj and 23 


are non-integers, arbitrarily choose the former to generate a new constraint. 
Then our Programme 2 becomes, 


gee gel W250! 
x =%24+5%4=— = = 
ee a eee 3 
1 
a a a 
aa a 
Ay >@4— > 
go ee 
ree 
—X:. —2x pi 
igi Mi 
2%. +24>1 


and our new programme becomes 


maximise: z = 2%, +2%2.+ 0234 024 
11 2 31 


bject to: a2 — <a, = = 
subject to 3 02 3°4 3 
4s 1 10 
t+ =% — -44 = — 
sues Aieaaae ame 
289 +24 >1 


with: all variables non-negative and integral 
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SO wiHl[wle wR ow 


we NIE x5 ww 


Then our first approximation of Programme 2 is z{ = 3, 23 S, x3 ur 


zy = xf = 0, and z* = 2. Arbitrarily choose #3 to generate the new 
constraint. 
ree 

2 gu gs 2 

oy es oe 

hie oe a a 

1 1 1 a 

aan ae 
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Then our Programme 3 becomes, 


maximise: 


subject to: 


with: 


2 = 24, + 224+ 0x3 + 0x4 + O25 


1 
fit teS 3 
Pe re ae 
3 Q°4 6 nD) 
42 | 
LQ g%4 g75 5 
G4—-%5>1 


all variables non-negative and integral 


We draw our table for this programme. 


This becomes 


X2 x3 L4 Xs x6 L7 | 
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X41 3 
X2 0 
£3 8 
4 1 

6 


The optimum point for Programme 3 is then, z] = 3, 23 = 8, 24 = 1, 
x5 = xf = x2§ = 0 and z* = 6. Therefore the solution to the original problem 
Programme 1 is z] = 3, x5 = 0 at the objective value z* = 6. 


# 


278 14 January, 2007 God’s Ayudhya’s Defence 


Kit Tyabandha, PhD Business Mathematics, notes and projections 


Students’ scores 
Business Mathematics 
2005-6 


Kit Tyabandha, PhD 
14% January, 2007 


Introduction 


This is a report of study and examination results for the course Business 
Mathematics taught by me during the second semester, 2005-6 academic 
year. It has been my concern from the beginning of 2006 that students had 
not done their homework themselves, therefore I have had since then practices 
and quizzes in class, with the hope that they would obtain some problem- 
solving skill. Everything, which includes exams, quizzes and practices, was 
done in an opened-book manner because that is how the real world works and 
that is where the progress lies. 


Midterm exam 


There were six questions in the midterm exam, namely from M1 to M6. MT 
is the total marks and MM the actual contribution toward the students’ final 
score. Table 6 shows the six marks, M1—M6, for each student ID, together 
with the total mark and the same in per cent for each. 


ID M1 M2 M3 M4 M5 M6 Total (60) Per cent (20) 


4661 0.5 83 5.2 0 2.2 1.5 17.7 5.9 
4801 6 7.5 5.2 0 2.7 3 24.4 8.1 
4802 7 64 5.2 0 2 3 23.6 7.9 
4803 65 76 5.2 35 3 1.5 27.3 9.1 
4804 

4805 0 5) 0.2 0 0 4.7 14.9 5) 
4806 7 3.3 5.2 0 1 2 18.5 6.2 
4807 7 0.8 52 1 1 0 15 5) 
4808 1 4.2 5.2 1 1 0 12.4 4.1 
4809 1 5) 0.2 1 1 2 15.2 5.1 
4810 1.5 5 0.2 3.5 1 0 16.2 5.4 


Table 6 Midterm marks, Business Mathematics, 2005-6, contribution of 
20 per cent toward the overall points 
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ID 

4811 
4812 
4813 
4814 
4815 
4816 
4817 
4818 
4819 
4820 
4821 
4822 
4823 
4824 
4825 
4826 
4827 
4828 
4829 
4830 
4831 
4832 
4833 
4834 
4835 
4836 
4837 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 
4848 
4849 
4850 
4851 
4852 


Table 6 (continued) Midterm marks. 
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M2 
7.6 


M3 
5.2 


5.2 
5.2 


M4 
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Total (60) Per cent (20) 


19.8 6.6 
21.3 7.1 
32.4 10.8 
25.3 8.4 
10.9 3.6 
24.7 8.2 
16.8 5.3 
10.1 3.4 
12.9 4.3 
12.8 4.3 
30.1 10 
29.6 9.9 
20.1 6.7 
32.3 10.8 
3.2 10.8 
30.2 10.1 
11.3 3.8 
13.9 4.6 
12 4 
13 4.3 
12.8 4.3 
14.5 4.8 
30.4 10.1 
31.1 10.4 
31.6 10.5 
33.1 11 
38.1 12.7 
26.6 8.9 
37.6 12.5 
11.5 3.8 
17.3 5.8 
18.7 6.2 
32.3 10.8 
36 12 
31.4 10.5 
25.6 8.5 
25.3 8.4 
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ID 

4853 
4854 
4855 
4856 
4857 
4858 
4859 
4860 
4861 
4862 
4863 
4864 
4865 
4866 
4867 
4868 
4869 
4870 
4871 
4872 


M1 


M2 
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M4 


M5 


M6 


1.5 
1.2 
1.5 
1.5 


Table 6 (continued) Midterm marks. 


Table 7 shows mark and rank of each student. 
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Total (60) 


22.6 
28.4 
25 
28 
30.1 
28.8 


26.3 
23.6 
22.7 


26.3 
24.3 
21.8 
21.4 


26.3 
31.1 
29.2 


Per cent (20) 
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ID 


4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 


661 
801 
802 
803 
805 
806 
807 
808 
809 
810 
811 
813 
814 
815 
816 
817 


ID 


KRR RR RRB BRR BRR BRB A 


818 
819 
820 
821 
822 
823 
824 
825 
826 
827 
829 
830 
831 
832 
833 
835 


836 
837 
838 
839 
840 
84 
84 
84 
84 
84 
84 
84 
84 
850 
851 
852 


onNraoao hk WN eS 


Table 7 Mark and rank of students’ midterm scores. 


ID 
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854 
855 
856 
857 
858 
859 
861 
862 
863 
865 
866 
867 
868 
870 
871 
872 
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Score distribution, Business Maths, Midterm 2005-6 
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Figure 16 Distribution of students’ midterm scores. 


The mean score for the midterm exam is 7.78, median 8.18, minimum 3.37 
and maximum 12.7. The standard deviation is 2.53. 
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Quiz 1 score distribution, Business Maths, 2005-6 
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number of students 
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Figure 17 Distribution of students’ scores from Quiz 1. 
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ID Score Rank ID Score Rank ID Score Rank ID Score Rank 


4661 5.2 18 4818 8 3 4836 5.7 15 4854 8.2 2 
4801 8.6 1 4819 4.2 22 4837 4 23 4855 7.8 5 
4802 6.1 13 4820 4 23 4838 7.2 9 4856 8.2 2 
4803 8.2 2 4821 4 23 4839 3.7 24 4857 8.2 2 
4805 6 14 4822 7.9 4 4840 5.2 18 4858 7 10 
4806 6 14 4823 7.7 6 4841 6.2 12 4859 8.2 2 
4807 7.9 4 4824 8.2 2 4842 5.4 17 4861 6.9 11 
4808 6.9 11 4825 5 20 4843 8.2 2 4862 7.6 7 
4809 3.4 25 4826 8.2 2 4845 4.2 22 4863 8.2 2 
4810 8.2 2 4827 8.2 2 4846 7.2 9 4865 5.1 19 
4811 7.2 9 4829 8.2 2 4847 7 10 4866 4.2 22 
4813 8.2 2 4830 4.2 22 4848 8.2 2 4867 5.2 18 
4814 8.2 2 4831 8.2 2 4849 7.7 6 4868 7.3 8 
4815 8.2 2 4832 6 14 4850 8.2 2 4870 7.7 6 
4816 6.2 12 4833 4 23 4851 5.1 19 4871 4.7 21 
4817 4 23 4835 7.8 5 4852 5.5 16 4872 8.2 2 


Table 8 Quiz II marks, Business Mathematics, 2005-6, contribution of 10 per cent toward the 
overall points 
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Quiz 2 score distribution, Business Maths, 2005-6 
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Figure 18 Distribution of students’ quiz 2’s scores. 
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Quiz 3 
We had our Quiz 3 on 14 February 2006, Valentine Day. From the results, 


the mean is 8.15, median 10, minimum 2, maximum 10, and the standard 
deviation 3.14. 
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Quiz 3 score distribution, Business Maths, 2005-6 
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Figure 19 Distribution of students’ quiz 3’s scores. 
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Final Exam 


The final examination was held on 20 February 2006. There were five ques- 
tions, with marks respectively from the first to the last, 10, 15, 10, 10 and 10, 
totalling 55. 


ID Q1 Q2 Q3 Q4 Q5 Total (55) Percent Scaled (30) 
4661 2 4.9 3 0 0.5 104 18.91 5.67 
4801 5 0 0 4 2.8 11.8 21.45 6.44 
4802 2 6.4 2 4 3 17.4 31.64 9.49 
4803 2 0 11 4 1 8.1 14.73 4.42 
4805 4 6.5 0 6.1 4 20.6 37.45 11.24 
4806 1 12.5 0 3 2 18.5 33.64 10.09 
4807 5 6.2 0 2.1 58 19.1 34.73 10.42 
4808 3 0 0 0 18 4.8 8.73 2.62 
4809 2 6.4 0 0.1 O 8.5 15.45 4.64 
4810 3 10.4 0 3.1 2 18.5 33.64 10.09 
4811 0 1.7 0 41 1 6.8 12.36 3.71 
4813 2 6 0 0 2 10 18.18 5.45 
4814 5 6.5 0 4 2.8 18.3 33.27 9.98 
4815 1 6.9 0 2.1 3 13 23.64 7.09 
4816 3 6.7 0 0.1 O 9.8 17.82 5.35 
4817 2 0 2 41 0.2 8.3 15.09 4.53 
4818 3 7 2.2 95 3 20.2 36.73 11.02 
4819 2 0 0 4 0 6 10.91 3.27 
4820 0 0.2 0 0 1 1.2 2.18 0.65 
4821 2 0 0 0 6.8 8.8 16 4.8 
4822 0 10 3.1 4 8 25.1 45.64 13.69 
4823 1 9.3 0 1 1.5 12.8 23.27 6.98 
4824 5 8.3 41 61 63 29.8 54.18 16.25 
4825 4 6 2.1 4 0 16.1 29.27 8.78 
4826 2 0.2 0 6.1 3.5 11.8 21.45 6.44 
4827 1 0 0 4 2 7 12.73 3.82 
4829 0 10.5 0 0 2 12.5 22.73 6.82 
4830 0 2 0 0 0 2 3.64 1.09 
4831 0 13.8 9 1 0 23.8 43.27 12.98 
4832 2 0 0 4 0 6 10.91 3.27 
4833 0 6.5 2 6 2 16.5 30 9 


Table 10 Final marks (continued) 
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ID Q1 Q2 Q3 Q4 Q5~ Total (55) Percent Scaled (30) 


4835 1 6.5 2 4 0 13.5 24.55 7.36 
4836 3 0 41 51 3.3 15.5 28.18 8.45 
4837 1 2 2.2 4 a 10.2 18.55 5.06 
4838 0 0.2 0 0 1 1.2 2.18 0.65 
4839 5 12.5 1.1 2.1 5 25.7 46.73 14.02 
4840 3 12.9 4 3.1 2 25 45.45 13.64 
4841 4 8.8 0 4 3.8 20.6 37.45 11.24 
4842 6 13.9 3.2 4 0 27.1 49.27 14.78 
4843 3 13.5 62 4 6.8 33.5 60.91 18.27 
4844 7 1.7 0 0 1 9.7 17.64 5.29 
4845 2 5.5 24 0.1 0 10 18.18 5.45 
4846 0 0 0 0.1 0 0.1 0.18 0.05 
4847 0 0.5 0 0 0 0.5 0.91 0.27 
4848 1 10.1 1.1 3 3 18.2 33.09 9.93 
4849 2 10.1 2.2 51 48 24.2 44 13.2 
4850 3 7.3 0 4 5) 19.3 35.09 10.53 


Table 10 Final examination marks, Business Mathematics, 2005-6, con- 
tributing 30 per cent towards overall points 


ID Q1 Q2 Q3 Q4 Q5- Total (55) Percent Scaled (30) 


4851 2 0 0 4 0 6 10.91 3.27 
4852 2 0 0 4 4 10 18.18 5.45 
4854 6 0 0 6 3 15 27.27 8.18 
4855 2.2 8 1 4 4 19.2 34.91 10.47 
4856 6 0 4 4 0 14 25.45 7.64 
4857 4.1 2.2 10 «6 2 24.3 44.18 13.25 
4858 0 6.8 9 6.1 2 23.9 43.45 13.04 
4859 3 6.2 0 6.1 0 15.3 27.82 8.35 
4860 4.5 13.5 3.2 41 0 25.3 46 13.8 
4861 2 6.1 0 4 3 15.1 27.45 8.24 
4862 2 9.9 0 4 4.8 20.7 37.64 11.29 
4863 3 126 28 7 3.9 29.3 53.27 15.98 
4865 4 0 0 4 2.8 10.8 19.64 5.89 
4866 3 3 2 0 0 8 14.55 4.36 
4867 0 0 0 4 0 4 7.27 2.18 
4868 0 0.2 1 0 0 1.2 2.18 0.65 
4870 0 4.2 0 0 0 4.2 7.64 2.29 
4871 4 8.7 3.1 4 5) 24.8 45.09 13.53 
4872 1 0 0 7 0 8 14.55 4.36 


Table 10 (continued) Final exam marks. 


From this we scale the total down to 30. The mean resulted is thus 7.74, 
the median 7.23, minimum 0.05, maximum 18.27, and the standard deviation 
4.46. 
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661 
801 
802 
803 
805 
806 
807 
808 
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810 
811 
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4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 


819 
820 
821 
822 
823 
824 
825 
826 
827 
829 
830 
831 
832 
833 
835 
836 
837 


4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 


Score 
0.65 
14.02 
13.64 
11.24 
14.78 
18.27 
5.29 
5.45 
0.05 
0.27 
9.93 
13.2 
10.53 
3.27 
5.45 
8.18 
10.47 


Score 
7.64 
13.25 
13.04 
8.35 
13.8 
8.24 
11.29 
15.98 
5.89 
4.36 
2.18 
0.65 
2.29 
13.53 
4.36 


Table 11 Final examination marks and ranks, Business Mathematics, 2005-6, totalling 30 per 
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Final exam score distribution, Business Maths, 2005-6 
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Figure 20 Distribution of students’ final examination’s scores. 


Total exam score 


All two exams and three quizzes total to 80 per cent. Of this the average score 
is at 31.63, the median 32.03, minimum 11.02, maximum 51.81 and standard 
deviation 8.85. Then the scores and their ranks are shown in Table 12 and 
the distribution of scores in Figure 21. 
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814 
815 
816 
817 
818 
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ID 

819 
820 
821 
822 
823 
824 
825 
826 
827 
829 
830 
831 
832 
833 
835 
836 
837 
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Score 
24.22 
32.25 
41.47 
44,24 
35.95 
51.81 
20.92 
19.22 
24.49 
23.24 
42,23 
27.6 

41.19 
19.51 
31.09 
38.22 
33.74 


KR RRR RBBB RRR BREA 


872 


Table 12 Total exam score, Business Mathematics, 2005-6, 80 per cent of total 


Score 
30.17 
45.19 
41.77 
40.15 
13.8 

29.5 

38.36 
44.65 
31.16 
28.46 
27.35 
26.69 
32.76 
40.49 
36.3 
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Total exam score distribution, Business Maths, 2005-6 
T T T T T T T 


a 
T 
1 


aN 
T 
1 


io) 
T 
1 


number of students 


Figure 21 Distribution of total scores from all three quizzes and two 
exams. 


Homework 


For the marks for homeworks the mean is 8.08, median 8.00, minimum 7.00, 
maximum 9.00 and standard deviation 0.44. 
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Attendance 


The scores for attendance have as their average 9.33, median 9.5, minimum 
5.00, maximum 10.00 and standard deviation 1.02. Scores and their ranks are 
given in Table 13. 
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661 
801 
802 
803 
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819 
820 
821 
822 
823 
824 
825 
826 
827 
829 
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832 
833 
835 
836 
837 


Table 13 Score and rank of attendance. 
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ID Midterm Quiz1 Quiz2 Quiz3 Final AW & Atd 


4661 5.9 0 5.2 10 5.67 12 
4801 8.13 1.6 8.6 10 6.44 16 
4802 7.87 2.1 6.1 10 9.49 17.5 
4803 9.1 4.5 8.2 10 4.42 17.5 
4805 4.97 1.6 6 10 11.24 19 
4806 6.17 2.4 6 10 10.09 17.5 
4807 5) 1.6 7.9 10 10.42 16.5 
4808 4.13 1.4 6.9 2 2.62 18.5 
4809 5.07 1.5 3.4 2 4.64 17.5 
4810 5.4 2.4 8.2 10 10.09 18 
4811 6.6 4 7.2 10 3.71 17 
4813 7.1 4 8.2 2 5.45 15.5 
4814 9.4 1 8.2 10 9.98 17.5 
4815 10.8 2.4 8.2 10 7.09 18 
4816 8.43 - 6.2 10 9.39 17.5 
4817 3.63 2.3 4 2 4.53 18 
4818 8.23 6.3 8 10 11.02 18 
4819 5.27 1.9 4.2 10 3.27 18 
4820 3.37 1 4 2 0.65 17.5 
4821 4.3 1.8 4 - 4.8 17.5 
4822 4.27 1.6 7.9 10 13.69 16.5 
4823 10.03 1.6 7.7 10 6.98 16.5 
4824 9.87 1.6 8.2 10 16.25 17.5 
4825 6.7 1.6 5) - 8.78 14.5 
4826 10.77 4.2 8.2 10 6.44 18.5 
4827 10.77 3.9 8.2 10 3.82 18 
4829 10.07 2.2 8.2 10 6.82 18 
4830 3.77 1 4.2 10 1.09 17 


Table 14 All results compiled into a single table, Business Mathematics, 
2005-6 
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ID Midterm Quiz1 Quiz2 Quiz3 Final AW & Atd 


4831 4.63 4 8.2 2 12.98 17 
4832 4 1.5 6 10 3.27 17.5 
4833 4.33 1.4 4 10 9 17.5 
4835 4.27 1.3 7.8 5) 7.36 18 
4836 4.83 2.2 5.7 9 8.45 18 
4837 10.13 1 4 10 5.06 16.5 
4838 10.37 4 7.2 2 0.65 17 
4839 10.53 2 3.7 2 14.02 18 
4840 11.03 1.6 5.2 10 13.64 17 
4841 12.7 4.1 6.2 10 11.24 18 
4842 8.87 1.9 5.4 5) 14.78 18 
4843 12.53 2.8 8.2 10 18.27 19 
4844 3.83 2.8 - 9 5.29 16.5 
4845 5.77 1.8 4.2 2 5.45 18 
4846 6.23 1 7.2 10 0.05 17.5 
4847 10.77 1.2 7 4 0.27 19 
4848 12 2.1 8.2 10 9.93 19 
4849 = 1.7 7.7 5) 13.2 17.5 
4850 10.47 2 8.2 10 10.53 17.5 
4851 8.53 2.6 5.1 a 3.27 17.5 
4852 8.43 1.7 5.5 10 5.45 17.5 
4854 7.53 4.3 8.2 10 8.18 18 
4855 9.47 1 7.8 5) 10.47 17 
4856 8.33 4 8.2 2 7.64 17.5 
4857 9.33 4.4 8.2 10 13.25 18 
4858 10.03 1.7 7 10 13.04 19 
4859 9.6 4 8.2 10 8.35 18 
4860 = - - 3 13.8 12 
4861 8.77 1.6 6.9 4 8.24 19 
4862 7.87 1.6 7.6 10 11.29 15.5 
4863 7.57 2.9 8.2 10 15.98 18 
4865 8.77 2.4 5.1 9 5.89 17 
4866 8.1 1.8 4.2 10 4.36 18 
4867 7.27 2.7 5.2 10 2.18 18 
4868 7.13 1.6 7.3 10 0.65 17.5 
4870 8.77 4 7.7 10 2.29 18 
4871 10.37 1.9 4.7 10 13.53 18 
4872 9.73 4 8.2 10 4.36 18 


Table 14 (continued) All results compiled into a single table, Business 
Mathematics, 2005-6. 


The total score has the mean 49.04, median 49.53, minimum 25.8, maximum 


70.81 and standard deviation 9.31. The score and ranking for every student 
are shown in Table 15. 
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ID Score Rank ID Score Rank ID Score Rank ID Score 


661 38.77 55 4819 42.64 49 4838 41.22 54 4856 47.67 
801 50.77 30 4820 28.52 65 4839 50.25 33 4857 63.19 
802 53.06 25 4821 32.4 64 4840 58.47 11 4858 60.77 
803 53.72 24 4822 53.96 21 4841 62.24 5 4859 58.15 
805 52.8 27 4823 52.82 26 4842 53.95 22 4860 25.8 
806 52.16 28 4824 6342 2 4843 70.81 1 4861 48.5 
807 51.42 29 4825 36.58 60 4844 37.42 56 4862 53.86 
808 35.55 61 4826 58.1 13 4845 37.22 57 4863 62.65 
809 34.1 63 4827 54.68 18 4846 41.99 53 4865 48.16 
810 54.09 20 4829 55.28 17 4847 42.24 52 4866 46.46 
811 48.51 36 4830 37.06 58 4848 61.23 7 4867 45.35 
813 42.25 51 4831 48.82 34 4849 45.1 46 4868 44.19 
814 56.08 16 4832 42.27 50 4850 5869 9 4870 50.76 
815 56.49 14 4833 46.23 44 4851 37.01 59 4871 58.49 
816 47.48 41 4835 43.73 48 4852 48.59 35 4872 54.3 
817 34.46 62 4836 48.19 38 4854 56.22 15 

818 61.55 6 4837 47.2 42 4855 50.74 32 


Table 15 Total score and ranking, Business Mathematics, 2005-6. 
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ID Score 
661 38.77 
801 50.77 
802 53.06 
803 53.7 
805 52.8 
806 52.1 
807 51.4 
808 35.5 
809 34.1 
810 54. 
811 48.5 
813 42.2 
814 56: 
815 56.4 
816 52.7 
817 34.4 
818 
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829 
830 
831 
832 
833 
835 
836 
837 


Score 
42.64 
28.52 
36 

53.96 
52.82 
63.42 
40.65 
58.10 
54.68 
55.28 
37.06 
48.82 
42.27 
46.23 
43.73 
48.19 
47.20 


Score 
41,22 
50.25 
58.47 
62.24 
53.95 
70.81 
41.58 
37.22 
41.99 
42.24 
61.23 
56.38 
58.69 
41.12 
48.59 
56.22 
50.74 


Table 16 Score and rank after adjustment for absence from tests. 
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Score 
47.67 
63.19 
60.77 
58.15 
51.60 
48.50 
53.86 
62.65 
48.16 
46.46 
45.35 
44,19 
50.76 
58.49 
54.30 
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Total exam score distribution, Business Maths, 2005-6 


a 
T 


number of students 
w B 
T T 


70 80 90 100 


Figure 22 Distribution of the total score of tests, attendance and home- 


work. 


From the boundaries of each group portrayed in Figure 22 we arrive at the 


grading scheme in Table 17. 


Range 
[66.81, 100] 
[59.95, 66.81) 
[49.66, 59.95) 
[39.38, 49.66) 
[31.38, 39.38) 
[28, 31.38) 

[0, 28) 


Table 17 Grading scheme derived from Figure 22. 


And then the grades, together with total scores, are shown in Table 18. 
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ID 
4661 
4801 
4802 
4803 
4805 
4806 
4807 
4808 
4809 
4810 
4811 
4813 
4814 
4815 
4816 
4817 
4818 
4819 
4820 
4821 
4822 
4823 


Table 18 Total score and grade, Business Mathematics, 2005-6. 


Total 
43.07 
50.77 
53.06 
53.72 
52.80 
52.16 
51.42 
35.55 
34.10 
54.09 
48.51 
42.25 
56.08 
56.49 
52.75 
34.46 
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ID 
4824 
4825 
4826 
4827 
4829 
4830 
4831 
4832 
4833 
4835 
4836 
4837 
4838 
4839 
4840 
4841 
4842 
4843 
4844 
4845 
4846 
4847 


Total 
63.42 
40.65 
58.10 
54.68 
55.28 
37.06 
48.82 
42.27 
46.23 
43.73 
48.19 
47.20 
41.22 
50.25 
58.47 
62.24 
53.95 
70.81 
41.58 
37.22 
41.99 
42.24 


ID 
4848 
4849 
4850 
4851 
4852 
4854 
4855 
4856 
4857 
4858 
4859 
4860 
4861 
4862 
4863 
4865 
4866 
4867 
4868 
4870 
4871 
4872 


Total 
61.23 
56.38 
58.69 
41.12 
48.59 
56.22 
50.74 
47.67 
63.19 
60.77 
58.15 
51.60 
48.50 
53.86 
62.65 
48.16 
46.46 
45.35 
44.19 
50.76 
58.49 
54.30 


Grade 
Bt 
B 
B 
ct 
ct 


Here student 4860 who had been absent from too many exams was given a 
Ct instead of a B. Also, she said she would come to Quiz 3, but did not. 
That makes her marks fairly 51.60(0.9) = 46.44, that isa CT. 
Student 4661 had not come to Quiz 1. According to our procedure, therefore, 


the total of his score becomes 38.77/0.9 = 43.07, and a Ct at that. 
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Bangkok, 14° January, 2007 


Judge not, that ye be not judged. 
Matthew 7:1 
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Kinder is a german word that means ‘children’. So the name of the place 
should more correctly be written ‘Kinderscout’. But then again there are 
Kinder Downfall and Kinder Low. To do the same thing everywhere would 
probably result in something seemingly out of place sitting in an English 
context. Therefore the name is normally written ‘Kinder Scout’. 


Kit Tyabandha 
Bangkok, April 2006 


806 14 January, 2007 God’s Ayudhya’s Defence 


